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Abstract 

When light originating from a laser diode driven by non-fluctuating electrical cur- 
rents is incident on a photo-detector, the photo-current does not fluctuate much. Pre- 
cisely, this means that the variance of the number of photo-electrons counted over a 
large time interval is much smaller that the average number of photo-electrons. At non- 
zero Fourier frequency Vl the photo-current power spectrum is of the form Vl^ / (1 -f Vl?') 
and thus vanishes as ^ 0, a conclusion equivalent to the one given above. The 
purpose of this paper is to show that results such as the one just cited may be derived 
from a (semi-classical) theory in which neither the optical field nor the electron wave- 
function are quantized. We first observe that almost any medium may be described 
by a circuit and distinguish (possibly non-linear) conservative elements such as pure 
capacitances, and conductances that represent the atom-field coupHng. Configurations 
involving a single electron are considered. The theory rests on the non-relativistic ap- 
proximation, that is, c is set as infinite. Nyquist noise sources (in which the Planck 
term ftti;/2 is being restored) are associated with positive or negative conductances, 
and the law of average-energy conservation is enforced. We consider only first and 
second-order correlations of photo-electric currents in stationary regimes. 

Detailed semi-classical treatments of various quiet sources were listed in the first 
version of this paper. Only the general principles are presently considered, detailed 
applications being postponed. 
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1 Introduction 



"Comprendre", c'est comprendre autrement ("Comprehend" means compre- 
hend differently). 

In the present introduction we outhne our objective, main concepts, approximations 
employed, key results, and describe how the paper is organized. 

Scope of the paper. Laser noise impairs the operation of optical communication systems 
and the measurement of small displacements or small rotation rates with the help of optical 
interferometry. Even though laser light is far superior to thermal light, minute fluctuations 
restrict the ultimate performances. Signal-to-noise ratios, displacement sensitivities, and 
so on, depend mainly of the spectral densities, or correlations, of the photo-currents. It is 
therefore important to have at our disposal formulas enabling us to evaluate these quantities 
for configurations of practical interest, in a form as simple as possible. 

We are mostly concerned with basic concepts leaving out detailed practical calculations. 
Non-essential noise sources such as mechanical vibrations are ignored. Real lasers involve 
many secondary effects that are presently neglected for the sake of clarity. For example, 
because of the large size of the cavity in comparison with wavelength, lasers tend to oscillate 
on more than one mode. Even if the side-mode powers are much reduced with the help 
of distributed feed-backs or secondary cavities, small-power side modes may significantly 
influence laser-noise properties, particularly near the shot-noise level. Side-mode powers 
should probably be less than 40 dB below the main mode power to be insignificant. In 
the case of gas lasers, multiple levels, atomic collisions, thermal motions, and so on, may 
strongly influence noise properties, but these effects are neglected here. 

The main purpose of this paper is to show that, contrary to what most previous works 
imply, the properties of quiet lasers may be understood on the basis of a simple semi-classical 
theory, that is, a theory in which neither the optical field nor the electron wave- function are 
quantized. The electrons may be uncoupled to one another (dilute atom gases) or strongly 
coupled as is the case in semiconductors, through the Pauli exclusion principle. This theory 
(proposed by one of us in papers from 1986 on, and in book form in 1989 p|) is accurate 
and easy to apply, yet little known. The physical concepts are hopefully better explained in 
the present paper than in previous ones. Once the necessary assumptions have been agreed 
upon, laser noise formulas for various configurations follow from elementary Mathematics. 
In particular, operator algebra is not needed. Previous Semi-classical theories and Quantum 
Optics theories may be found in [2]. 

In this introductory part, the principles are presented but application to simple circuits, 
to the noise of lasers incorporating multilevel atoms or having spatially varying phase- 
amplitude coupling factors, the linewidth of inhomogeneously-broadened lasers, and the 
role of electrical feed-backs, is postponed. 
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Our interpretations of the basic mechanisms behind quiet-laser operatiorQI, the role of 
the Petermann K-factoJl and non-linear gain, on some respects differ from those found 
elsewhere in the literature. For simple laser models we observe however exact agreement 
between our results and those derived from Quantum Optics, whenever a comparison can 
be made. 

Theories found in most of the Optical Engineering literature rest on the concept that 
the classical oscillating field is supplemented by a random field. They involve parameters 
whose values are difficult to establish before-hand, because it is difficult to evaluate quan- 
tities such as the non-linear gain coefficient separately from the complete laser operation. 
Accurate models of real lasers (optical frequencies) are more difficult to establish than for 
masers (microwave frequencies) because of the small light wavelength. As a consequence, 
agreements reported between theory and experiments may not be conclusive as far as fun- 
damental concepts are concerned. Discrepancies between our semi-classical theory and 
previous semi-classical theories occur near the shot-noise level, and are therefore difficult 
to establish experimentally. Detailed comparison with observations may require specially 
constructed lasers. 



Main concepts of the theory and approximations. The important conclusion of 
the present discussion is the following. When an oscillator is driven by a constant power 
source and dissipates energy by amounts occurring at random times, the power delivered 
to the outside is nearly non-fluctuating. That is, the oscillator radiation is quiet. On many 
respects, lasers are akin to any oscillator, whether mechanical or electrical. We consider 
below the so-called "grand-mother pendulum" for the purpose of showing that the main 

^Prom our view-point a constant pump current entails a constant photo-current under ideal conditions. 
In a recent book [3J the basic mechanism behind quiet-laser operation is described as follows: "Although the 
noise generated in the external resistor is far below the shot-noise level, this does not mean that the carrier 
injection into the active region is regulated[...]. The carriers supplied by the external circuit are injected 
stochastically across the depletion layer before they reach the active region" (Presumably, by "stochastic" 
the authors means "Poisson-distributed"). The authors then introduce potential fluctuations to explain 
the observed quiet radiation. It may be, however, that these authors description is just another way of 
describing the same Physics as in the present paper. 

^It was observed early by E.I. Gordon [4] that in the linear regime laser line-widths are enhanced above 
those given by the well-known Schawlow-Townes formula for various circuits involving lumped elements or 
transmission lines, see O p. 120]. This linewidth-enhancement, which relates to the fact that gain and loss 
regions occur at different locations, may be described alternatively in terms of resonant complex potentials 
(or integrals of resonant complex fields). The line- width enhancement factor may be observed in strictly 
single-mode laser or maser oscillators. The K-factor effect discovered by Petermann [6] is of great practical 
importance for some laser diodes. It is not in our opinion fundamentally different from the Gordon effect 
just described. Further, contrary to a wide-spread belief, the K-factor cannot be applied directly to above- 
threshold lasers. The law of average-energy conservation tells us that quiet light should be observed with 
a quiet pump, leaving aside non-unity quantum efficiency, current leakage, and so on, irrespectively of the 
value of the K-factor value. An entirely different view point has been recently expressed [7]: "Because of 
the non-orthogonality of the laser eigenmodes, the noise from other modes is homodyned into the lasing 
mode, leading to an excess noise in the lasing mode, which is the Petermann excess noise" (our emphasis). 
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features of quiet-radiation oscillators may be understood on a purely classical basis, and are 
in fact, most common. This conclusion is based on elementary mechanical considerations 
that have nothing to do with Quantum Optics. (Note that in order to turn a "grand- 
mother pendulum" into a grand-mother dock one would need introduce a period-counting 
mechanism that would perturb the system as it is presently described). In similar laser 
oscillators the Planck constant h enters because the coupling of the oscillator to absorbers 
is effected with the help of electrons that may reside in only two states separated in energy by 
e = fkv. The mechanical model discussed below differs from laser models in that dissipation 
events remain Poisson-distributed. The regulation then originates from the fact that the 
energy lost at a "dissipation event" (to be defined later on) is proportional to the time 
interval that occurred from the previous event to the present one (this is called a "marked" 
Poisson process with marks depending on the Poisson-process realization). 

The basic element of a grand-mother pendulum is a mass M suspended at the end of a 
weightless bar of length L, in the earth gravitational field with acceleration g. As was first 
shown by Galileo, the oscillation period T = 27r-^ L/g does not depend on the oscillation 
amplitude as long as this amplitude remains small. Note also that the period does not 
depend either on the mass M. The parameters L, g, M (and later on m and 6z), and thus 
the period T, are fixed quantities in the present discussion. For simplicity, we assume that 
T = Is. If we label successive periods by A; = 1,2... we identify time t = kT = k. Let us 
call z{t) the pendulum mass height, with ^(1/4) = -z(3/4) = 0, the lowest mass position. 
We set z{0) = /iq, -^(1/2) = hi and z{l) = h2, see Fig. [H If the pendulum motion is not 
disturbed, the energy E = Mgh^ = Mghi = Mgh2, since the pendulum mass possesses no 
kinetic energy when its maximum height is being reached. 

Because of dissipation, and ignoring for the moment fast fluctuations, the pendulum 
energy decays in the course of time according to a law of the form E{t) = E{0) exp(— t/rp), 
where E{0) denotes the energy at t = and the constant Tp is called the pendulum life- 
time. The power Pd{t) dissipated in the medium surrounding the pendulum is viewed 
as a measurable quantity, perhaps through temperature increments. The lifetime Tp may 
therefore be measured, perhaps after averaging over a large number of identical pendula 
having the same initial heights. This life-time will be subsequently evaluated on the basis 
of a microscopic theory of damping. 

In order to obtain stationary oscillations some energy must be supplied to the pendulum. 
This is achieved with the help of a mass M suspended at the end of a cord (for simplicity 
this mass is assumed to be equal to the pendulum mass). An escapement mechanism allows 
the suspended mass to drop by a fixed height 5z at each swing of the pendulum, say at 
the end of every period, and then to increment the pendulum energy by 6E = Mg6z, so 
that the pendulum mass reaches a height exceeding the previous one by 5z. The pendulum 
period T = 1 being constant, this implies that the suspended mass delivers a constant power 
P = 6E, if we ignore the detailed process that may occur during a single period. Such a 
constant power supply may be called a "quiet pump" and is equivalent to the constant- 
current drive of lasers . Consider now both the power P supplied by the suspended weight, as 
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Figure 1: Variation of the spectral density of the power dissipated by a "grand-mother 
pendulum" as a function of the Fourier frequency il. The smooth curve is theoretical, 
while the slightly wavy curve results from a numerical simulation of the corresponding 
marked Poisson process. Inset, schematics of the pendulum of mass M. With probability 
pr <^ 1 the pendulum picks up a molecule of mass m <C M at 2; = and releases it at 
z = hi, causing damping. A suspended weight (also of mass M) drops by a length 5z at 
each swing of the pendulum, thereby incrementing its mass height by 5z. The escapement 
is represented symbolically. 
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discussed above, and the dissipation mechanism characterized by the empirical constant Tp. 
The equation of motion of the pendulum energy reads dE{t)/dt = P~P(i{t) = 5E—E(t)/Tp. 
In the steady-state (dE{t)/dt = 0) the average energy is therefore (E) = 5E Tp. 

We now model the damping mechanism at the microscopic level. We select a positive 
number pr < 1, and consider a random number x uniformly distributed between and 
1. If a; < 1 — no damping occurs. If ho denotes the pendulum mass height at the 
beginning of such an uneventful period, its height just after the period is /i2 = /io + as 
we discussed above. If, however, the random number x > 1 — pr, the pendulum picks up at 
t = 1/4, z = 0, a molecule of mass m, initially at rest, and releases it at t = 1/2, z = hi, 
again with zero velocity. If ho denotes the pendulum initial height, the maximum height 
hi reached is given by the law of conservation of energy by Mgho = [M + m,)ghi. Note 
that the mass increment does not affect the period (or the quarter of the period). At the 
end of the period we have h2 = hi + Sz = ho/{l + m/M) + Sz. The energy delivered to 
the molecule collection is then AE = mghi = mgho/{l + m/M). To summarize, given the 
height ho at the beginning of a period, the height at the end of the period is + with 
probability 1 — pr and /io/(l + m/M) + 6z with probability pr. 

Let the molecule-picking events occur at times ki, where i = 0,1,2.... We denote by 
hifi,hi^i,hi^2 the pendulum-mass heights as defined above, but relating to the ith-event. 
According to the previous description, molecule-picking events are Poisson distributed. 
At event i we attach a "mark" equal to the dissipated energy mghi^i. The previous 
description shows that jii = /Lti_i/(1 -|- m/M) + {mg/{l + m/M))6z{ki — ki-i) is a linear 
deterministic function of the previous mark and of the previous time interval ki — ki-i. 
Prom now on we assume that pr <^ l,m <^ M . The average /;.- value is then given by the 
average-power balance as {h) = 6z/pr, and the average dissipated energy for a molecule- 
picking period is e = mg {h) = mg6z/pr. The pendulum lifetime may now be evaluated as 
Tp = M/{pr m). As far as fluctuations are concerned, note that if a molecule-picking event 
occurs at t = and the next one at an anomalously-large time t » 1/pr, the pendulum 
mass height is much incremented, and the event at time t dissipates a large energy. This is 
how one can explain in a qualitative manner the mechanism behind dissipation regulation 

The energy dissipated from k = up to k = K may be written E^ = Ylk=o l^kS{k — ki), 
where 5{k — ki) = I'lik = hi and zero otherwise. The spectrum of Pd{k) = Ed{k) — Eii{k — 1) 
may then be obtained from a Pourier transform with respect to k. The result is compared 
to the analytical result given below. The dissipated power may be measured by allowing 
the molecules to fall back to the z = level and letting them dissipate their kinetic energy 
into a bolometer, a device that measures powers through temperature changes with a time 
constant of a few seconds. The masses m raising is analogous to the charging of an optical 
detector battery. Regulation is observed only when an averaging is performed over a large 
time duration. 

We have performed a numerical simulation of the marked Poisson process outhned above 
for M = 1kg, L = g/(^TxY ~ 0.25m, where g =9.81 ms^^, so that T = Is. Further, we 
set dz = 1/um, so that 6E = Mg6z ^ lO^^j, m = 0.001kg, pr = 0.01. We obtain for 
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the average pendulum mass height h = 6E/{pr m g) ^ 10cm, (E) = Mgh Ij and 
Tp = M / {pr m) 10^ s. The numerical result for the spectral density of the dissipated 
power fluctuations is shown in Fig. [TJ and compared to the theoretical result given below, 
in which we set Tp = lO^s, e = mgh « 10~^j, {Pd) = pr e = P = 6E 10/uW. 

An analytical formula for the spectral density of the dissipated power fluctuation APd{t) 
may be obtained when the fluctuations about the average value are small. On the basis of 
calculations analogous to those given later on for lasers, we obtain 



where denotes the Fourier frequency. In the present numerical application, we have 
e = 10~^j, Tp = lO^s, and {Pd) ~ 10/iW. As the figure shows, there is a very good agreement 
between the above analytical formula and the numerical simulation of the marked Poisson 
process considered. 

Let us consider now generators of electromagnetic waves. The only difference that 
exists between a microwave oscillator such as a reflex klystron and a laser relates to the 
different electronic responses to alternating fields. In a microwave tube the electron motion 
is usually not harmonic and its coupling to a single-frequency electromagnetic field may be 
understood accurately only through numerical calculations. In contradistinction, masers 
and lasers employ basically two-level molecules or atoms, and this results in simplified 
treatment^. The phenomena of stimulated emission and absorption are essentially the 
same for every oscillator. The noise properties are also similar. Let us quote the Nobel- 
prize winner W. E. Lamb, Jr. ^ p. 208]: "Whether a charge q moving with velocity v in 
an electrical field £ will gain or loose energy depends on the algebraic sign of the product 
evS [...]. If the charge is loosing energy, this is equivalent to stimulated emission. [...] In 
the domain of electronics, a triode vacuum-tube radio-frequency oscillator was developed 
by L. de Forest in 1912. This was in fact the first maser oscillator made by man". One may 
go one step further and assert that any sustained oscillator is a laser. The act of counting 
oscillations, however, requires specific arrangements. 

The lasers considered oscillate in a single electromagnetic mode in the steady state. 
Only stationarjifl fluctuations of the currents driving the active elements are allowed. The 

^For two-level atoms, upward electron jumps (stimulated absorption) and downward jumps (stimulated 
emission) may be treated symmetrically according to the time-dependent Schrodinger equation. Strictly 
speaking, the two-level approximation holds rigorously only for electrons immersed in a magnetic field, the 
lower energy state corresponding to the case where the electron magnetic moment points in the direction 
of the field and the higher energy state corresponding to the electron magnetic moment pointing in the 
opposite direction. In the case of atoms the electron energy is bounded from below but may extend to 
arbitrarily large values. The symmetry between stimulated emission and stimulated absorption therefore 
rests on the approximation that two levels only are important. In particular, the scattering states are 
ignored. The two-level approximation may cause apparent violation of oscillator-strength sum rules. These 
difficulties are un-consequential in the present theory. 

■^A fiuctuation is called "stationary" when correlations of all orders are independent of the initial time. 
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system elements are supposed not to depend explicitly on time. 
• Basic set-up. 

An optical set up involves three basic components. First a light source driven by an 
electrical current (called the pump). Second, an optical circuit involving slits, lenses, 
beam-splitters, resonators, and so on, which we view as being conservative, that is, free 
of loss or gain. Third, Ught detectors deUvering photo-currents. Light sources deliver 
optical power while light detectors absorb optical power. Ideally, the detector photo- 
currents could be employed to pump the light sources so that the complete system 
could operate in an autonomous manner. Such equilibrium configurations will be 
discussed in subsequent parts. A two-state electronic system, with energy-separation 
hujo, exhibits a positive conductance and absorbs power if the applied potential U is 
slightly lower than hujo/e, and exhibits a negative conductance and emits power if the 
applied potential slightly exceeds fiujo/e. In practical configurations there is therefore 
a slight irreversible loss of energy, which may be carried away by acoustical waves 
from the electronic configuration. This energy loss, in the meV range, should not 
be confused with the large irreversible loss of energy that, according to the Quantum 
Optics view-point, occurs when a photon is radiated away into vacuum from an atom. 
The latter is on the order of 1 eV. 

In many experiments, we only need to know time-averaged photo-currents. This 
information suffices for example to verify that light passing through an opaque plate 
pierced with two holes exhibits interference patterns. The experiment is performed 
by measuring the time-averaged photo-currents issued from an array of detectors 
located behind the plate. Other experiments involving the transmission of information 
through an optical fiber require that the fluctuations of the photo-current about its 
mean be known. A light beam carries information if it is modulated in amplitude 
or phase. Small modulations may be obtained from the present theory by ignoring 
the noise sources, but they are not discussed explicitly for the sake of brevity. The 
information to be transmitted is corrupted by natural fluctuations (sometimes referred 
to as "quantum noise"). 

We restrict ourselves to stationary non-relativistic configurations. That is, the free- 
space permeability is set equal to zero, or, equivalently, the speed of light in 
free space, c, is set at being infinite. These quantities therefore nowhere enter into 
the theory, and questions having to do with special relativity are irrelevant. More 
precisely, electron velocities v are much smaller than c and transition frequencies u 
are much smaller than mc^ /h, where m denotes the electron mass. We acknowledge 
that under these conditions some atomic properties are being overlooked. Relativistic 
effects are, for example 1) the apparent increase of the electron mass, 2) the value of 

This adjective is employed difFercntly in the expression "stationary states" where "stationary" means that 
the electron wave-function modulus is time independent. 
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the electron magnetic moment fj. = fJ^B derived from the Dirac equation, and 3) the 
spin-orbit energy splitting. This splitting, which results from the fact that, crudely 
speaking, atomic electrons perform circular motions at velocity v in nuclei electrical 
fields and thus "see" magnetic fields, is in fact small in hydrogen atoms, but becomes 
important for heavier atoms because v/cis not negligible. Quantum Electrodynamics 
enables Physicists to evaluate: 4) the (Lamb) energy splitting between 2Pi/2 and 
251/2 hydrogenic states, 5) the correction l+a/27r... to the electron magnetic moment, 
where the fine-structure constant a 1/137 is set equal to zero in the non-relativistic 
approximation, 6) the radiative decay of excited-state atoms (in that case, however, 
indirect approximate methods based on Statistical Mechanics or the Classical Maxwell 
Equations with retarded potentials may be employed), and 7) the Casimir force. 

On the other hand, in previous semi-classical theories the spontaneously emitted 
field is considered to be the fundamental source of noise. The classical optical field 
is supposed to be incremented by the field spontaneously emitted by upper-state 
atoms with a phase uniformly distributed between and 2tt (hence the randomness). 
Instead, we view noise as basically originating from stimulated electron jumps from 
one state to another, and spontaneous electronic decay is neglected for simplicity in 
the major part of the paper. 

• Non-fluctuating driving currents. 

We almost exclusively consider laser diodes driven by constant (non-fluctuating) elec- 
trical currents. Such currents may be obtained from a battery or a large charged 
capacitance and a cold series resistance. This conclusion follows from the Nyquist 
formula derived from Classical Statistical Mechanics that says that at material ab- 
solute temperature Tm = OK no fluctuations are involved in an equilibrium state. 
Detailed analysis shows that the Nyquist formula holds when a steady current flows 
through the resistance as long as the Ohm law remains applicable. Alternatively, we 
may generate non-fluctuating currents from space-charge-limited cathodic emission. 
It is now-a-day possible to inject in a device one electron at a time. A discrete realistic 
picture of a non-fluctuating current is accordingly the regular injection of electrons, 
say one every nano-second, if only small Fourier frequencies are being considered. 
This discrete picture is employed in numerical simulations. We may also employ a 
very large inductance with a current flowing through it. 

• Non-fluctuating radiation. 

It occurred as a surprised to the physics community when Golubev and other [9] 
proved theoretically in 1984 on the basis of the Quantum Optics laws that lasers driven 
by a quiet pump (e.g., a non-fluctuating current) deliver sub-Poissonian (or "quiet") 
photon streams. From our viewpoint, this observation would be better expressed by 
saying that when a laser is driven by a non-fluctuating current and the output light 
is incident on a photo-detector, the photo-current does not fluctuate much. In the 
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latter formulation the notion of laser-light statistics is being by-passed. The above 
prediction then may be viewed as a strictly classical result, resulting from the law of 
conservation of the average energy, as we discuss below. What is non-classical (i.e., 
quantum in nature) are the shot-noise fluctuations. This so-called "Schottky effect" 
has been observed long ago in vacuum tubes. This is perhaps for this historical reason 
that the Schottky effect is often referred to as being a "classical effect". But because 
it originates from the discreteness of the electric charge, it should be viewed instead 
as an intrinsically quantum effect. If one considers integration times large compared 
with the duration between successive photo-electrons, the discrete character of the 
electrical charge flow tends to be washed out, the theory becomes classical in nature, 
and accordingly a non-fluctuating photo-current is obtainecil. 

• Law of average- energy conservation. 

Let us explain in some detail how the law of average energy conservation is being 
employed. The electrical pump raises electrons initially in the absorbing state at rate 
Re{t) and thus supplies a power hu!eRe{t), where huje denotes the electrons transition 
energy. (Note that electronic states are often referred to as "atomic states". Because 
there exists now truly atomic lasers, the distinction is important if confusion is to 
be avoided). When these electrons decay back, the energy in the optical resonator is 
incremented. The concept of "light energy" is understood here only in a restricted 
sense. In order to determine the energy contained in a laser resonator at some time, 
say t = 0, one may cut-off the pump and measure the number of subsequent photo- 
detection events. It should be noted, however, that semiconductors (incorporated in 
particular in laser diodes) contain some energy of their own that cannot easily be 
separated out from the field energy. 

Even though the optical field is not quantized here, the word "photon" is employed 
occasionally as another name for the energy of loss-less resonators. Precisely, the 
resonator energy is written as {m-\-l/2)hL0o, where m is called the number of photons 
in the resonator, cOq the resonator (angular) frequency, and h the Planck constant 
(divided by 27r). Likewise, the word "photon rate" is another name for electromagnetic 
power divided by hwo- Conversely, the energy in the optical resonator is employed 
to raise detecting electrons initially in the absorbing state, with transition energy 
huja, to the upper state at a rate Rait). This power is delivered to the external load, 

'it is interesting to note that similar concepts (relating this time to the conservation of the average 
angular momentum rather than to the average energy) were recently advanced by C.S. Unnikrishnan [lOJ. 
That author shows that if a pair of electrons in the singlet state is emitted and their magnetic moments 
are detected at separate locations at angles differing by 6, the only correlation consistent with conservation 
of the average angular momentum is the quantum result -cos(^)/4, if the readings are normalized to ±1/2. 
What is strictly "quantum" is the discreteness of the electron spin. In the large-spin S limit the correlation 
— cos{9)S{S + l)/3 approaches the correlation — cos(S)S'^/3 evaluated from classical considerations. B. 
d'Espagnat, though challenging Unnikrishnan's interpretation of Bell's results, seems to agree with his 
factual conclusions 
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perhaps followed by an electronic amplifier. Ideally, we have huja ~ fuve ~ huJo, in 
which case the source-detector configuration may be viewed as being reversibl^. The 
law of conservation of energy then says that the integral from t = to i = T of the 
power difference [Re{t) — Ra{t))hiOo is equal to the system energy increment from 
t = to t = T, consisting of electronic and field energy increments. These, however, 
are finite. It follows that in the limit T ^ oo, we must have Re{t) ~ Ra{t)- More 
precisely, limT^oo {Re{t))j^ — {Ra{t))rp = 0, where the substricts T refer to average 
values taken over a time duration T. In the Fourier domain, this means that the input 
and output rate spectral densities must be the same in the limit Q ^ 0. 

In the case of laser diodes, the number n of electrons in the conduction band fluctuates 
as a consequence of the laser-diode dynamics. The Fermi-Dirac law then tells us that 
the potential applied to the current-driven diode fluctuates, and should be written 
as U{t). If J denotes the non- fluctuating pump current the input power JU{t) is 
no longer a constant. However, detailed calculations show that the fluctuations of 
U have a negligible effect on the energy balance, so that the previous argument still 
holds. 

• Circuit representation. 

The configurations investigated in this paper are described in terms of conservative 
elements such as capacitances and inductances, whose values may be obtained from 
separate classical measurements, as is done is conventional electronic^ll, and positive 
and negative conductances. 

• Average conductances. 

To define conductances we will consider first a single electron treated according to 
the Schrodinger equation. In atoms, electrons are submitted to the static field of 
nuclei. We suppose instead that the electron is located between parallel perfectly- 
conducting plates at the same potential, pierced with holes, with plates at a negative 
potential outside. This configuration provides a zero field, while nuclei fields are of 
the form 1/r^ at a distance r from the nucleus. The difference affects the details of 
the wave-functions, but not the principles. 

^This situation is analogous to that of the reversible heat engines discovered by Carnot in 1824. Re- 
versibility occurs when bodies are contacted only when their temperatures are nearly the same. As Carnot 
acknowledged, a non-zero power (non-zero heat flow) occurs only when there is some temperature difference 
between the contacted bodies. However, it is legitimate to consider the limit in which this temperature 
difference tends to zero. If this is the case, the mechanical energy delivered per cycle tends to a well deflned 
limiting value. Cycles are then very slow and the power generated is very small. 

'^The evaluation of a capacitance from its geometric dimensions is straightforward. If one insists in 
evaluating inductances from their geometric dimensions one needs suppose that they contain electrons. 
The latter have magnetic moments and the magnetic permeability ^ may be much larger than /io just 
above the Curie temperature Tc, being of the form (T — Tc)~^'^ ■ In that case a non-zero inductance is 
obtained even though is set equal to zero. 
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Let us sketch the way the Schrodinger equation is being appHed. The one-electron 
wave-function \l^{x) describes an ensemble of identically-prepared systems. Accord- 
ing to Born, \ilj{x.t)\'^ denotes the probability density of finding the electron at x if 
a position measurement is performed at time t, and \ipip,t)f denotes the probabil- 
ity density of finding the electron momentum as p if a momentum measurement is 
performed at time t, where ip{p) is essentially the Fourier transform of ip{x). The 
electrons are submitted, besides static electrical fields, to electrical fields oscillating 
at some optical frequency lu. We only consider the resonant case, and obtain the 
usual Rabi oscillations. 

The Quantum Mechanical (QM)-averaged induced current {i{t)) is proportional to 
the QM electron average momentum {p{t)) obtained by integrating p\'i/;{p,t)f . The 
QM-average current {i{t)) is the product of a sinusoidal variation at the optical fre- 
quency and a much more slowly-varying envelope. Because both v{t) and {i{t)) vary 
essentially at the optical frequency, the QM-average power {P{t)) = v{t) {i{t)) may 
be further averaged over an optical period, introducing a factor 1/2. This power, 
denoted simply as P{t), delivered by the alternating potential, serves to increment 
the electron energy and possibly the static potential source energy. Note that these 
energies have signs, so that the energy "delivered" may in fact be an energy received. 

We first consider configurations in which the electron interacts with the field only for 
a finite time r, perhaps because it is fiying between the plates at some fixed speed. If 
the electron enters in the interaction region in the lower state, it has some probability 
of being in the upper state when it leaves the interaction region. If this is the case, 
the electron energy may be delivered to some external collector. Such considerations 
lead to an expression for the average conductance G "seen" by the optical potential 
applied to the plates (ratio of the induced current and applied potential), which is 
positive in the case just described, but may be negative if the electron enters in the 
interaction region in the upper state. For small values of the transit time r compared 
with the Rabi period, the conductance G does not depend on the strength of the 
optical potential. That is, the system is linear. 

The fiying electron configuration just described is however not the one that we con- 
sider in the major part of this paper. We consider instead an electron present all the 
time in the interaction region, but ascribe to it a probability 272P22 of spontaneous 
transition from the upper to the lower level, and a probability 271P11 of spontaneous 
transition from the lower to the upper level, where pii,p22 are the probabilities that 
the electron be in the lower or upper states, respectively, and 71,72 are non-negative 
parameters. Supposing that 71 = 0,72 = 7, the power P{t) delivered by the optical 
potential serves to increment the electron energy as before, and to supply energy to 
the static potential. Given the 7-parameters, we may evaluate the conductance "seen" 
by the optical potential in the steady-state, that is, for t ^ 00. This conductance is 
independent of the optical potential when 7 3> ^r, where Qr denotes the so-called 
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Rabi frequency. To avoid a confusion, let us note that spontaneous decay is usually 
viewed an irreversible loss of energy from an atom electron. In the present configu- 
ration, energy is simply transferred from the optical potential to the static potential, 
or the converse, the electron playing an intermediate role. To treat conveniently the 
situation presently described one must introduce mixed-state density matrices. 

To conclude, the system considered is described by a circuit, consisting of intercon- 
nected (conservative) capacitances, inductances, and positive and negative conduc- 
tances. More generally, we may allow conductances to depend on parameters such 
as the optical frequency w, the number n of electrons in the conduction band for 
semiconductors or the number rig of atoms in the emitting (upper) state for atomic 
lasers, on the emitted photon rate R, or on the medium temperature T^, or strain in 
solids. 

• Nyquist noise at optical frequencies. 

As hinted above, optical set-ups are viewed as black boxes characterized by in-going 
and out-going photo-currents, whose statistical properties are either prescribed or 
sought for. Once a medium has been described by a circuit we are concerned with 
potentials and currents varying at, or near, some optical frequency Wq. These will 
be called "optical potentials", V{t), and "optical currents", I(t), respectively, to 
distinguish them from static potentials, U, and slowly- varying currents, j(t). We 
introduce optical potentials (or electric fields) and optical currents (or magnetic fields) 
for the sole purpose of ensuring that the photo-currents conserve the average energy 
in the sense explained above. 

One needs the quantum form given by Nyquist in his celebrated paper, supplemented 
by the term huJo/2 previously suggested by PlanciJl. The complete formula will be 
referred to as the "Nyquist-like" formula. An experimental verification of that formula 
at a temperature of 1.6K is illustrated for example at the beginning of Gardiner's book 
on Quantum Noise [12|- As a matter of fact, only the Planck term is important in 
the major part of this paper because we suppose that absorber atoms are all in the 
lower state (T=OK) while all the emitter atoms are in the higher state (complete 
population inversion). Various methods will be presented showing that, in that case, 
the induced-current-fluctuations spectral density is equal to HuJoG, where G denotes 
the absolute value of the conductance. Current-noise sources are independent of one 
another. The Nyquist current noises may be supposed to be normal, i.e., jointly 
gaussian distributeciHi. 

^Planck actually treated atoms rather than field resonances as quantized oscillators, so that the nature 
of his contribution on that subject remains unclear. 

^Within our linear or linearized approximations the noise currents are therefore also normal and cross- 
correlations of any order may be obtained from second-order cross-correlations. It follows that measurable 
noise currents are time reversible, and thus do not refiect the fact that the optical circuit elements are 
causal. 



14 



Dependence of G on frequency lj. 

In general, the conductance G depends on frequency, as recalled above. A well-known 
theorem says that in the linear regime the Nyquist formula is applicable to frequency- 
dependent conductances G{ui), as long as the temperature is uniform. 

Dependence of G on the number n of electrons. 

Both the average conductance G and the spectral density of the induced- current 
fluctuations are proportional to the number n of electrons as long as these electrons 
are not coupled directly to one another. In semiconductors the electrons are directly 
coupled to one another through the Pauli exclusion principle, and the conductance 
G depends non-linearly on the number n of electrons. At an electron temperature 
Tm = OK (roughly equal to the semi-conductor temperature) a gain proportional 
to m}/^ would be appropriate. However, because the fluctuations of n are small, a 
linearized conductance of the form G{n) = Go+ [clG/dn) (n — rio) may be employed. 
The spectral density of the Nyquist-like noise sources are supposed to be unaffected 
by the small deviations of n from Ug- 

Dependence of G on the emitted (or absorbed) rate R. 

In semiconductors the conductance G may depend significantly not only on the num- 
ber n of electrons in the conduction band, but also, explicitly, on the emitted power. 
This effect is called here "gain compression" (another name is "non-linear gain"). It 
plays a significant role in laser-diodes operation, increasing in particular the laser- 
diode relaxation-oscillation damping. The usual Nyquist-like formula remains appli- 
cable in special circumstances only. 

Dependence of G on temperature. 

The Fermi-Dirac statistics tells us how the electron energy distribution broadens at 

the electron temperature increases, an effect that tends to lower the conductance. 
The temperature depends in turn on the pump power, the material heat capacity and 
on thermal conductances. 

Linear and linearized regimes. 

Only two limiting cases will be considered, namely the linear regime and the linearized 
regime. In the linear regime optical potentials and currents are proportional to the 
fundamental noise sources. The response of linear systems to specified sources is 
straighforward, but dispersion effects need investigation. This regime is applicable 
to lasers below the so-called "threshold" driving current and, usually, to attenuators 
and amplifiers. 

The linearized regime is applicable to well-above-threshold lasers. In the linearized 
regime one first needs evaluate average optical potentials and currents ignoring the 
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noise sources. This is the so-called "steady state". Next, one supposes that the 
deviations of the optical potentials and currents from their average values, denoted 
by A, are proportional to the fundamental noise sources. The latter enter again when 
flowing powers are being evaluated, that is, current noise sources are not given for 
free, so to speak, but they do enter in the power balance. This is because we take this 
effect into consideration that our theory differs from previous semi-classical theories. 

The intermediate situation in which the system is neither linear nor can be linearized 
that may occur for closed-to-threshold lasers is not considered. As said above, we 
treat only the stationary regime found when a laser is driven by a constant current, 
possibly supplemented by stationary fluctuations, and no element is explicitly time- 
dependent, in which case photo-detection events are stationary as well. 

We assume that the atomic polarization may be adiabatically eliminated, so that our 
equations involve only the optical field, proportional to the optical potential V, and 
the numbers of electrons in various levels. The latter derive from rate equations that 
may sometimes be simplified further by neglecting electron-population time deriva- 
tives ("slaving principle"). Because spontaneous decay plays only a secondary role in 
our theory it is ignored for the sake of simplicity in this introductory part. 

• Potential fluctuations and correlations. 

In laser diodes employing semi-conducting materials, a constant-current drive J en- 
tails a static potential U across the diode that slightly exceeds Eg/e, where Eg denotes 
the semiconductor energy gap, because the bottom of the conduction band is filled 
up with n electrons, according to the Fermi-Dirac distribution. Likewise, there are 
n holes at the top of the valence bands. The rate equations that we shall introduce 
later on involve random fluctuations of n, and thus fluctuations Af/ of the potential 
U . This fluctuation is very small, yet measurable. One may also measure the correla- 
tion between AC/ and the detected current fluctuation AJ. This correlation may be 
defined in such a way that it is independent of any linear optical loss that may occur 
between the laser and the detector. From our view-point, AC/" is a small secondary 
effect that may, initially, be neglected. 

• Light spectrum. 

The light spectrum is a well defined quantity. To observe it, is suffices to insert 
between the laser and the photo-detector a narrow-band, cold and linear filter whose 
response is centered at some frequency uim- The average photo-current {J{t)) is 
proportional to the light spectral density y{ujm)- 

In the linear regime, the light spectrum may be evaluated from the modulus square 
of the system response to the Nyquist-like noise sources. Instead of using a narrow 
linewidth filter as said above, the light spectrum may be derived from the photo- 
current spectrum, the latter being the auto-convolution of the light spectrum. 
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In the linearized regime, the hght spectrum may be evaluated by first neglecting 
amplitude fluctuations and considering frequency fluctuation Au;(i), the latter being 
defined from the time derivative of the phase fluctuations of the optical wave incident 
on the photo-detector, which may be evaluated from the linearized-system response to 
the Nyquist-like noise sources. Experimentally, frequency noise may be converted to 
photo-current noise if the detector is preceded by frequency-selective optical circuits. 
A dual-detector arrangement is advisable. 

Key results. When light originating from a laser diode driven by non-fluctuating elec- 
trical currents is incident on a photo-detector, the photo-current does not fluctuate much, 
as we emphasized earlier. Precisely, this means that the variance of the number of photo- 
electrons counted over a large time interval is much smaller that the average number of 
photo-electrons. As we shall see, this is a consequence of the law of average energy conser- 
vation. Lasers having that property are called "quiet lasers". Viewed in another way, at 
high power, the photo-current reduced power spectrum (the adjective "reduced" meaning 
that the spectrum singularity at 17 = has been removed) is of the form 17^/(1-1-17^), where 
Q denotes the Fourier frequency, and thus the spectral density vanishes at 17 = 0. This con- 
clusion is equivalent to the one given earlier concerning the photo-count variance. We will 
say that light is sub-Poissonian when the spectral density of the photo-current is less than 
the average rate at small Fourier frequencies. Note that some authors say instead that light 
is "sub-Poissonian" when its normalized correlation (to be later defined) at zero time delay 
is less than unity. The two definitions are in general non- equivalent. The Quantum-Optics 
view-point is that quiet light is "non-classical", while, for the reasons explained earlier, we 
view quiet radiation as being entirely classical. 

The conclusion that for quiet lasers photo-current spectral densities vanish at zero 
Fourier frequency holds as long as the elements involved are conservative. Accordingly, 
the conclusion holds irrespectively of dispersion (related to the so-called "Petermann K- 
factor"), of the value of the phase-amplitude coupling factor (introduced independently in 
1967 by Haken and Lax and usually denoted by a), and of the amount of gain compres- 
sion (introduced by Chanin and alternatively called "non-linear gain"). These effects do 
affect however the photo-current spectral density at non-zero Fourier frequencies, the laser 
linewidth, and other laser properties. Note that conventional vacuum tubes with space- 
charge-limited cathode emission such as refiex klystrons should also emit quiet electromag- 
netic radiation. We do not know whether this has actually been observed, nor whether it 
can be observed in consideration of the klystron modest efficiency, and of thermal or fiicker 
noises. 

Organization of the paper. Besides the introduction and the conclusion, this paper 
consists of six sections. The first one gives an account of the most relevant results in Physics. 
The second one lists mathematical results relating to deterministic or random functions. 
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The third one is a discussion of the Circuit Theory. The fourth outhnes the Classical 
and Quantum equations of electronic motion. The fifth presents a more general theory 
of electron-field interaction, based on spontaneous electronic transitions. The sixth offers 
alternative methods of establishing that the spectral density of Nyquist-like noise sources 
associated with a conductance is proportional to the absolute value of the conductance. 

Aside from historical works, many citations relate to our own work (1986-2006). It 
is our intention to provide in later versions a more comprehensive list. Some important 
references, not cited here, may be traced back, however, from the more recent papers and 
books cited. 



18 



2 Physics 



According to the latin poet Lucretius, a follower of Democritus, there are no forbidden 
territories to knowledge: "...we must not only give a correct account of celestial matter, 
explaining in what way the wandering of the sun and moon occur and by what power things 
happen on earth. We must also take special care and employ keen reasoning to see where 
the soul and the nature of mind come from,...". And indeed, the three most fundamental 
questions: what is the origin of the world? what is life? what is mind? remain subjects of 
scientific examination. Needless to say, the present paper addresses much more restricted 
questions. 

We will first recall how Physics evolved from the early times to present, no attempt be- 
ing made to follow strictly the course of history. The theory of light or particle motion and 
the theory of heat followed independent paths for a long time. The Einstein contributions 
proved crucial to re-unite these two fields early in the 20th century. We may distinguish 
"pictures" based on our in-born or acquired concepts of space and time that may not an- 
swer all legitimate questions nor be accurate in every circumstances, and complete theories. 
Quantum theory is considered by most physicist as being accurate and complete, although 
some questions of interpretation remain hotly debated. We will consider in some detail 
the theory of waves and trajectories that are essential to understand the mechanisms be- 
hind vacuum-tube and laser operation. We also offer view-points concerning the Quantum 
Theory of Light. 

2.1 Early times 

Prom the time of emergence of the amphibians, earth, a highly heterogeneous stuff, is our 
living place. On it, we experience a variety of feelings. We feel the pull of gravity, breath 
air, get heat from the fire and the sun, and feed on plants growing on earth and water. Our 
experience, both as human beings and as physicists, is based on these living conditions. 
One may presTime that natural selection led human beings to an intuitive understanding of 
geometrical-physical-chemical quantities such as space, time, weight, warmth, flavor, and 
so on. At some point in the evolutionary process a degree of abstraction, made possible 
by an enlarged brain, facilitated our fight for survival. An example of abstract thinking 
is the association with space of the number 3, corresponding to the number of perceived 
dimensions. People "in the street" may however wish to distinguish the two horizontal-plane 
dimensions and the vertical dimension, considering that, for the latter, up and down are 
non- equivalent directions. It was not appreciated in the ancient times that the distinction 
between "up" and "down" is caused by the earth gravitational field, and that people living 
on the other-side of the earth have the same feelings as we do in their every-day life, even 
though, with respect to our own reference frame, they are "up-side-down". As we shall see, 
analogous considerations may apply to time, according to Boltzmann. 

Another naturally evolving concept is indeed the distinction between past and future 
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and physical causality: matter acts on matter only at a later time. The so-called "arrow of 
time" is a much debated subject. According to Boltzmann, in an infinite universe, there may 
be large-scale spontaneous fluctuations of the entropy (that one may crudely describe as 
expressing disorder). Past — > future would correspond to the direction of increasing entropy. 
There may be times where the entropy decreased instead of increasing. But the distinction 
is purely a matter of convention (in analogy with the "up and down" distinction mentioned 
above). This view point is consistent with the fact that the fundamental equations of 
Physics are (with the exception of the rarely occurring neutral-kaon decay) invariant under 
a change from t to —t. There has been objection to the Boltzmann view-point, however, 
and most recent authors would rather ascribe the time arrow to cosmic evolution, with the 
universe starting at the "big-bang" time in a state of very low entropy. 

In contrast with the rational view concerning causality, the magic way of thinking 
presupposes the existence of causal relationships between our desires, fears, or incantations, 
and facts. Now-a-days, magic thinking co-exists with rational thinking probably because 
it gives people sharing similar beliefs a sense of togetherness and helps a few individuals 
acquire authority and power. The consequences of irrationality are often too remote to be 
of concern to most. 

The control of fire by man some 500 000 years ago and drastic climatic changes that 
occurred, mainly in Europe, some 23 000 years ago, trigerred evolutionary events. Likewise, 
the practice of growing crops made possible a population explosion some 10 000 years ago, 
particularly in Egypt, and gave an incentive for measuring geometrical figures, precisely 
accounting for elapsed times, and measuring weights. Let us now consider more precisely 
what is meant by matter, space and heat. 

Empedocle (~500 BC) viewed the world as being made up of four elements, namely 
earth, water, air and fire. These elements remain a source of inspiration for poets and 
scientists alike, but they are not considered anymore as having a fundamental nature. 
Democritus (~400 BC) pictured reality as a collection of interacting particles that cannot be 
split ("a-toms"). Aristotle wrote in his Metaphysics VIIL "Democritus apparently assumes 
three differences in substances; for he says that the underlying body is one and the same in 
material, but differ in shape, position, and inter-contact". This picture may still be viewed 
as being basically accurate. 

The present work is not concerned with the cosmos per se. Yet, one cannot ignore that 
observations of the sky have been a source of inspiration in the past and remain very much 
so at present. Early observers distinguished stars from planets, the latter moving apparently 
with respect to the former. The ancient Greeks (Ptolemeus) conceived a complicated system 
of rotating spheres aimed at explaining the apparent motion of these celestial objects. 
Aristarque (310-230 BC), however, realized that the earth was rotating about itself and 
about the sun, the latter being considered to be located at the center of the universe. This 
heliocentric system was rediscovered by Copernic (1473-1543) and popularized by G. Bruno 
(burned at stake in Rome in 1600 for heresy). Next came the establishment of the three laws 
of planetary motion by Kepler, the dynamical explanation of these laws by Newton, which 
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involves a single universal constant, namely G. The deeper theory proposed by Einstein in 
1917 (General Relativity) appears to be in good agreement with observations. It involves 
two fundamental constants, namely G and c. 

When two bodies are in thermal contact they tend to reach the same temperature. 
Thus, two differently constructed thermometers may be calibrated one against the other 
by placing them in the same bath and comparing their readings. In the case of thermal 
contact the hotter body loses an amount of heat gained by the colder one but the converse 
never occurs. It may well be that the condition of heat-engine reversibility, discovered by 
Carnot in 1824, could have been made at a much earlier time and could have served as 
a basis for subsequent developments in Physics. The present attitude is rather that one 
should derive the laws of Thermodynamics from Classical or Quantum theories. It may 
be however that, to the contrary, the latter theories cannot be formulated unambiguously 
without the former. 

2.2 How physicists see the world now-a-day 

Beyond a qualitative understanding of the nature of heat, early observers were able to 
perform measurements of temperature and gas pressure with fair accuracy. Temperatures 
were measured through the expansion of gases at atmospheric pressure, linear interpolation 
being made between the freezing (0°C) and boiling (100°C) water temperatures. The con- 
cept of absolute zero of temperature emerged through the observation that extrapolated 
gas volumes would vanish at a negative temperature, now known to be -273.15°C=0 kelvin. 
The Classical Theory of Heat was established in the 18th and 19th centuries mainly by 
Black, Carnot and Boltzmann. The major contribution is due to Carnot (1824) who intro- 
duced the concept of heat-engine reversibility. The fact that hot bodies radiate power was 
known very early (some reptiles possess highly-sensitive thermal-radiation detectors). It is 
however only in the 19th century that the proportionality of the total radiated power to the 
fourth power of the absolute temperature was established. DifHculties relating to the theory 
of blackbody radiation led Einstein around 1905 to the conclusion that Classical Physics 
ought to be replaced by a more fundamental theory, namely the Quantum Theory, even 
though important conclusions may be reached without it. Another motivation for studying 
in some detail the theory of heat is that lasers are in some sense heat engines. They may 
be "pumped" by radiations originating from a hot body such as the sun. But, just as is the 
case for heat engines, a cold body is also required to absorb the radiation resulting from 
the de-excitation of the lower atomic levels. Lasers are able to convert heat into work in 
the form of radiation, but their efficiency is limited by the second law of thermodynamics. 
Output-power average values and fluctuations may be similar for lasers and heat engines. 

The grand picture we now have is that of a world 13 billions years old and 13 billions 
light-years across containing about 10^^ galaxies. Apparently, 80 % of matter is in a dark 
form, of unknown nature, that helped galaxy formation. Our own galaxy (milky way) 
contains about 10^^ stars and possesses at its center a spinning black hole with a mass 
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of 3.5 millions solar masses. Eight planets (mercury, venus, earth, mars, jupiter, saturne, 
uranus, neptune) are revolving around our star (sun). Penzias and Wilson discovered in 
1965 the cosmic background microwave radiation, which accurately follows the Planck law 
for a temperature of 2.73 kelvins. This cosmic black-body radiation is almost isotropic. Yet, 
minute changes of intensity according to the direction of observation have been measured, 
which provide precious information concerning the state of the universe some 300 000 years 
after the "big-bang". Numerous observations relating to ordinary stars such as the sun, 
neutron stars, quasars, black holes are particularly relevant to high-energy physics. It is 
expected that gravitational waves emitted for example by binary stars or collapsing stars 
will be discovered within the next ten years or so. Their detection may require sophisticated 
laser interferometers operating in space. In such interferometers, laser noise plays a crucial 
role. Reactors aim at creating on earth conditions similar to those occurring in the sun 
interior, i.e., temperatures of millions of kelvins, and to deliver energy, perhaps by the 
year 2050. An alternative technique employs powerful lasers shooting at a deuterium- 
tritium target. A reduction of the laser-beam wave- front fluctuations are essential in that 
application. For a review of the present views concerning the universe, see for example |13| . 

2.3 Epistemology 

Epistemology is the study of the origin, nature, methods and limits of knowledge. Un- 
doubtedly, Physics is an experimental science. Its purpose is to predict the outcome of 
observations, or at least average values of such observations over a large number of similar 
systems, from few principles using Mathematics as a language. Observations are required 
to set aside as much as possible human subjectivity. This is done by performing a large 
number of "blind" experiments, the same procedure being repeated again and again in 
independent laboratories. A physical theory should be "falsifiable", that is, one should be 
able to realize, or at least conceive, an experiment capable of disproving it. 

The average value (a) of a quantity a is calculated by summing aP{a), where P{a) is 
the probability density of a. It is apparently difficult to provide an unambiguous definition 
of the word "probability". Let us quote Dose [14] "There is a fundamental mistrust in 
probability theory among physicists. The need to extract as comprehensive information as 
possible from a given set of data is in many cases not as pressing as in other fields since active 
experiments can be repeated in principle until the obtained results satisfy preset precision 
requirements. [In other fields] the available data should be exploited with every conceivable 
care and effort". As data comes in our estimate of the probability P{a) improves, and 
eventually approaches an objective value, defined according to the frequentists view-point. 

In practice, most scientific progresses were accomplished with the help of intuitively- 
appealing pictures, describing how things happen in our familiar three-dimensional space 
and evolve in the course of time. These pictures are supposed to tell us how things are be- 
hind the scene, or to suggest calculations whose outcome may be compared to experimental 
results. Let us quote Kelvin: "I am never content until I have constructed a mechanical 
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model of the subject I am studying. If I succeed in making one, I understand; otherwise I 
do not". But many models, helpful at a time, are often discarded later on in favor of more 
abstract view-points. The Democritus picture of reality has been worked out in modern 
time by Bernouilli, Laplace and a few others. Given perfectly accurate observations made 
at some time, called "initial conditions", the theory is supposed to predict the outcome of 
future observations if the system observed is not perturbed meanwhile. (Poincare, how- 
ever, pointed out that for some systems, e.g., three or more interacting bodies in Celestial 
Mechanics, the error grows quickly in the course of time when the initial conditions are 
not known with perfect accuracy). The equations that describe ideal motions are time 
reversible, so that when the system is known with perfect accuracy at a time its state in 
the past as well as in the future is predictable. Predictions for earlier times (perhaps a 
misnomer) make sense if measurements were then made but not revealed to the physicist. 
What we have just described is sometimes referred to as the "Classical Paradigm". 

Reality is surely a concept of practical value. Anyone wishes to distinguish reality, as 
something having a degree of permanency, from illusions or dreams that are transitory in 
nature. On some matters, the opinions of a large number of people are sought, supposing 
that their agreement would prevent individual failures. In that sense, reality may exist 
independently of observers and be revealed by observations. But according to Bohr the 
purpose of Physics is not to discover what nature is, but to discover what we can say 
consistently about it. We stick to the Bohr view-point that observations relate only to 
complete set ups, including the preparation and measurement devices, the latter being 
considered classical. A specific measurement device is described in p5]. In effect, the 
object to be measured should be able to switch another object involving a large number of 
degrees of freedom from one metastable state to another. 

Advanced notions are not needed in this paper. It seems nonetheless that some under- 
standing of the Physics conceptual difficulties is useful. Seemingly reasonable pictures may 
fail to agree with observations in special circumstances. Indeed, consider a source and two 
measuring apparatuses, one located on the left of the source and the other on the right. 
Using Stapp |16j terminology, apparatuses may be set up to measure either size (large or 
small) or color (black or white) but not both at the same time. It is observecl^ that (l,b), 
(w,w) and (b,l) never occur, where large, small, white and black have been abbreviated 
by their first letters. The first term in these expressions correspond to the left-apparatus 
outcome and the second term to the right-apparatus outcome. In writing "large", for ex- 
ample, we of course imply that the apparatus has been set up to measure size, while in 
writing "black", for example, we imply that the apparatus has been set up to measure 
color. Let us now attempt to explain the above observations on the basis of the following 
picture: Assume that there exist four kind of particles, namely (Iw), (lb), (sw) and (sb). 
The source is supposed to shoot out one of these particles to the left and one to the right 

^"in reality, we are referring here to Quantum Mechanical (QM) predictions rather than to real obser- 
vations. There has been, however, so many experimental observations that agree with QM, that one may 
overlook the fact that observations have perhaps not been made for the system presently considered. 
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according to some probability law (there are all-together 16 probabilities summing up to 
unity, but only 4 of them will be considered). However, the fact that (l,b) never occurs 
implies that pr(lw,lb)=0. Here, "pr(lw,lb)=0" means that the source is not allowed to shoot 
out a particle of the kind Iw on the left and a particle of the kind lb on the right. Indeed, if 
it were allowed to do so, the left apparatus, set up to measure size, would sometimes give 
"1", while the right apparatus, set up to measure color, would sometimes give"b", contrary 
to observation. For the same reason, the source is not allowed to shoot out "lb" on the 
left and "lb" on the right, a condition that we write as pr(lb,lb)=0. Next we note that 
the observation that (w,w) never occurs implies that pr(lw,lw)=0. Finally, the observation 
that (b,l) never occurs implies that pr(lb,lw)=0. Accordingly, the probability that "1" be 
found on both sides is, considering the four possible combinations, pr(l,l) = pr(lb,lb) + 
pr(lw,lb) + pr(lb,lw) + pr(lw,lw) = 0, where the probabilities obtained above have been 
employed. Observations reveal, however, that if "1" is found on the left side, the probability 
that "1" be found also on the right side is equal to 0.065, i.e., is non-zero. It follows that the 
picture of a source shooting out two particles disagrees with observations. Of course, the 
non-zero probability quoted above (6.25 per cent) applies to elementary particles having 
only two attributes, each of them exhibiting only two possible values, not to macroscopic 
objects that may have other, measurable, attributes. It is frequently the case that an ef- 
fect deemed impossible according to Classical Mechanics, for example the transmission of 
a particle through (or above) a barrier of greater energy, is in fact observed (tunneling). 
This is because it is considered impossible, even in principle, to measure particle energies 
on top of the barrier. 

We are not concerned in the present paper with Physics in general but only with sta- 
tionary configurations, so that the epistemology of that part of Physics could perhaps be 
made more precise. The system is allowed to run in an autonomous manner, that is without 
any external action impressed upon it, and there is a continuous record of some quantities, 
particularly the times at which photo-electrons are emitted or absorbed. Systems on which 
we may act from the outside are not considered. Photo-electrons may be accelerated to such 
high energies by static fields that no ambiguity occurs concerning their occurrence times. 
The question asked to the physicist then resembles the one asked to people attempting to 
recover missing letters from impaired manuscripts: can you determine the missing letters 
from the known part of the text? In the present situation one would like to be able to tell 
whether an event occurred during some small time interval, given the rest of the record. Or 
at least give the probability that such an event occurs in the specified time interval. In other 
words, given a large collection of similar systems, on what fraction of them does an event 
occur? Instead of being given impaired records, we may be given information concerning 
the various components that constitute the system, such as lenses, semi-conductors, and 
so on, characterized by earlier, independent measurements. These measurements are deter- 
ministic in nature because they are performed in the classical high-field (yet, usually, linear) 
regime. In view of the observed uncertainty, spontaneous noise sources must obviously be 
introduced somewhere in the theory. W'e consider that the noise sources are located solely 
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at emitters and absorbers, viewed as being similar in nature. Pound [T7] described earlier 
lasers in terms of a Nyquist theorem extended to negative temperatures. 



2.4 Waves and trajectories 

Physics courses usually first describe how the motion of masses may be obtained from the 
Newtonian equations. But it might be preferable to let students get first familiarity with 
classical waves, for example by observing capillary waves on the surface of a mercury bath. 
Such waves are described by a real function of space and time that one may denote tpi^x, t) in 
one space dimension. One reason (to be explained in more detail subsequently) to consider 
waves as being of primary interest is that the law of refraction follows in a logical manner 
from the wave concept, but does not from the ray concept. Once wave concepts have been 
sufficiently clarified, the many-fold connections existing between waves on the one hand, 
and particles or light rays on the other hand, may be pointed out. Note that, historically, 
the motions of macroscopic bodies and light rays were established first (around 1600) and 
the properties of waves later on (around 1800 for light and 1900 for particles). Few precise 
results concerning waves seem to have been reported at the time of the ancient Greece. 
Yet, casual observation of the sea under gently blowing winds suffices to reveal important 
features. Had such observations been made, the course of discoveries in Science would 
perhaps have been quite different from what actually occurred. 

Waves at the surface of constant-depth seas propagate at constant speed u. In realistic 
conditions there is some dissipation and the wave amplitude may decrease but the wave 
speed remains essentially unchanged. This is a striking example of a physical object whose 
speed does not vary, no force being impressed upon it. The only condition required is 
that the medium parameters (the sea depth in the present situation) do not vary from one 
location to another. 

In the 1630s Galileo observed that macroscopic objects move at a constant speed when 
no force is exerted on them, in contradiction with the then-prevailing Aristotle teaching. 
A related finding by Galileo is the principle of special relativity: The laws of Physics 
established in some inertial laboratory are the same in another laboratory moving at a 
constant speed with respect to the first. In the year 1637 Descartes proposed the following 
interpretation for the refraction of light rays at the interface between two transparent media 
such as air and water. Descartes associates with a light ray a momentum that he calls 
"determination" having the direction of the ray and a modulus depending on the medium 
considered but not on direction. He observes that the rc-component of the momentum 
should not vary at the interface as a consequence of the uniformity of the system in that 
direction, justifying this assertion by a mechanical analogy, namely a ball traversing a thin 
sheet. The law of refraction asserting that 008(6*1)/ 003(^2), where the angles are defined 
with respect to the x-axis and the subscripts 1,2 refer to the two media, does not depend 



^^The law of refraction is most commonly written as nisin(ii) = n2sin(i2), where the n are refractive 



on the ray direction, follows 




Note that Descartes was only 
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concerned with trajectories in space, i.e., he was not interested in the motion of light pulses 
in time, so that questions sometimes raised as to whether light pulses propagate faster or 
slower in air or in water are not relevant to his discussion. 

No one at the time suggested that there may be a connection between particles or light 
rays on the one hand, and waves on the other hand. The wave properties of light were 
discovered by Grimaldi, reported in 1665, and explained by Huygens in 1678. The wave 
properties of particles were discovered much later by Davisson and Germer in 1927. In 
modern terms the Galileo, Descartes (and later Newton) concepts imply that particles and 
light rays obey ordinary differential equations. But without the wave concept the law of 
refraction for light or for particles relies on observation and intuition rather than logic. 

A wave packet has finite duration but includes many wave crests. A key concept is that 
of group velocity defined as the velocity of the peak of a wave packet, or short pulse. In 
particular, what is usually called the "velocity" of a (non-relativistic) body is the group 
velocity of its associated wave. But usually wave packets spread out in the course of time. 
In the non-linear regime though, wave packets, called solitons, may exhibit particle-like 
behavior in the sense they do not disperse. Bore-like solitary waves created by horse-drawn 
barges were first reported by Russell in 1844. 

Let us be more precise about waves. As said above, waves are very familiar to us, 
particularly gravity waves (not to be confused with the Einstein gravitational waves) on 
the sea generated by wind, or capillary waves generated on the surface of a lake by a 
falling stone. Simple reasoning and observations lead among other results to the law of 
refraction. Waves are defined by a real function tp{x,t) for one space coordinate x, and 
time t, obeying a partial differential equation. If the wave equation is unaffected by space 
and time translations we may set tp{x, t) = f{x — ut) for arbitrary speeds u. This results into 
an ordinary differential equation for the function f{x) which in general admits solutions. 
Let us begin our discussion with monochromatic (single-frequency) waves propagating in 
the X direction in a conservative linear and space-time invariant medium. The wavelength 
A is the distance between adjacent crests at a given time. We define the wave number 
k = 2it/X. The wave-period T is the time it takes a crest to come back, at a given location. 
We define the frequency uj = 2'k/T. It follows from the above definitions that the velocity 
of a crest, called the phase velocity, is u = oj/k. Such waves propagate at constant speed 
without any action being exerted on them. For linear waves there is a definite relationship 
between oj and k independent of the wave amplitude, called the dispersion equation. For 
gravity waves in deep inviscid (non- viscous) waters we have, for example, uj = ^/gk, where 
g « 9.81m/s^ is the earth acceleration. When the water depth h is not large compared 
with wavelength (shallow water), the dispersion relation involves h as a parameter [HI see 
ref. 8]. 

The above considerations may be related to mechanical effects. Indeed, if a wave carry- 
ing a power P is fully absorbed, the absorber is submitted to a force F satisfying the relation 

indices and the i angles are defined with respect to the normal to the interface. 
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P/uj = F/k. This ratio, called "wave action", depends on the nature of the wave but does 
not vary if some parameter is changed smoothly, either in space or in time. For a wave 
of finite duration r, the energy collected by the absorber is £^ = Pt and the momentum 
received (product of its mass and velocity) is p = Ft. 

If the water depth h is changed at some time t = from, say, Im to 2m, it is observed 
that k is unchanged as a consequence of the wave continuity. But invariance of k implies a 
frequency change since the dispersion equation depends on h. In that case, the wave speed 
changes at time t. Conversely, If the water depth h changes at some location x = from, 
say, Im to 2m, it is observed that u is unchanged as a consequence of the wave continuity. 
But invariance of uj implies a wave number change since the dispersion equation depends 
on h. In that case the wave speed changes at a; = 0. 

Consider now a monochromatic wave (fixed frequency cv) propagating in two dimensions 
with coordinates x, y. The direction of propagation is defined as being perpendicular 
to the crests and the wavelength A = 2Tr/k is defined as the distance between adjacent 
crests at a given time. But one may also define a wavelength A^; in the direction x as 
the distance between adjacent crests in the x-direction at a given time. Let the wave 
be incident obliquely on the interface between two media, the x-axis. For gravity waves 
the two media may correspond for example to h{y) = Im, y > and h{y) = 2m, y < 0. 
Because of the continuity of the wave, A^; is the same in the two media. If we further 
assume that the propagation is isotropic, that is, that k does not depend on the direction 
of propagation of the wave in the x,y plane, the law of refraction follows, namely that 
kx = ki cos(0i) = k2 cos(^2), where the subscripts 1,2 refer to y > and y < respectively, 
and the angles 6 are defined with respect to the interface, that is, to the x-axis. The law 
of refraction therefore follows from wave continuity and isotropy alone. 

Questions relating to the velocity of light pulses are important for the transmission of 
information. A wave-packet containing many wave crests moves at the so-called "group 
velocity" v = duj/dk, which often differs much from the phase velocity u defined above. 
Considering only two waves at frequency lo and u) + du), the relation v = duj/dk may be 
visualized as a kind of Moire effect. Wave crests move inside the packet, being generated 
at one end of the packet and dying off at the other end. For waveguides we have uv = 
c^, V < c, u > c. For matter waves associated with a particle the group velocity v coincides 
with tlie particle velocity. Since the energy E = p^ /{2m) and p = mv, a previous relation 
reads p^ / (2muj) = p/k. It follows that u = u/k = p/2m = v/2. For gravity waves the 
dispersion relation gives instead u = 2v. A general result applicable to loss-less waves is 
that the group velocity v is the ratio of the transmitted power P and the energy stored per 
unit length. It never exceeds the speed of light c in free space. 

Wave solutions of the form i/j{x,t) = tp{x — ut), where tp{x) is some given function and 
u a constant, exist also for non-linear wave equations. When the tpix) function is localized 
in X, the invariant wave-form is called a solitary wave. In some cases, solitary waves exhibit 
transformations akin to those of particles when two waves collide and are called "solitons" 
in the sense that the soliton integrity is preserved. 
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As said before, most continuous media may be modeled by discrete circuits. For exam- 
ple, a transmission line may be modeled by series inductances and parallel capacitances. 
Free space may be modeled by electrical rings in which electrical charges move freely and 
magnetic rings in which (hypothetical) magnetic charges would move freely. If each elec- 
trical ring is interlaced with four magnetic rings and conversely, the Maxwell equations in 
free space obtain in the small-period limit. 

Under confinement along the x-direction, waves at some fixed frequency w may be 
viewed as superpositions of "transverse modes". For a transverse mode the wave-function 
factorizes into the product of a transverse function tpi^x; co) and a function of the form 
exp{ik{uj)z — iut). Another connection between waves and rays rests on the representation 
of transverse modes by ray manifolds. These are not however independent rays. A phase 
condition is imposed on them that leads to approximate expressions of ip{x;Lo) and k{uj). 
Note the analogy with Quantum-Mechanics stationary states, z and t being interchanged. 

Thus the wave-particle connection is many fold. First the medium in which the wave 
propagates may be approximated by a discrete sequence of elements, for example a periodic 
sequence of springs and masses for acoustical waves and electrical inductance-capacitance 
circuits for electromagnetic waves, with a period allowed to tend to zero at the end of the 
calculations. One motivation for introducing this discreteness is that computer simulations 
require it anyway. A more subtle one is that some divergences may be removed in that 
way. We have mentioned above capillary waves on a mercury bath. They may be treated 
by considering the forces binding together the mercury molecules and their inertia, ending 
up with equations of Fluid Mechanics. Like-wise, acoustical waves in air may be described 
through the collision of molecules in some limit (isothermal or adiabatic). Second, wave 
modes may be described approximately (WKB approximation) by ray manifolds. Third, 
one may consider the behavior of wave packets in the high-frequency limit and liken these 
wave packets average trajectories to those of macroscopic bodies. 

We have described above the motion of light and particles in terms of waves. Semi- 
classical theories such as the one employed in the present paper rest indeed on such wave 
concepts, namely Quantum Mechanics for describing electrons, and Circuit Theory for 
describing the relationship between potentials and currents. The speed of light in free- 
space is irrelevant in that theory. Quantization then only means that electrons are identical 
point particles. When particles such as electrons are electrically charged they may be 
accelerated to arbitrarily large energies by static electrical potentials. Being then in the 
classical domain there is no ambiguity concerning their arrival time. Uncharged point 
particles such as neutrons could conceivably be accelerated similarly by gravitational fields, 
even though this may turn out to be difficult in practice. 

2.5 Atoms and elements 

Around 1927 it was discovered theoretically by de Broglie and subsequently verified exper- 
imentally that a wave of wave- number k = mv/h should be associated with electrons of 
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velocity v. An approximate solution for the motion of an electron following a closed classical 
path in the neighborhood of a positively charged nucleus thus amounts to prescribe that an 
integral number n of wavelengths 2Tr/k fits along the closed classical path. These discrete 
solutions are called "stationary states" and n is essentially the principal quantum number. 
According to the Pauli principle, at most two electrons (with spin ±h/2) may be ascribed 
to each of these states. At Tm = OK and without excitation by other particles, only the 
lowest-energy states are filled with electrons. Different elements (H, He, Li...) differ by the 
number Z of protons in their nuclei. 

In the next paragraph we recall how the chemical and electronic properties of the various 
elements found in nature follow from the above principle, and describe what happens when 
atoms get closer and closer to one another to form crystals. Then we recall the basic 
properties of semi-conductors. 

2.6 Electron states 

We summarize below the most basic concepts concerning elements found in nature and 
their electron states. The simplest element is the hydrogen atom consisting of a proton 
with an electrical charge e and a mass much larger than the electron mass m, so that 
for most purposes the proton may be considered as being fixed in space. This proton 
attracts one electron of charge — e so that the assembly is neutral. According to Classical 
Mechanics the electron may circle around the proton at a distance r with a velocity v 
such that the centrifugal force be balanced by the attraction from the proton, namely 
m,v^/r = e^/ (47reor^) . Prom this view-point any distance r may occur, the velocity v being 
appropriately chosen. According to Quantum Theory a wave-length 2Trh/mv is associated 
with electrons moving at velocity v. The resonance condition is that an integral number n 
of wavelengths fits within the electron path perimeter 27rr. According to this model, due 
to Bohr, there is only a discrete sequence of allowed electron energies, corresponding to 
n = 1,2.... The more exact theory due to Schrodinger leads to symmetrical ground states, 
called s-states, and anti-symmetrical 3-times degenerate first-excited states, called p-states. 

The elements found in nature (roughly 100) were classified by Mendeleiev in 1869 on 
empirical grounds. Helium nuclei consist of two protons, lithium nuclei of three protons, 
and so on, with an equal number of electrons, so that atoms are electrically neutral. (There 
may be various numbers of neutrons bound to the protons, which depart from the number 
of protons by a few units, corresponding to different isotopes, some of them being unstable. 
Neutrons are not considered in the present discussion). Most elements have an outer layer 
consisting of a number of electrons going from 1 (e.g., sodium) to 8 (e.g., neon). Particularly 
important are 3-5 crystals, such as gallium-arsenide. 
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2.7 Semi-conductors 



Our purpose here is to give readers unfamiliar with solid-state physics an overview of 
the most important phenomena. For silicon, the number of outer electrons is 4. Two 
silicon atoms (or more) may bind to one another by exchanging electrons of opposite spins 
(covalent binding). When two atoms are approaching one another, their electronic states get 
perturbed. As it happens, the isolated-atom electron s-state acquires an energy greater than 
the isolated-atom electron p-states. For a large number N of atoms, the atomic separation 
a sets up at a value that minimizes the total energy. The original s-states then split into 
N states that are so-closely spaced in energy that they form an almost continuous band 
of states called the conduction band. The original 3-fold degenerate p-states split into 
states that are so-closely spaced in energy that they form three almost continuous band of 
states called the valence bands. Because the degeneracy is lifted these three bands should 
be distinguished. They are called respectively the heavy-hole band, the light-hole band and 
the split-off band. For our purposes, only the heavy-hole band needs be considered. 

The separation in energy between the bottom of the conduction band and the top of 
the valence band is called the band gap Eg, often expressed in electron-volts. At r=OK, 
the lower-energy valence band is filled with electrons while the higher-energy conduction 
band is empty. At that temperature the electrons are unable to respond to an external 
field because no state is available to them (except perhaps at extremely-high fields). If an 
electron is introduced in the conduction band by some means it moves in response to an 
electrical field with an apparent mass nic smaller than the free-space mass m. If, on the 
other hand, an electron is removed from the valence band one says that a "hole" has been 
introduced. This hole is ascribed a positive charge e and a mass usually larger than m. 

When two materials having different band gaps are contacted the band gap centers 
align approximately, and potential steps occur both in the conduction and valence bands. 
In the case of a double-hetero-junction the lower-band-gap material is sandwich between 
two higher-band-gap materials. The potential steps tend to confine both free electrons and 
free holes in the central low-band-gap material (e.g., GaAs). Being confined in the same 
volume electrons and holes easily interact. 

As the band-gap decreases electrons may undergo virtual transfers from one band to the 
other more easily. As a consequence the material is more easily polarized by external (static 
or optical) electrical fields. In other words, the material permittivity e(a;) increases as the 
band gap decreases. This is why the permittivity (or refractive index) of the low-band- 
gap gallium-arsenide is significantly larger than the permittivity (or refractive index) of 
the large-band-gap aluminum arsenide. When a small-band-gap semi-conductor (GaAs) is 
sandwiched between two higher-band-gap semiconductors (AlAs), the higher-index material 
may guide optical waves. This fact is important for the guidance of optical waves in 
laser diodes employing double-lietero-junctions. An happy circumstance is therefore that 
electrons, holes, and light, may all get confined in the central part of the double-hetero- 
junctions considered. 
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Gallium possesses 3 electrons in the outer shell and arsenide possesses 5 electrons. Equal 
numbers of these atoms may associate to form a crystal of gallium-arsenide (Ga-As), a ma- 
terial particularly important in Opto-Electronics. The reason for this importance is that, 
unlike silicon, this is a "direct band-gap" material. In direct band-gap materials the mini- 
mum of the conduction-band energy and the valence-band maximum energy correspond to 
the same electron momentum. Accordingly, electrons lying at the bottom of the conduction 
band may get easily transferred to the top of the valence band, and conversely, the law of 
momentum conservation being then fulfilled (photon momentum is negligible). In such a 
process, an energy Eg is absorbed by light through stimulated or spontaneous emission 
processes. Unfortunately, this energy may also be absorbed by another electron (Auger 
effect) that subsequently cascades down, its energy being converted into heat. 

Finally, one should say a word about doping, considering as an example a silicon crystal. 
When a small number of silicon atoms are replaced by arsenic atoms, these atoms, referred 
to as "impurities", easily deliver an electron (n-doping). Conversely, when a small number 
of silicon atoms are replaced by gallium atoms these atoms easily capture electrons (p- 
doping). A p-n diode consists of two contacting semi-conductors, one with p-doping and 
one with n-doping. Electron currents may be injected into p-n diodes, and in particular 
into double-hetero-junctions. This is the current referred to in this paper as the laser-diode 
driving current or the photo-current of a quantum photo-detector. It is denoted by J. 

The above discussion hopefully provides the essential concepts that one needs to get 
some understanding of the electrical behavior of laser diodes. Note that we denote by z the 
coordinate along which the optical wave propagates (junction plane) and by x the direction 
perpendicular to the semiconductor layers. Guidance along the transverse y direction is 
also considered. 

2.8 Detectors and sources 

In subsequent sections we discuss sources of electromagnetic radiations and ways of detect- 
ing them. It is appropriate to consider first detectors because there exist natural sources of 
radiation such as the sun, and the difficulty was initially to detect such radiations rather 
than to generate them. Detectors convert high-frequency radiation into slowly varying 
currents. The mode of operation of some detectors, called "classical detectors", may be 
explained on the basis of the Classical Equations of Electron Motion. For others, called 
"quantum detectors", the Quantum Theory of Electron Motion is required. Conversely, 
sources convert slowly varying currents into high-frequency radiation. The mode of op- 
eration of some sources, called "classical sources", may be explained on the basis of the 
Classical Equations of Electron Motion. For others, called "quantum sources", the Quan- 
tum Theory of Electron Motion is required. 

Historically, the first light detector was of course the human eye. Modern detectors were 
first vacuum tubes operating with a low-work function cathode and accelerating potentials. 
There exist now quantum detectors whose mode of operation is based on the phenomenon 
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of stimulated absorption. An early man-made generator of high-frequency radiation is a 
vacuum-tube called the "reflex klystron". The main light sources are hot bodies and lasers. 

2.9 Classical detectors 

Classical detectors are diodes that exhibit non-linear current-potential characteristics. If 
a sinusoidal potential is applied to the diode the current then exhibits a non-zero average 
value, which is a measure of the applied sinusoidal-potential amplitude. 

Let us recall the basic mode of operation of conventional electronic diodes, photo- 
detectors and photo-multipliers. Conventional electronic diodes are made up of two parallel 
plates (labeled in what follows the lower and upper plates) separated by a distance d in 
vacuum. The lower plate, called "anode" is at zero potential by convention, and the upper 
plate, called "cathode", is raised at the potential —U with [/ > 0. Suppose that at time 
t = an electron is freed from the upper plate and attracted by the anodJ^. Considering 
only absolute values, the electron momentum increases linearly with time t according to 
the law p{t) = elJt/d, where — e denotes the electron charge, until it reaches the anode at 
time r = dy^2m/eU, where m denotes the electron mass. The electron kinetic energy is 
then converted into heat. In the following, r is neglected, that is, it is set equal to zero. 
Fig. [2] illustrates in a), the photo-current, represented as a function of time. Because the 
output circuit capacitance is taken into account each electron arrival corresponds to an 
exponentially-decaying pulse of the form exp(— f/rc). In b), photo-current spectrum for 
the case where the output circuit is a resonating circuit tuned at some Fourier frequency 
Qo = I /Vic. 

In the case of photo-multipliers, the electron kinetic energy, instead of being dissipated 
into heat, is employed, at least in part, to free two or more electrons from the anode 
(secondary-emission effect). The latter are accelerated by a third plate, and so on, so 
that each electron freed from the cathode by light gets converted into an electron bunch 
containing n electrons, for example, n = 10^. The situation is the same as if the absolute 
value e of the electron charge had been multiplied by n. The purpose of photo-multipliers 
is not to improve the signal-to-noise ratio, which may only degrade. It is to raise the signal 
to such high levels that the thermal noise of subsequent electronic amplifiers is rendered 
negligible. 

^^To achieve this, the cathode "work function" energy must be overcome by heat (thermo-ionic emission), 
high electric fields (field emission), electrons (secondary emission), or light (photo-electric emission). Elec- 
trons may be freed by thermal motion provided k^Tm be of the order of the metal work-function. If nickel 
is coated with barium oxide, a temperature of 1000 kelvin may suffice. Field emission occurs with kilo- volt 
potentials if the cathode has the shape of a needle. Electrons may be freed by light provided huo exceeds 
the metal work-function, where h denotes the Planck constant and a; the light frequency. Visible light for 
example is adequate when the cathode is coated with cesium. The non-zero initial electron velocities are 
presently neglected, that is, the initial electron momentum p(0) — 0. Electrons in a metal are bound to it 
because they are attracted by their image charge. They may escape, though, because of a tunneling effect 
whose understanding requires Quantum Mechanics. But once the electron is sufficiently far away from the 
cathode, the Classical Equations of Electron Motion are appropriate. 
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(a) (b) 

Figure 2: In a), on top the photo-current is represented as a function of time, the output 
circuit capacitance (lower part) being taken into account. In b) the photo-current spectrum 
(on top) for the case where the output circuit (lower part) is a resonating circuit tuned at 
some Fourier frequency f^o- 

In temperature-limited thermo-ionic diodes the electronic density is so low that the 
Coulomb interaction between electrons may be neglected. The electrons are emitted inde- 
pendently of one-another and their emission times are Poisson distributed (see the math- 
ematical section). In that case the diode current fluctuation Aj(f) obeys the so-called 
shot-noise formula, with a (double-sided) spectral density equal to e{j{t)). But when the 
electron flow is space-charge limited the current is sub-Poissonian, i.e., the spectral density 
is much smaller than the one just given. This effect was discovered in 1940 by Thompson, 
North and Harris [T9] . 

The detectors considered above have been idealized for the sake of simplicity neglecting, 
e.g., dark currents and thermal noise. As said before, the current flowing out of photo- 
detectors may be viewed as a sum over the positive integer k of delta-functions of the form 
e5{t — tfc), where the are occurrence times. If this current is transmitted through a 
low-pass filter such as the one shown in Fig.[2l individual pulses may overlap, however, and 
not be distinguishable any more from one another. The current fluctuation then resembles 
gaussian noise irrespectively of the event-times statistics. 

2.10 Quantum detectors 

Quantum photo-detectors (sometimes referred to as "narrow-band" photo-detectors), in- 
volve two electron energy levels, coupled to continua, and operate through the process of 
stimulated absorption. Ideally, the device is reversible in the sense that the electrical en- 
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ergy may be converted back into light energy through the process of stimulated emission. 
In contradistinction, conventional photo-diodes necessarily dissipate energy in the form of 
heat. Of particular interest is the visible-light-photon-counter, see [H p. 181], which has a 
quantum efficiency of 0.88, a gain of 30 000, a time response of 2ns, but unfortunately a 
rather large dark count of 20 000 counts per second. 



2.11 Classical sources 

The first high-frequency oscillator was probably a triode, with a feed-back mechanism from 
the anode to the grid controlling the current fiow. We are particularly interested in the 
refiex klystron discovered by the Varian brothers in 1937, which may deliver electromagnetic 
radiation up to a frequency of about 10 GHz, and employs space-charge-limited cathodic 
emission (see e.g. |20j). As was recalled above, Thompson and others discovered in 1940 
that the current emitted by space-charge-limited cathodes is strongly sub-Poissoniar^l. 
Refiex klystrons were mostly employed as low-noise local oscillators in radar heterodyne 
receivers until they were superseded by solid-state devices. 

Refiex klystrons involve two anodes made up of grids presumed to be transparent to 
electrons and separated by a distance d. Two cathodes are located just outside the anodes. 
The lower one emits electrons, while the upper one plays the role of a refiectoi0 so that 
the electron motion as a function of time is a zig-zag path. The two anodes are part of a 
resonator. When the device oscillates an alternating potential v{t) appears between the two 
anodes. If the electron emission time is appropriate, the electron looses its energy giving it 
up to the oscillating potential through an induced current. But since the electron emission 
times are uniformly distributed along the time axis, the net interaction with the field 
vanishes. Accordingly, initially, the electrons do not deliver any energy to the oscillating 
potential. It is as a result of the field action on the electron trajectories that non-zero 
energy exchanges between the field and the electrons may occur. This effect is called 
"bunching". Once the electron has lost most of its energy it gets captured by the anodes, 
and instantaneously jumps from the anode to the cathode through the static potential^. 
The static-potential energy then gets reduced by ell. In some sense, the electron plays 
an intermediate role. Indeed, the net effect of a complete electronic cycle is that, for each 



^^One of us (J. A.) set up in 1954 an experiment demonstrating that below threshold reflex klytrons 
behave as negative conductances, and that, in confirmation of the Thompson and others discovery just 
cited, the currrent is strongly sub-Poissonian (unpublished). The expected sub-Poissonian nature of the 
emitted radiation, however, was not observed. 

^''Usually the refiector is raised at a potential slightly lower than that of the emitting cathode to prevent 
electrons from being captured. It also helps finely tune the klystron oscillation frequency. 

^^In conductors the number of electrons is essentially equal to the number of atoms, a huge number. 
It follows that electrons are moving at very low speeds, on the order of l^m/s, even for large currents. 
Accordingly, the above statement that jumps are instantaneous may seem surprising. As a matter of fact, 
electrons appearing on one plate are not the same as the electrons hitting the other plate. In conductors 
the electrical charge should be best viewed as a continuous incompressible fiuid. 
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electronic event, the static potential source delivers an energy eU to the oscillating potential 
source v{t). 

A phenomenon akin to stimulated absorption may be understood similarly. This time, 
we suppose that the electron is emitted by one of the anodes. Without an alternating 
field this electron would remain permanently in the neighborhood of the anodes. However, 
a resonance with the alternating potential may force the electron to oscillate along the 
X axis with increasing amplitude until its energy reaches the value eU, in which case it 
gets captured by one of the cathode. The net effect of this electronic process is that the 
alternating source gives energy to the static source, the opposite of what was discussed in 
the previous paragraph. 

2.12 Quantum sources 

The first man-made quantum oscillators involving discrete matter levels were masers, oper- 
ating at microwave frequencies. Subsequently maser action was discovered to occur natu- 
rally near some stars. The first laser, generating visible light, was discovered by Maiman in 
1959. The fact that space-charge-limited cathodes generate light with sub-Poisson statistics 
was first demonstrated by Teich and Saleh in 1983 |21] . 

The best-known light source is thermal radiation. A hot body like the sun radiates 
energy. The energy inside a closed cavity at absolute temperature Tm contains an energy 
given by the law discovered in 1900 by Planck. An important feature of this law is that 
it involves a previously unknown universal constant h with the dimension of action or 
angular momentum (energy x time). Prom the Quantum Mechanical view point, heat excites 
electrons to atomic levels higher in energy than the ground state energy. These electrons 
then may decay spontaneously to the ground state by emitting a quantum of light. A similar 
mechanism is at work in the so-called "light-emitting diodes" (LED) but the spectrum of 
LED is narrower than that of thermal sources, though much greater than that of lasers. 

A conventional neon tube generates light because the electric discharge excites neon 
atoms that subsequently decay to the ground state, thereby emitting ultraviolet light (sub- 
sequently converted into visible light) by the process of spontaneous emission, similar to 
what happens in thermal sources, but with a narrower spectrum. The so-called Helium- 
Neon laseillj, radiating light at a free-space wavelength of 0.63 ^m, differs from conventional 
neon tubes in that two mirrors located at both ends, and facing each others, force the emit- 
ted light to move back and forth in the tube. Light gets amplified by the process of stim- 
ulated emission, and damped by the process of stimulated absorption. The former exceeds 
the latter when there are more atoms in the higher state than in the lower state (population 
inversion). To achieve this condition the lower-level population must be reduced through 
spontaneous decay to even-lower levels (3-levels lasers). Eventually a steady state of os- 
cillation is reached. The emitted light spreads out in free space as little as is allowed by 
the laws of diffraction, and the laser light is nearly monochromatic (single frequency). The 

^^Helium plays the role of a "buffer gas", allowing the lower neon level to get depopulated. 
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laser linewidth, though small, is of major importance in some applications. Laser diodes 
(also called injection lasers) employ a semi-conductor with a doping that deliver electrons 
(n-type) and a doping that absorbs electrons (p-type) 

To summarize, lasers essentially consist of single-mode resonators containing three-level 
atoms or other forms of matter with a supply of energy called the pump and a sink of 
energy, perhaps an optical detector. As said before, the latter converts the light energy into 
a sequence of electrical pulses corresponding to photo-detection events. When the pump is 
non-fluctuating the emitted light does not fluctuate much. Precisely, this means that the 
variance of the number of photo-detection events observed over a sufficiently long period 
of time is much smaller than the average number of events. Light having that property is 
said to be "sub-Poissonian", or below the "standard quantum level" (SQL). 

2.13 Quantum Electrodynamics 

Historically, the concept that light should be quantized appeared around 1905 on the basis 
of thoughts expressed by Planck and Einstein. Measurements on black-body spectra were 
performed around 1900 with the help of gratings of appropriate periods, and described by 
a formula that involves the universal constant h. On the other hand, the wave properties 
of electrons were discovered only decades later because the concept that electrons might 
possess wave-like behavior ought to wait for the observation that atoms emit light at well- 
defined frequencies, and because of the technical difficulty of sending electrons emitted from 
a small-area source on crystals (playing the role of gratings) in a very good vacuum. The 
interpretation of the observed diffraction patterns involves the constant h. It is perhaps not 
preposterous to suggest that these two key discoveries could have occurred in the reversed 
order. Had this be the case, the Planck constant would have been considered as being 
fundamentally related to atomic behavior, and the subsequent appearance of the same 
constant in black-body radiation would have been viewed as a consequence of the atomic 
theory. 

The non-relativistic approximation. Some ancient philosophers thought that the 
speed of light was infinite. It would then be immaterial to say that light propagates from 
the sun to the eye (say) or the converse. The non-relativistic approximation is applicable to 
an hypothetical world in which the speed of light in free space would be arbitrarily large, the 
other constants {h, Co, e, m, defined in Section [2. 141) remaining as they are. Out of the latter 
constants we may define the electron spin along some quantization axis ±h/2, the Bohr ra- 
dius Co = 'ineoh'^ / me^ ^ 0.53 10"^'' meters, a speed unit Vq = /A-ntoh = h/aom ^ 2.19 10^ 
meters/second, a nominal metal plasma frequency ujp = Vojao « 4.13 lO^^rad/s, and the 
electron magnetic moment, equal to the Bohr magneton eh/2m. A large part of Physics 
may be obtained on the basis of such a non-relativistic approximation. 
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The finite speed of light. Romer, however, discovered in 1676 through a kind of Doppler 
effect, using the motion of a Jupiter satelhte as a clock, that light propagates at a finite speed 
c « 300 000 km/s. The law of causality, as it is presently understood, then implies that light 
propagates from the sun to the eye, for example. The Maxwell theory of electromagnetic 
waves suggests that radiated heat, as well as light, consists of electromagnetic waves of 
some sort. In 1862 Maxwell wrote "[electromagnetic waves travel] at a speed so nearly 
that of light that it seems we have strong reason to conclude that light itself (including 
radiant heat and other radiations) is an electromagnetic disturbance in the form of waves 
propagated through the electromagnetic field according to electromagnetic laws." 

Note that if Vq, defined above, is written as ac, where a 1/137 denotes the fine- 
structure constant, the non-relativistic approximation amounts to setting a as equal to zero. 
Because of the finite speed of light, a number of small corrections to the non-relativistic 
theory, on the order of a, were observed. In particular, let us cite the spin-orbit splitting 
of electron states in atoms, the electron anomalous magnetic moment, the particle-like 
behavior of 7-rayS, the Casimir effect. A theoretical explanation of these effects is offered 
by Quantum Electrodynamics, in which the optical field is quantized. 

Casimir p2l [23] discovered in 1948 that two conducting plates of area A separated by 
a distance d attract each other with a force 

TT^ Ahc _ 7T^ A 

240 d4 240a d2 47reod2 ^ > 

The first expression was obtained theoretically by Casimir by ascribing an energy huj /2 to 
each electromagnetic mode of frequency uj. However, the same result was later obtained by 
Lifchitz from the Nyquist-like noise currents associated with small losses, letting the losses 
go to zero at the end. The second expression in ([2]) suggests an interpretation of the force 
as the product of the number of cells in the area A and the force that an electron would 
exert on a hole separated from it by the distance d. The Casimir prediction agrees with 
experimental results with an accuracy of one per cent, with d typically equal to Ifim. 

At face value, the Casimir force would become arbitrarily large in the non-relativistic 
approximation, c — > 00. This is not the case, however, because the expression given above 
is valid only if d S> c/ujp, where ujp denotes the nominal plasma frequency in metals defined 
earlier. It follows that in the non-relativistic limit the Casimir force vanishes. Conductors 
attract each other only through the van der Waals forces that binds molecules together. 
Prom now on we only consider the non-relativistic approximation. 

Quantum Optics. Most physicists opinion is that the photon concept is essential to 
understand "non-classical" (e.g., sub-shot-noise) states of light, even in the limit where 

^'^ 7-rays have been observed with an energy fttj « 1 micro-joule. They appear as point particles because 
they are detected at rather precise times, locations, energies and momenta, independently of the nature of 
the absorbing material, whether it be steel or lead, say, and look like ultra-relativistic charged particles. 
But in the non-relativistic regime, light and electrons behaviors are quite distinct. 
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c ^ oo. For Quantum Theories of Light, see e.g., [21 ESj EHJ [27l EHJ EU [30] . Let us recall 
some of the arguments given in favor of the photon concept, which evolved into the modern 
second-quantization procedures. 

Light-quanta (later on called "photons") were introduced by Einstein on the basis of 
the following argument. Consider a collection of two-level atoms in a state of thermal 
equilibrium with the black-body radiation field. When an atom in the upper state decays 
to the lower state by emitting light spontaneously it recoils if the light emission is directed 
but would not if light were radiated (almost) isotropically. Einstein calculations indicate 
that a directed emission is required if the Maxwellian atomic velocity distribution is to be 
recovered. The Quantum Optics theory considers that the light emitted spontaneously by 
an atom may be expanded into quantized spherical waves. This concept, though accurate, 
does not fit well with the above Einstein picture. But the Einstein picture fits well with 
the view that light is emitted only if it is directed toward some absorbei^. 

The Quantum-Optics view point is well expressed by Carmichael (op.cit., p. 1213) in 
those terms: "What role does photo-electric detection actually play in the return of the 
atom to its ground state [in resonance fluorescence] after each photon emission? Indeed, 
what does it means to speak of photon emissions as realized events separate from photo- 
electron detection? my viewpoint is that photo-electric detection does not cause atomic- 
state reduction. Projection of the atom into its ground state is caused by the dissipative 
nature of the atomic dynamics, and reoccurs, on average, at the mean spontaneous-emission 
rate, with complete indifference to the presence or absence of the observer. Photo-electric 
detection merely monitors emitted (realized) photons. It does not intrude into a coherent 
Quantum Dynamics in the manner implied by a measurement-induced wave-packet reduc- 
tion; it is the irreversible decay into the vacuum that interrupts the coherence of the source 
dynamics. No doubt, in general outline, this view-point is already widely held". One may 
wonder whether the conceptual ambiguity outlined in the above quotation is not a conse- 
quence of the requirement that energy be conserved exactly (one photon being absorbed 
whenever a photo-detection event occurs) instead of being conserved on the average only. 

Another argument often given in favor of the concept that light consists of lumps of 
energy Huj is the observation that when a light beam of constant small intensity is incident 
on an ideal photo-detector (i.e., free of dark current and thermal noise) photo-current events 
sometimes occur long before the required optical energy Hu) has been collected, in apparent 
violation of the law of conservation of energy. However, one should require that only the law 
of average energy conservation be enforced. For a single system, there exists no independent 
way of measuring the "light intensity" as a function of time. The only information one may 
obtain concerning the intensity of a light beam is through the output of photo-detectors, 
and this brings us back simply to the observation made. A number of authors have shown 

^*For an arbitrary distance R between an emitting and an absorbing atom, see |31j . These authors employ 
orthodox Quantum Electrodynamics and comment that "In a sense, every photon is virtual, being emitted 
and then, sooner or later, absorbed", and "virtual photons are messenger particles that cannot be directly 
detected". 
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that many effects that were at a time supposed to prove the reality of the photon concept 
such as the photo-electric effect, may in fact be interpreted in a semi-classical manner. 

Most Quantum Theories of laser action begin with a discussion of the statistical oper- 
ator p of the optical field in empty loss-less resonators, which are treated in analogy with 
mechanical oscillators. Pumping and losses are subsequently introduced in an approximate 
manner by enforcing the preservation of the commutation relations. For the case presently 
considered (stationary linearized laser) the results of such calculations exactly coincide with 
ours. One must thus presume that the approximations are the same. 

Other Quantum Optics treatments consider atoms in either their upper (pumping 
atoms) or lower (detecting atoms) states introduced at specific times into the optical cavity 
and spending there a fixed time r. Whether the atoms leaving the cavity are in their lower 
or upper state may in principle be measured, and the corresponding probabilities may be 
evaluated. In that way some properties of stationary lasers may be predicted. However, 
the flying-atoms configuration is quite different from the one discussed in the major part 
of this paper. 

Another kind of Quantum-Optics theory considers instead continuous resonant photo- 
detection processes. That is, atoms in the absorbing state are coupled to the optical 
cavity at all times. This is the configuration considered in this paper. A recent treat- 
ment [32] employs the operators introduced by London in 1926 having the property that 
E± \m) = |mibl), where \m) denotes a state with exactly m photons. To our knowledge, 
this continuous-detection theory has not been applied to quiet lasers. This is probably the 
kind of Quantum Optics theory that one should employ to establish in a logical manner a 
link with the present Semi-Classical theory (i.e., rather than just comparing final formulas). 

As far as electrons are concerned, physicists "second-quantize" the Schrodinger wave- 
function. Electron second-quantization is a "book keeping" method, convenient when the 
optical field itself is quantized, but presently unnecessary. The statistical properties of non- 
interacting electrons that may exchange heat with a reservoir (this is the canonical ensemble, 
to be distinguished from the grand-canonical ensemble) have been obtained from second 
quantization of the electronic wave function. However, a much simpler direct solution has 
also been found, based on the partition of integers [33] . 

Physical paradigms. A physical "paradigm" rests on a number of universal constants, 
on particles parameters, and on recipes to relate the theory to observations. For example 
Newtonian Celestial Mechanics employs a single universal constant, namely G, point par- 
ticles have as sole parameter their mass (referred to that of a particular piece of platinum, 
called the kilogram). Given the position and speed of the particles at a given time, the 
theory provides their positions and speed at all times. 

The present theory employs as universal constants h, and, in some circumstances, 
k-Q. However, G and c are not used. Particles (electrons) are characterized by their mass, 
electrical charge, and, in some circumstances, by spin and magnetic moment. Given the 
elements constitutive of a particular device, we determine the statistics of photo-electrons. 
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which may, in principle, be measured with unlimited accuracy. 
2.14 Units and notations 

Our notations and conventions may differ from those employed by engineers, physicists, or 
experimentalists, which are not fully consistent. We attempted to follow the majority rule 
unless this leads to confusion. To simplify formulas we sometimes set as unity quantities 
such as the characteristic conductance of transmission lines. Otherwise, SI units (see below) 
are employed throughout. 

Numerical values. Useful numerical values in the realm of Non-Relativistic Physics are 

G (Newton gravitational constant) ~ 6.67 10"''^^ SI 
g (earth gravitational acceleration at see level) ^ 9.81 meters per second squared 
e (absolute electron charge) ^ 1.60 10~^^ coulombs 
m (electron mass) ^ 9.10 10"^^ kilograms 
h (Planck constant divided by 27r) ^ 1.05 10"^^ joules x second 

k-Q (Boltzmann constant) ~ 1.38 10"^^ joules/kelvin (3) 

and 

^— = 10-^(2.99792458 10^)^ farads/meter (4) 
47reo ^ ^ ' 

which is exact, i.e., not subjected to revision as a consequence of later measurements, and 
involves a finite number of digits. The constants e, m, ^, 47reo, A;b are the only ones that 
enter into the present theory because, in agreement with the non-relativistic approximation, 
we set the free-space permeability /Ho = or equivalently c = oo. We sometimes employ as 
energy unit the electron-voltRi 1.60 lO-"*^^ joules. Note that capacitances may be evaluated 
from Eo and their geometric dimensions, namely the electrodes areas and their spacings. As 
far as inductances are concerned we suppose that they are known from measurement. 

Notation. We list here only some of the notation employed in this paper. Different func- 
tions are distinguished by explicitly writing out their arguments. For example the Fourier 
transform of a function is denoted by '0(/u), i.e., with the same symbol, even though 
these are different functions. When the arguments are similar tlie same function, however, 
is intended. For example w{t) and w{t) represent the same function. One should not con- 
fuse a constant U (no argument) with a function U{x), for example. To avoid confusion 
between the arguments of a function and products we sometimes employ parentheses of 
different size, e.g., f{x/d) is a function / of x/d, while f(x/d) represents the product of / 
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and x/d. As usual, cos^(x) means (cos(x)) , and likewise for other trigonometric functions. 
iIj*{x) is equivalent to (■(/'(x))*. 

For a two-state electron the lower and upper levels ("working levels") are denoted 1 and 
2, or "a" (absorbing) and "e" ("emitting") levels, respectively. For 4-level electrons, levels 
of increasing energy are denoted 0,1,2,3. Pumping occurs from to 3. The spontaneous- 
emission rate of an electron from 2 to 1 is denoted by l/r^ = 27. It is often denoted in 
other works by 2/3, 7, or F. The decay time of an oscillator with loss is denoted by Tp (the 
subscript "p" stands for "photon", even though the word is not always appropriate). 

Probabilities are denoted by pr{a > b) and probability densities by P{x), for example, 
but the letter P may also refer to power. We employ double-side spectral densities, so that 
the usual shot-noise formula 2e \J\ is written as e | J|, i.e., without a factor of 2 (— e denotes 
the electron charge and J the average current). We will perform three main kinds of aver- 
aging. First a Quantum Mechanical (QM) averaging over an ensemble of similarly prepared 
systems (if these systems exhibit macroscopic differences, these are supposed to be too small 
to affect the QM averages). Second, statistical averaging. We consider here an ensemble of 
similarly-prepared systems that differ from one another significantly. Third, an averaging 
over an optical period. These averagings could be denoted as \\{P)qm) ) 1 

\ \ V / ^jmj / Opt.period 

where, as an example, P represents the power supplied by an optical potential source. This 
quantity may still be a slowly- varying function of time, which we denote for brevity P{t). 
In some cases, for stationary systems, we may further introduce an averaging of P{t) over 
an arbitrarily large period of time. Other notations will be outlined in Section [3l 

Decibels. We recall here a notation commonly employed in Electrical Engineering. Usu- 
ally, an amplifier is loaded with a nominal conductance such as 20 milli-siemens (resistance 
of 50 O) and the input impedance is equal to that of the load. If the input power of an 
amplifier is Pin and the output power is Pout, the amplifier gain in decibel (abbreviation 
"dB") is defined as 101ogxo(^out/^m)- If the amplifier is linear the gain does not depend 
on the input power. In terms of the potentials Vin,Vout at the input and output ports, 
the gain reads 20logiQ{\Vout/Vin\) dB, because powers are proportional to the modulus- 
squares of the potentials in the situation considered. Likewise, an attenuation is defined 
as WlogiQ{Pin/ Pout) = — 101ogio(-Poitt/-Pm)- A gain of 3dB means that the input power is 
multiplied by a factor close to 2. Note that dBm means decibels above a power of 1 mW. 
For example, 30dBm represents approximately a power of 1 watt. 

Electrical circuit schematics. In schematics, current sources are represented by circles 
with an arrow in them, while potential sources are represented by a circle with + and - signs, 
to define what is meant by positive current or positive potential, as shown later in Fig. [H 
By potential (or current) sources we mean potentials (or currents) that do not depend on 
the current (potential) delivered. These are called in Quantum Optics "prescribed classical 
sources". 
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3 Mathematics 



Ideally, the present section should derive from axioms all the mathematical results sub- 
sequently employed. This goal is not accomplished for lack of space, time (and ability). 
We therefore provide references to text books. A difficulty is that many text-books con- 
sider general situations from which it is not always easy to retrieve the desired information. 
Often, physical intuition or numerical calculations help. For example the celebrated Wiener- 
Khintchine theorem that relates spectrum and correlation was first obtained by Einstein in 
an intuitive manner. 

We will clarify in the present section the notation employed, recall elementary mathe- 
matical formulas, and give the main properties of random point-processes. We then offer a 
simple picture of quiet light generation, and show that random deletion of photo-electrons 
leaves unaffected the reduced spectrum. This section is entitled "Basic Mathematics" be- 
cause most results are mathematical in nature, although rigor is overlooked and the physical 
motivation is often pointed out. 

3.1 Complex numbers 

A complex number is denoted either as z = R{z}+\Q{z} or as z = z' +\z" , and z* = z' — iz" 
denotes the complex conjugate of z. We denote = zz* = z'^ + z"'^. A complex notation 
is often employed for describing quantities that vary sinusoidally in time. According to 
the complex notation, the function i{t) = ^/2\I\ cos{ujt + (/>), where the frequency uj and 
the phase 4> are real constants, is written as i(t) = V2^{I exp{—iu:t)} , where the complex 
number I is defined as I = |/| exp(— Similar definitions apply to potentials v{t) varying 
sinusoidally in time, that is v{t) = \/25R{y exp(— icjt)}. We have chosen to introduce the 
factor \/2 so that the average optical power, defined as the time average of the current- 
potential product v{t)i{t) be simply equal to the real part of the product VI*, i.e., without 
the factor 1/2 that would otherwise occur. V and I are called rms (root-mean-square) com- 
plex potentials and currents, respectively, or more briefiy, optical potentials and currents. 
The minus sign in exp(— iwt) is employed in optics because waves propagating forward in 
space then involve a term of the form exp(ifcx), where k denotes the wave-number, that 
is, with a plus sign. But for slow variations the function j{t) = V^\J\ cos(Ot -|- (f)) is de- 
noted j{t) = \/23f?{ Jexp(jr2t)}, where the complex number J is defined as J = \ J\ exp(j(/)), 
as is usually done in electrical engineering. The complex notation considerably simplifies 
calculations for real, causal, linear and time-invariant circuits submitted to sinusoidal po- 
tentials or currents. Even though the squares of i and j are both equal to -1 these two 
numbers should be distinguished. A "bi-complex" notation describes in an exact manner 
sinusoidally-modulated sinusoidal signals. It is recalled in 13.151 

Second and third-degree equation. Elementary algebraic operations are concisely 
recalled, for example, in [3l]. We give below the solution of second and third-degree equation 
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in the form appropriate to our intended application. 

Consider the polynomial ap^ + hp + c = Q. The solutions are 



—h± \/}p- — 4ac , , 

^± = 2a • 



For the third-degree polynomial + a2]? + a\-p + ao = we evaluate sequentially 



oi aia2 oo 0,2 



si = {s + r)^/^ exp(i7r/3) S2 = {s - r)^/^ exp(i27r/3) s ± r > 0. (6) 

The three roots are 

Pi = S1 + S2- Y 

P2 = si exp(i27r/3) — S2 exp(i7r/3) — ^ 

P3 = -Si exp(i7r/3) + S2 exp(i27r/3) - y. (7) 
3.2 Vectors and matrices 

W'e consider a fixed coordinate system and define a vector (bold-face letter) as a collection 
of complex numbers denoted for example 

Transposition, indicated by a "t" in upperscript, interchanges lines and columns. For 
example 

a* = ( ai 02 ) . (9) 
The scalar product of two vectors a and b is defined as 

a.b = ai 61 + 02 62 (10) 
The modulus square of the length of a complex vector is 

laf = a.a* = ai a\ + 02 > 0. (11) 
A matrix is denoted for example 

\^ M21 M22 y ^ ^ 
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The trace of a square matrix is the sum of the diagonal elements 

trace{M} = Afn + M22, (13) 

and the determinant 

det{M} = M11M22 - M12M21. (14) 

The sum of two matrices is obtained by summing their elements. The product L of two 
matrices M and AT" is for example 

Lij = MiiNij + Mi2M2j i, J = 1, 2. (15) 

We have(ML)* = L*M*. M'^ denotes a matrix such that M'^M = 1 = 

When det{M} = the matrix is singular, and cannot be inverted. A matrix is said to be 
symmetrical when = M, Hermitian when M^* = M, unitary when MM^* = 1. We 
have trace{AB} = trace{BA}. The trace of the product of two Hermitian matrices is 
real. 

Cauchy- Schwartz inequality. Consider two unit vectors a, 5, that is, such that \a\ = 
\b\ = 1, The Cauchy-Schwartz inequality reads 

\a.b*f = {a.b*){b.a*) < 1. (16) 

This relation is obtained by replacing a in (fTTIl by a — (a.6*)b. 

Density matrices. Let p denote an Hermitian matrix of trace 1. 

The pure-state density matrix is constructed from the vector a with \a\ = 1 as 

P = f "0 («t «2) = ( "^1 "^1 ) . (17) 

We readily find that = p, and thus trace^p^^ = trace{p} = 1. Note that p is unaffected 
by a change of the phase of a. As examples, we may have 

Let Pa, Pb be two such matrices. By explicit calculation we find that trace{pa Pb} = 
{a.b*){b.a*). Thus, from (fTBIl we have trace{pa Pb} < 1- 
The mixed-state density matrix is defined as 

P = ^PkPk Pk>0 ^pfc = l, (19) 

k k 
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where the pk may be called weights. Since trace{pj. pi} < 1 we obtain 

trace{p^} = traceiJ^Pk^PlPkPl} <J2p^^^ '^^^^ 

k I k 

If trace{p'^} = 1 we have (a.6*)(6.a*) = 1. The density matrix p is then of the pure-state 
form in ([I7| . 

Let us consider a pure-state density matrix pj^, and suppose that the QM (Quantum 
Mechanical) average of some quantity, such as the power Pk, may be obtained from the 
formula 

iPk)QM = trace{pkP] (21) 

where P denotes some known Hermitian 2x2 matrix. Next, suppose that the pure-state 
density matrix p^ occurs with probability prk, k = 1,2.... The statistical and QM-average 
of the power, denoted by a double bracket is, using the properties of the trace 

\'\P)qm) =^prktrace{pkP} = trace{pP} p = ^prkPk- (22) 

\ / statistical — ' — ' 

k k 

Note that the mixed-state density matrix is of the form in (fT9]) with p^ = pr^- Prom now on, 
the double bracket is replaced by a simple bracket. The diagonal elements pii,p22, called 
"populations", are non-negative and sum up to 1. The off-diagonal elements pi2 = P21 ^'^^ 
called "coherences". 

Setting X = 2/>'^2) U = "^P'vi^ ^ = /'22 — Pii (unrelated to coordinates in space), the density 
matrix may be written as 

1/ \ — z X + \y 



It follows that 



• , • (23) 

2\x — \yl + z' 



trace{p'^} = ]^{l + x^ + y'^ + z^), (24) 

and trace{p'^} < 1 <;=^ x'^ + y'^ + z"^ <l. 

If p depends on time, ■^trace{p'^} = x^ + y^ + z^. It follows that if at some 
time we have a pure-state density matrix, x'^ + y'^ + z'^ = 1, we must have at that time 

+ + - ^- ^^^^^ suppose that 



dx 

'dt 

dy 

dt 

dz 

'dt 



—2a^x 
—z — 2a^y 

y-2-fz + 267, (25) 
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where 0,7 are non-negative parameters and — 1 < 6 < 1. The above condition holds 
provided that 

{a-l)z^ + bz-a<0 \z\<l, (26) 

that is, 2a > 1 — V 1 — • When 5 = ±1 we must have 2a > 1. Physical arguments given 
later on show that for a one-electron model 2a is in fact unity. 

3.3 Fourier transforms 

The Fourier transform ^l^{k) of the function ip(x) and the reciprocal relation are 

POO 

ip{k) = / dx exp{—ikx)ip{x) (27) 



-00 
00 



dk 

^/;(x) = / — exp{xkx)i;{k), (28) 



00 



where = —1, and k is called the wave- number. Note that the position of the 2tt factor 
varies from one author and another, without of course affecting the end results. Obviously, 
the Fourier transform of [ik)^ip{k) is equal to the nth derivative of ip{x) with respect to x. 
If ip{x) is real, we have ip*{k) = ip{—k). 

Relations similar to ^7} and ([28]) hold with x changed to t, k to u and (to be consistent 
with our conventions for optical signals) i changed to — i. Then the element of integration 
in ([28]) is dv = du/2TT^ where denotes as usual the optical frequency 

/oo 
dtexp(\u;t)ij{t) (29) 
-00 

/oo 
— exp(-iu;t)V^(w), (30) 

In the Fourier-frequency domain, changing uj ^ Q,\ ^ — j, we write 

/oo 
dtexp{-jnt)Tp{t) (31) 
-00 

f°° dfl 

m = J — expOQi)V^(O), (32) 

Note the following physical application. For particles moving in time-independent po- 
tentials V{x), stationary states tp{x) are real functions of x (to within an arbitrary over-all 
phase factor that we set equal to 1). If furthermore V{x) is an even function of x, ip{x) is 
either an even or odd function of x. It follows from the above considerations that the tp{k)- 
functions are, respectively, real even or imaginary odd. In the present mathematical section 
we set h = 1 and do not distinguish the electron momentum p from the wave-number k. 
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We will need the following expression of the Dirac 5-distributio 

/°° dk 
— ex-p{ikx), (33) 

implying that its Fourier transform is unity. Using this expression one may prove that 

+00 ^r, r+00 

— \i^{k)f= dxmx)\\ (34) 

Let \tp{x)\'^ dx be interpreted as the probability of finding the position of a particle 
between x and x + dx if & measurement is performed, and P{p)dp = \ip{p)\'^ dp be the 
probability of finding the electron momentum between p and p + dp ii a measurement is 
performed (this latter measurement may be accomplished by letting the particle free at 
some time t and observing its position on some far-away screen). We are led to define the 
wave function in momentum space as 

]^ r+00 

ip{p) = — / dx e'x.p{—\px)'il^{x), (35) 
v27r J-00 



according to (f34ll . In general the wave-function depends on time, and the average value of 
p, evaluated at some time i, depends on time and is denoted {pit)). 

For two functions tp{x) and (j){x) and their respective Fourier transforms ip{k) and (^(fc) 
we obtain from the expression in (f33]l of the d{.) distribution the identity 

-^{k)(t)'{k) = / dxi){x)(t)*{x). (36) 



—00 J —00 

Transformations. If we set in (f36]l V'(^) = '0i(^) and 0(A;) = k^ip2{k), we have, using 
the observation following (f28ll that the Fourier transform of (ik)"\l:{k) is equal to the nth 
derivative of ^p{x) with respect to x 

—k^i,,(k)r2{k) = / Vi(x)(— )>2*(x), (37) 



'12 



provided the integrals exist. 

Through two integrations by parts we obtain that 



V(x)(-^)V(x) = / (/)(a;)(— )V(2;). (38) 



di'jc J (yc) dx 

On the other hand, for any derivable function ^p{x 



^^The 5(t)-distribution may be viewed alternatively as a function equal to 1/h for —h/2 < t < h/2 and 
otherwise, so that the area under the function is unity, letting h go to zero at the end of the calculations. 
Many other forms of the 5-function may be used, with less-singular derivatives than for the one just given. 
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Average wave-number. Setting n = 1 in ([371) we find that 

—kMmik) = -iJ dx (40) 
Using ([39]) with tpix) = ^^2(2;), the above expression may be written as 

W12 = dx Mx) - x{^)'^ r^ix). (41) 

Finally, employing ([38]) with il){x) = tpi{x) and (/)(x) = X'i/'2(^) we obtain 

^^^12 = "2 y ^ (^2^^^)^^ " ^^^"^^-^ j ■ ^^^^ 

In particular, we have from (j42]) . suppressing the subscripts, the relation 



00 



(A;) ^ / _A:V;(/c)^*(A:) = - / :e ( V(x)^^ - V'*(a:)^^ ) , (43) 



which will prove useful in establishing the first Ehrenfest equation. 

Eigen- functions. If tjji{x) and ■02(2;) are solutions of the eigen-equations d'^ipi^2{x)/dx'^ + 
ei 2fpi 2 = 0, we obtain from (j42]) that 



ik)i2 = 2 J ^ Mx^ix) ^ ^ \ ^' xi2. (44) 

On the other hand, setting n = 2 in ([37[) we obtain that (i,j = 1,2) 

/+00 jr, /•+00 d^ib*(x) r+00 

—^,{k)er,{k) = - dx^x)^^ = e] I dxMx)i^-ix). 
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The above expressions will be employed in Section [5^ in relation with the time-dependent 
Schrodinger equation. 

3.4 Convolution and Laplace transforms 

Let us consider a real, causal, linear and time-invariant system. These conditions imply 
that for a potential source v{t) the current i{t) (or more generally the response to a source 
of any kind) is given by 

r+00 

i{t) = I du h{u)v{t — u) = h*v = h*y, (46) 
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where the kernel h{u) is real, equal to for u < 0, and middle stars denote convolution 
products. Convolutions are associative, so that parentheses in convolution products are 
unnecessary, and commutative. For example, for a conductance G, we have i{t) = Gv{t), 
and thus h{t) = G6{t), where 6{.) denotes the Dirac distribution. If v{t) = V{p) ex.p{pt) 
where p denotes a complex number (not to be confused with particle momenta), i{t) = 
/(p)exp(pt), where I{p) = H{p)V{p) and 

H{p)= / dt exp{-pt)h{t), (47) 
Jo 

a Laplace transform, defines H{p) for complex p. In most of this work we set p = — iw, and 
H{—iu!) is denoted Y{uj) and called the admittance. 

In particular the Laplace transform of exp(At) is l/(p — A). It follows that if the 
reciprocal of a polynomial in p may be written as a sum of terms of the form l/{p—pk) the 
inverse Laplace transform is easily obtained. More generally, the Heaviside theorem says 
that if f{p) is a polynomial with distinct roots (not to be confused with probabilities) pk, 
k = 1,2. ..n {f{pk) = 0), the inverse Laplace transform of 1/ f{p) is 

r-ir 1 1 _ V exp(pfct) 

The Laplace transform of the convolution of any number of functions is the product of 
their Laplace transforms. For example, the Laplace transform of Gfc(t) = w{t)*w{t)...*w{t) 
(fe-times) is the fcth power w{p)^ oiw{p), where w{p) denotes the Laplace transform oiw{t). 
Thus, the Laplace transform of G{t) = w{t) + 'w{t) * w{t) + .... is the sum of an infinite 
geometric series fSS^, p. 53] 

Gip) = (49) 

1 — w{p) 

Let us define an average waiting time 

{t)^fdttn^it) = -{^)^^^. (50) 

If (l49]l holds, G{t — > oo) is finite and the other terms are decaying exponentials, so that 
G{p — > 0) « G{t = oo)/p, we obtain 

l- = G{t = ^). (51) 
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3.5 Relative noise, normalized correlation and variance. 

In the present section we clarify the notation employed, and explain why it differs from 
the one used in other works. More precise definitions and result are given in a subsequent 
section. 

We are mostly concerned with photo-currents j{t) = —eV{t). Here, V{t) is the sum 
over k of 5{t — ifc)-distributions where the tk form a stationary point (or "event") process 
of average rate D. We call "reduced spectrum" S^AoiP) the spectral density S^t>{^) of 
that process with the singularity at = removed. The "relative spectrum" ^ad/£)(^) is 
obtained by dividing ^Aoi^) by -D^. Finally, the "relative noise" M{^) = ^ad/d(^)-1/^ 
vanishes for shot-noise, also referred to as the "standard quantum limit". We are mostly 
interested in circumstances where the relative noise is negative. In the engineering literature 
most authors call "relative-intensity noise" the quantity 2oS^ad/d(^)5 and express it in 
decibels/hertz. Our main objection to using it is that, unlike the relative noise employed 
in this paper, the relative-intensity noise does not enjoy the property of being independent 
of (cold, linear) attenuations. Further, the expression "decibel per hertz" is difhcult to 
comprehend. 

We call "normalized correlation" g{T) the auto-correlation of the V{t) process with the 
singularity at r = being removed, divided by the square of the average rate D. It may 
take any non-negative value. The reason why we do not use the Quantum Optics notation 
5^^)(t) is two-fold. One is that correlations of order other than the second are not employed, 
so that no confusion may arise. More importantly, (7^^-*(t) is usually defined in terms of 
the so-called "light intensity" I{t) according to ^'•^•'(t) = {I{0)I{t)) //^. For mathematical 
reasons a quantity so-defined cannot be less than unity. Since g^"^^ (r)-values less than unity 
are apparently measured, the Quantum-Optics view point is that I{t) should be considered 
as an operator instead of an ordinary function of time. In this paper the concept of "light 
intensity" does not enter and y(r) refers to photo-currents exclusively. 

We denote by V(r) the variance of the number of photo-detection events occurring 
within the time-interval T, divided by D, minus 1. In Quantum Optics, the Mandel Q- 
parameter is defined by a similarly looking expression, namely Q = variancc(m)/ (m) — 1, 
but the operator m sometimes refers to the number of photons in the cavity, rather than 
to the number of photo-electrons. Thus our V(r) and the Q-factor may have different 
physical meanings. Some authors employ the Fano factor = Q + 1. Sometimes in the 
literature, however, the Fano factor refers, not to the number of photons in a cavity, but to 
the normalized spectrum of electrical-current fluctuations. Again, because of this ambiguity 
we set the (double-sided) spectral density of the current driving a laser diode as ^ (J), where 
^ = for a quiet pump and ^ = 1 for a Poissonian pump, and consider that the Fano factor 
^ relates to the number of photons (ratio of energy and Hloo) in an optical resonator. 
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3.6 Random processes 

Instead of attempting to clarify the meaning of the word "probability", let us give the 
axioms. Given two events A and B, A + B denotes an event that occurs when A or B, or 
both, occur. A probability is denoted by "pr(.)". The theory rests on the two following 
axioms: 

• The probability pr{A) of an event ^ is a non-negative number, unity for a sure event. 

• If ^ and B are mutually exclusive (i.e., the occurrence of one at a given trial precludes 
the occurrence of the other) pr{A + B) = pr{A) +pr{B). 

The product AB of two sets of events A and B is the set of the events that are common to 
A and B. Two events are called independent \l pr{AB) = pr{A)pr{B). 

Random variable To every outcome ^ of an experiment (such as rolling a dice) we 
associate a number x{Q, called a random variable, {x < X} denotes the set of all outcomes 
C such that x(C) < X. The distribution function F{X) is the probability pr{x < X} that 
x{() be less than or equal to X. The probability density is defined as P{X) = dF{X)/dX. 
Since, in Physics, X has usually a dimension (e.g., time) the dimension of P{X) is the 
reciprocal of that of X. 

Stochastic process We are given an experiment such as rolling a dice, specified by its 
outcome C,. To every outcome we assign a time function x{t;C,). This is called a stochastic 
process. 

Wide-sense stationarity A real process x{t) is said to be "wide-sense stationary" when 

{x{t)) and the correlation R{T,t) = {x{t)x{t + t)) do not depend on t. It follows that for 
real wide-sense stationary random processes R{t) = R{—t). To prove it, set t = — r and 
remember that x{t) is an ordinary function of time so that x{t) and x(t + r) commute. 

Ergodicity Usually one observes the photo-current j{t) from a photo-detector in a single 
set-up. The mathematical treatment of noise, on the other hand, rests on the consideration 
of an arbitrarily large number of macroscopically-identical set ups, averaging referring to 
these many set-ups. The question thus arises as to what practical conclusions may be 
drown from formulas derived from the formalism. The answer is that, provided a stationary 
system is ergodic, the statistics may be obtained from a single set-up. If this is the case, 
statistical averages are equivalent to time averages. A system may be ergodic only if the 
correlation (j(0)j('7")) tends to (j(0))^ as r tends to infinity. More stringent conditions must 
be fulfilled, however, that we suppose met. In principle, once the statistical calculations 
have been performed, one should verify that the system considered is ergodic in order to 
be able to apply the results to a single system. In practice, this step is omitted. 
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3.7 Spectrum 

The spectrum ^ti^) of a process is a real non-negative even function of 0. It may be 
obtained by considering a finite duration T, evaluating of average of the modulus square of 
the Fourier transform of x(t), dividing by T, and letting T go to infinity. The motivation 
for introducing a subscript "i" is that we intend to introduce later on a spectrum denoted 
y{^l) (without a subscript) obtained from =^^(0) by removing the singularity at = 0. 
Alternatively, the spectrum may be expressed as the Fourier transform of R{t) (Wiener- 
Khintchine theorem) 



^t(O) = lim - 



+ 00 



dt x{t) exp(jr2t) 



dr R{t) exp(ji7T) 



+ 00 



dr R{t) cos(ilr). 



(52) 



The two above expressions agree in the mean if the integral from to oo of tR{t) is finite 
[361 p. 336]. We have employed above the electrical-engineering exp(jQt) notation. To 
compare with the previous notation in ([27l) change Q. to k, t to x, and j to — i. Evaluating 

^(x(t) ibx(O))^^ we notice that — i?(0) < R{t) < i?(0). This condition does not suffice 

however to make R{t) positive definite, that is, its Fourier transform could still be negative. 

In the special case where x{t) does not depend on time we have R{t) = {x'^)= constant. 
Substituting in ([52]) we find that ^t{^) = 27r (x)^ ^($1), where 5{.) denotes the Dirac 5- 
distribution. We are thus led to define a reduced spectrum S^{VL) = S^t{^) — 27r (x)^ 5{VL). 

Conversely, the correlation may be expressed in terms of the spectrum through the 
inverse Fourier transform according to 



Rir) 



+CX3 



dn 

2^ 



St{n) cos(J]t) 



(53) 



If we define y{t) = x{t) - {x{t)) we have {y{t)) = 0. The function C(t) = (y(O)y(r)) is 
called the (auto) covariance of the process x(t). 

We are mostly interested in the case where x{t) = J2k ^(^ ~ with < tk < T, where 
the tk are referred to as "points" with a constant density D (point processes are discussed in 
more detail in section [3^ . We have {x{t)) = D. The spectrum of x(i) exhibits a singularity 
2'kD'^5{^) at O = 0. The first expression in ([52]) give the reduced spectrum 



lim — 

T-,oo T 



where = 27rn/r, n = 1,2. 
as one wishes by setting n 



allowed k 



(54) 



Note that n = is not allowed, but O can be made as small 
1 and letting T go to infinity. This expression is useful to 
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evaluate spectra through numerical calculations that generate runs, each with a different 
tk sequence. Because of the assumed ergodicity, a single run suffices. Total energies are 
obtained by integrating over frequency from — oo to +00 the quantity ((ir2/27r)^(0). For 
two independent processes x{t) and y{t) of spectral densities and S^y, respectively, the 
spectral density of z{t) = ax{t) + by{t) is J>^z = \o\^ -S^x + |&|^ -S^y 

The moments of a quantity such as x, denoted (x"), are defined as the integrals over x 
from —00 to +00 of x"P(x), where n = 1,2... and P{x) denotes the probability density of 
X, or the sum from A; = 1 to 00 of k"'pr{k)^ where pr{k) denotes the probability of having 
the outcome k. Note that the above integrals or sums may not exist, even for well-behaved 
probability laws. It follows from (f52]l and ([27|l that 



2\ dn 



{^') = C(0) = / ^^(f^) (55) 



00 



if (x) = 0. Thus, if the spectrum shape is known, the variance of x determines the spectrum. 

The sources of noise in our theory are narrow-band current sources written as c{t) = 
^/2{C' {t) cos{u)o't) + C" {t)sui{u]ot)) ^ whcrc ujo denotes the average laser frequency, and 
C'{t), C"{t) are jointly-stationary slowly-varying real functions of time. It can be shown 
that c{t) is wide-sense stationary if and only if {C'{t)) = {C"{t)) = and the auto and cross 
correlations fulfill the conditions Rc'C'i'^) = Rc"C"{'^)-,Rc'C"{'^) = —Rc'C'i'^)- Further- 
more, we assume that the statistics is independent of a phase change, and this entails that 
Rc'C"{t) = 0- Let us recall the following result. If C'{t), C"{t) are uncorrelated and their 
spectra yduj) = 5^c"iy) vanish for |a;| > Wc, then 5^c^^) = -5^C'('^ ~ '^o) + •^C'iy' + ^o) 
pi p. 380]. 

Consider a constant-amplitude frequency-modulated signal x(t) = zo's.iuot -|- </)(t)). If 
Aa;(t) = d4){t)/dt is a stationary low-frequency gaussian process of (double-sided) spectral 
density ^S^Aoj, the spectral density yx{^) of x{t) is Lorentzian with a full-width at half power 
(FWHP) 5u = y/^^ |37I p. 140], that is, is of the form y^ioj) (xl/[l + {2{uj - uJo)/5^f]. 
Note that here Aw(t) denotes a process, while 5uj is a real positive number, the spectral 
width. 



3.8 Point processes 

For physical motivation note that, ideally, photo-currents in an experiment lasting from 
i = to t = T are of the form j{t) = -eV{t), V{t) = D + AD{t) = J^k^i^ - tk) 
with < tfc < T, where — e denotes the electron charge and 5{.) the Dirac distribution. 
Experimentally, we may measure a number of quantities relating to j{t) with the help of 
integrators, narrow-band filters or electron counters. 

Point processes are sequences of increasing positive real numbers t^, A; = 1, 2 . . . . Each k 
value corresponds to a "point" occurring at time tfc. We consider M such sequences, labeled 
by m = 1,2...M, called runs. As said above, averages denoted by the sign (.) refer to sums 
from m = 1 to m = M of some quantity defined for each run divided by M, letting M go to 
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infinity. Clearly, averaging is a linear operation, that is {a + b)={a) + {b) . Point processes 
are fully defined either by the so-called "non-exclusive multi-coincidence rates", or from the 
statistics of durations between successive points (related to the "exclusive multi-coincidence 
rates"). Stationary point processes have a time-independent density D. Let d{t) denote the 
number of points occurring up to time t, that is the number of k values such that tk < t. 
The (rather intuitive) result that {d{T)) = DT is demonstrated in Section [3.131 

Poisson processes The most important process is the Poisson process. Consider a Pois- 
son process of density 1. The density that the first point occurs at time r > is exp(— r), 
whether or not there is a point at t = 0. If there is a point at t = 0, w{t) = exp(— r) is 
the waiting-time density. Consider next an inhomogeneous Poisson process of density X{t) 
with a point at t = 0. It may be reduced to a Poisson process through a transformation of 
the time scale dr = X{t)dt. Thus the waiting-time density, that is, the probability that the 
next point occurs in the interval (r, r + dr), divided by dr, reads, 

w{t) = X{t) exp{- [ dtX{t)) = --^exp{- [ dtX{t)). (56) 
Jo «T Jo 

It follows from the second form above that the integral of w{t) from to oo is unity, 
provided X{t) does not tend to as i — > oo. This means that the point eventually occurs. 
The average duration between adjacent points is, after an integration by parts 

(r) = / drexpi- / dtX{t)). (57) 
Jo Jo 

If, for example, X{t) = 1, we obtain (r) = 1 as expected. Conversely, 

^(''^jffe—s'^Kr''™'*"')- '''' 

Ordinary renewal process With to = 0,ti = Ti,t2 = ti +T2i--- and the rj,i = 1,2... 
independent and distributed according to the same density w{t), one generates an ordinary 
renewal process t^, k = 1, 2.... Such a process is non-stationary, but it tends to be stationary 
for large times. 
Let us denote 

w{p) = / dt ex-p{—pt)w{t) (59) 
Jo 

the Laplace transform of the waiting-time density. Given that there is a point at i = 0, the 
probability G{t)dt that there is a point between t and t + dt, is the sum of the probabilities 
that this occurs through one jump, two jumps, ...Because the jumps are independent and 
have the same densities, we obtain the Laplace transform of G{t) 

/•oo 

G{p) = / dtexp{-pt)G{t) (60) 

^0 
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in the form, see Section (|3.8|) . 



G(p) = (61) 

1 — w{p) 

It follows that, given the density X{t) of an inhomogeneous Poisson process, viewed as a 
renewal process, we may in principle obtain the Laplace transform of the auto-correlation 
function G{t). As an example, suppose that X{t) = 1 <;=^ w{t) = exp(— t) whose Laplace 
transform is 'w{p) = 1/(1 + p). Thus, from (fHTl) . G{p) = l/p and G{t) = X{t) = 1 as 
expected. 



3.9 Event-rate spectrum 

The photo-current spectrum e'^^x>{^), where Q denotes the Fourier frequency, may be 
measured by letting the photo-current j{t) flow through a narrow-band filter with center 
frequency Jl. The function yx>{^) exhibits a 2ttD'^S{U) singularity, but the reduced spec- 
trum o5^a_d(^) is non-singular, where AD = V — D. In the theory to be subsequently 
presented, S^/\d{^) is obtained by setting in rate equations d/dt jQ. 
The relative noise is then defined as 

Af{n) ^ yAD/om - ^, (62) 

where ^AD/oi^) = ^ad(0)/D2. 

One may be interested instead in the normalized second-order correlation function g{T), 
a non-negative even function of the delay time r. Aside from normalization, g{T) is the 
correlation of 'D{t) with the singularity at r = being removed. Alternatively, Dg{T)dT 
may be defined as the probability that an event occurs between r and t + dr, given that 
an event occurred at t = r. As shown in Section 13.131 g{T) is related to the relative noise 
defined above by the integral relations 

/oo 
dT(5(r) - 1) exp(-jf^r) (63) 
-oo 

g{T) - 1 = y_ ^■^(^) expOOr). (64) 

The motivation for introducing g{T) — 1 in (f63]l is that this quantity tends to as r tends 
to infinity because widely separated events are in that limit independent for stationary 
processes. The above relations are closely related to the Wiener-Khintchine relations. They 
are established in Section [3.131 directly for point processes. Note that our definition of "sub- 
Poissonian" photo-currents is that AA(0) < 0. This does not necessarily imply that g(0) < 1. 
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3.10 Photo-count variance 

As shown in Section [3.131 the normalized variance V(T) of the number of events occuring 
during some time T and g{T) are related 

- .) . (ee, 

In the special case of a Poisson process we have g(r) = 1, V(T) = and AA(i7) = 0, that is, 
^AD = D. 

As an example let us consider a high-power laser driven by a non-fluctuating current. 
The relative noise will be obtained in the form 

M{n) = ^- ^ (67) 

D{l+{^r,f) 

where Tp is the so-called "photon life time" of the resonator. Prom this expression we obtain, 
setting Z) = 1 for simplicity, that 



+00 



dn 1 , „ ^ 1 , T 



9{r) = 1 - / — -2 exp(jl7r) = 1 - — exp(--) (68) 

1 



5(0) = 1 - ^ (69) 

iTp 

Of course g{T) ^ 1 if r ^ oo. In the present situation ^(0) < 1. 



3.11 Dark-room picture 

For the sake of illustration let us present a simple picture of regular point processes. The 
initial point process considered is periodic and consists of events occurring at t = 1,2... 
time units, that is tk = k. Under circumstances to be defined later on (delay times much 
larger than unity) this process may be viewed as being almost stationary. The density is 
clearly unity. 

In our picture, one person (representing an electron) enters into a dark room every time 
unit and wanders randomly in the room until he finds the exit. This picture may describe 
regularly-pumped lasers at high power because electrons entering the cavity then are quickly 
converted into photons. Photons wander in the optical resonator for some time and then 
get instantly converted into photo-electrons. The point process is written as = k + ^k, 
where the are independent of one-another and distributed according to the same density 
P{^k) = Pic)- An appropriate distribution would be the exponential one. 
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Let us treat a special case that may be solved almost by inspection, namely the case 
where -P(C) = I/t",- if < < and otherwise, and r, are large integers. Consider 
a pair i,j ^ i of k values such that i + may be in the first time slot (0, dt) and j + 
may be in the second time slot (r, r + dr). Inspection shows that this is possible only if 
—Tr <i<0,T — Tr<j<T. Ignoring first the restriction j 7^ i, we find that the probability 
we are looking for is the number of allowed i,j values, that is, the product of the i,j ranges, 
times namely r^/r^ = 1. This result is accurate if r > r,.. But if r < r,. one must 

subtract from the numerator of the previous expression the number of «, j-values that are 
equal, namely — r, so that the normalized correlation reads 

9{t) = 1, T>Tr 

9{r) = "^~^"r"^ =1-^, T<Tr. (70) 

Tr Tr 

In particular, g{0) = 1 — 1/Tr, indicating a modest amount of anti-bunching, remembering 
that Tr ^ 1. The same result is obtained for the laser model in (f69l) if we set = 2rp to 
make the average life-times the same in the two models. 

The reduced photo-events spectrum is obtained from ^(t) through a Fourier transform 
according to ([631) as 



J\f{n) = 2 [ dx{x-l) COs{nTrX) = 2 

Jo 



cos($7Tr) — 1 

2 



{QTrY 



^,^s COs(VtTr) — 1 , , 

where we have set x = r/r^, remembering that the density (average rate) D = 1. We note 
that =5^(0) = 0, as one expects from the fact that the primary process is regular and that 
no event has been lost or created. The spectral density of the process considered, given in 
(TTT]) . is illustrated in FigM 

Using (f65]l we obtain in the present model 

V(T) = -1 + ^, T>Tr 



V(r) = -- + ^, T<Tr. (72) 

Tr oT„ 



It is easy to see that the expression of g{T) in ([661) is verified in that example. 



3.12 Random deletion 

Random deletion of events (also called "thinning") means that each event is ascribed a 
probability 1 — of being deleted. For example, considering the first event of a given 
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Figure 3: Plain line: spectral density corresponding to the dark room picture, see (fTT]) . The 
dotted line corresponds to the shot-noise level. 



run, we flip a coin. If head, that event is preserved (probability 1/2). If tail, it is deleted. 
The same procedure is applied to the other events of the run and to the events of other 
runs, each time with a new coin flipping. Obviously the average rate D of the process is 
multiplied by pr. An important result is that the function g{T) and thus the other two 
functions defined above, and in particular the relative noise AA(0), are not affected. Indeed 
consider the case where there is one event in the time slot [0, dt] and one event in the time 
slot [r, r + dr], corresponding to a product of 1. In any other circumstances the product 
is 0. After thinning the probability of having again (1,1) is multiplied by pr'^. But the 
denominator in the normalized correlation ^(t) is also multiplied by pr"^, so that the result 
is unchanged. The average rate may be restored by an appropriate scaling of the time axis. 
But since in general g{T/pr) ^ ^(t), rescaled thinning affects the statistics with the sole 
exception of Poissonian processes, in which case ^(t) = 1. 

3.13 Point processes. Mathematical details 

Point processes are possibly unlimited sequences of increasing positive real numbers tk,k = 
... — 1, 0, 1, 2 . . . , collectively denoted as {ifc}. Each k value corresponds to an "event", or 
"point", occurring at time tfc. We consider M such sequences, labeled by m = 1,2...M.., 
called runs. Averages denoted by the sign (.) refer to sums from m = 1 to m = M of some 
quantity defined for each run divided by M, letting M go to infinity. The full specification 
of a point process requires correlations of all order. However we limit ourselves to first 
and second-order quantities. Note that the full time axis is considered, with past being 
distinguished from futurj^. 

^"For random points on a line it is often required to treat the two directions "left to right" and "right to 
left" symmetrically. In that case the complete intensity function has to be abandoned. 



58 



Point processes are defined by a complete intensity function 

p{t; Ht) = lim 5-^pr{N{t, t + 6) > Q\Ht}, (73) 

(5— >0+ 

where pr{N(t,t + 5) > 0\Ht} denotes the probabihty that at least one event occurs within 
the interval {t,t + 6], given the previous history Ht, that is the point process up to and 
including t (signs "(,)" indicate that the end point is not included while signs "[,]" indicate 
that the end point is included). The condition that the tk are increasing may be written 
as pr{N{t,t + S) > 1 = o{6)} for all t. This ensures that there are essentially no multiple 
occurrences at a given time. Such processes are called "orderly" for an obvious reason. In 
the special case of Poisson processes the previous history does not matter, and we only 
need to know the density of the process. Other important special cases are processes with 
independent (perhaps equally distributed) increments). 

Given a point process, new point processes may be obtained through the four following 
operations, whose use in the present paper will be indicated. 

• A change of time scale. In that way an inhomogeneous Poisson process may be 
converted into a (homogeneous) Poisson process. 

• Thinning, in which some of the points in the original process are deleted with a 
constant probability, independently of all the other points. This happens when the 
light incident on a photo-detector is attenuated. 

• Translation of individual points. See the section [3.111 

• Superposition, in which a number of separate processes are merged. Given the point 
process relating to a single electron submitted to a field, we shall consider the process 
obtained for independent electrons. Another example is that of a radio-active source 
involving many independent atoms. We assume that superposition of an arbitrary 
large number of orderly processes is a Poisson process [38] • O 

^^As a proof of this conclusion up to second-order, consider the superposition of M independent processes 
of density 1 (for simplicity) and the same second-order correlation g{r). These processes, of duration r, 
may be cut out from a single sample provided they are sufRciently far apart to be almost independent. 
As before, for the original process, (/(r)dr denotes the probability that, given that an event occurred at 
t = 0, another event occurs in the interval (r, r -I- rfr). Consider now the superposed process consisting of 
the k = 1,2...M independent processes. The density of the superposed process is obviously M. Suppose 
that for fc = 1 an event occurs at f = 0. The probability that an event of that same sample occurs in the 
(r, r -I- dr) interval is, as said above, g(T)dT. But events in the (r, r + dr) interval may originate from the 
M — 1 other processes, with probability (M — l)dT. It follows that, given that an event occurred at t = for 
the sample k = 1, the probability that another event occurs in the (r, r + dr) interval is ((?(t) + M — l)dr. 
What has been just said for = 1 applies to all the M samples, so that the probability that, given that an 
event occurred at t = 0, another event occurs in the (r, r + dr) interval is M{g{T) + M — l)dr. Dividing 
by for normalization, we obtain that (?**^^(t) = 1 + ^^^^"^ , which tends to unity in the limit that M 
goes to infinity. Remember that the result (;(r) — 1 for all r values characterizes Poisson processes up to 
the second order. It follows that if we consider a large collection of electrons submitted to the same optical 
field, the detection events tend to be Poisson distributed independently of the single-electron response. 
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Let d{t) be the number of events occurring up to time t, that is the number of k values 
such that tk < t. Obviously d{0) = since the are positive numbers. Let us prove that, 
for some measurement time T, {d{T)) = DT, where D is a constant called the intensity of 
the process. We introduce the (positive) number Dh(t) := d{t + h) — d{t) of events occurring 
between t and t + h. Because the process considered is stationary [Dhit)) does not depend 
on t. It is convenient to split the measurement time T into time slots of duration h = T/n, 
labeled by i = 1,2, ...n. Eventually, we let n go to infinity, so that it is unlikely that more 
than one event occur within any time slot. Thus, if = D]^{{i — \)h) denotes the number 
of events occurring during slot i = 1, 2...n, we have either = 1 or = and [Di) does 
not depend on i. For later use note that Di^ = Di. The number d{T) of events occuring 
during the measurement time T is the sum of the Di with i running from 1 to n, so that 
its average reads 



where we have set D = {Di) /h. 

Because the process is stationary its auto-correlation {Dfi{t + T)Dfi{t)) does not depend 
on t for every h > and every r > 0. The degree of second order coherence ^(r) is the 
limit of {Dh(t + T)Dh{t)) / {Dh{t))^ as h goes to 0. Let us set for j > i 




(74) 



{DiD^) = {Difgn{{j~i)-) 



(75) 



and evaluate 




n(A)+2(A)'i; 9n{{j-i)-) 



1=1 j=i+l 



{d{T))+2{D,fY,in-i)9n{-) 



1=1 
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(76) 



i=l 



In the limit n 



oo the sum may be replaced by an integral and gn by g, thus 




(77) 
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After slight rearranging the variance of d{T) may be written in the form 

^ var(d(r)) _ ^ {d{Tf) - {d{T)f _ 
^ '~ {d{T)) {d{T)) 



2D I dT{l--){g{r)-l) (78) 







since dT{l — t/T) = T/2. The motivation for introducing ^(r) — 1 in the integral is 
that this quantity usually tends to quickly as r tends to infinity. Intuitively, this is 
because widely separated events tend to be independent and consequently in that limit 
{DiDj) « (A) {Dj) = (A) . Setting L> = 1 for brevity, relation ^ may be written as, 

k=0 

where Pc{T)dtdT denotes the probability density of having an event between and dt and 
an event between r and r + dr or, equivalently, Pc{T)dT is the probability density of another 
event being registered during the time interval r and t + dr, given that an event occurred 
at i = 0. In ([79]) P{k, r) denotes the probability of k events being registered between t = 
and t = T. 

We define the event rate 'D{t) = D + AD{t) as the sum over k of S{t — t^), where 6{.) 
denotes the Dirac distribution. Going back to time slots of small duration h, the event 
rate is the sum over i of Di/h where, as before, A = or 1. The denominator h may be 
omitted because of the subsequent normalization. The calculations given below parallel the 
ones given above in relation with the photo-count variance and some details will therefore 
be omitted. The spectral density is defined in terms of the event times tk occurring during 
runs of duration T, to be later tend to infinity, see (f52]l . 



TyADm = {^Di exp{-jni) Dj exp(jOj) \ (80) 
\j=i j=i I 



rp n n rp 

[D})+D\-fY. E 9n{{3-^)-)co^{{3-m) (81) 



n ' 

" n' ^ ' " ' n' 

1=1 j=i+i 

= DT + D^- V(r - —)gnC—) cos(if^). (82) 
n ^-^ n n 

1=1 

The above expression may be transformed as was done earlier for evaluating the variance 
of d. On account of the fact that 

T 

Ei , ,2mn^ , , 

(l--)cos(^) + l = 0, (83) 

i=l 
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where n denotes any non-zero integer, converting the sum into an integral, and in the large 
T limit, we obtain 



This relation between the relative noise AA(r2) and the normalized correlation function 
g{T) has been established directly for point processes. It is however instructive to show 
how this relation may be alternatively derived from the Wiener-Khintchine (WK) theorem 
in (f52]) . We consider the event-rate process 'D{t) = 6{t — tk) which, in the limit /i — > 
corresponds to Dh{t)/h, where Dh{t) is the number of events between t and t + h, as defined 
earlier. For simplicity we suppose, without loss of generality, that the average rate D is 
unity. If R{t) denotes the correlation of T>{t)^ the WK theorem tells us that 



But R{t) presents a singularity at r = 0. We are thus led to define (?(t) = R{t) — 5{t). On 
the other hand, S^-r>{^) presents a singularity 2tt5{^) but ^ad{^) = o5^D(f^) — 27r(5(0) is free 
from singularity. Substituting these expressions in the above WK equation and using some 
form of the (5-distribution we recover (|84|) . Since presently the correlation -R(O) = oo, the 
restriction — -R(O) < R{t) < R{0) established in Section [2. 141 does not entail any restriction 
on g{T). 

3.14 Useful integrals 

A number of integrals from x = — oo to x = oo will be needed in subsequent papers, which 
may be evaluated by contour integration. The method is as follows. 

Complex numbers are denoted hy z = z' + \z" , where = —1. The complex conjugate 
of z is denoted z* = z' — \z". Let f{z) be a function of z whose only singularities are simple 
poles at zi, Z2 ■ ■ ■■ One calls residue at z^ the coefficient of {z — Zk)~^ in the (Laurent) series 
expansion of f{z) near Zk- The integral of f{z) along a closed counterclockwise contour is 
equal to 27ri times the sum of the enclosed pole residues. 

For example, closing the real axis by an upper half-circle of infinite radius we obtain 



Here we have a single enclosed pole at a; = i. The coefficient of l/(x — i) in the integrand 
is l/(2i) when x = i. We obtain similarly 




(84) 




(85) 




(86) 




(87) 
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where a denotes a non-zero constant. 

Further, for apphcation to inhomogeneously-broadened lasers, let us define a weight 
function 

that reduces to the Dirac (5-distribution when g tends to 1, and 

/.„„ = 8b' + ,r f '^r'.M^r (89) 



We obtain 



1 + 



ho = 8g (90) 

/i2 = + 8g{g - 1) (91) 

/20 = V(5-1)+ 4/(5 + 1) (92) 

hi = 8gy (93) 

ho = 3(5 - 1)/ + 65(5' - + 5' (35' + 35 + 2) (94) 

hi=2y[y\g-l) + g\g + 3)] (95) 

/32 = (5 - 1)/ + 25(5' + 3)/ + gHa - 1)(5 + 2) (96) 

h3 = 2y[{3g + l)y^ + 3g\g-l)] (97) 

/34 = (35 + 5)/ + 6g{g^ - l)y^ + g^Sg - 2){g - 1) (98) 



We also need for evaluating the gain of semiconductors 

da; y/x 1 



vr Jo {x + l){x-a) 1 + 



(99) 



if a is negative, and 1/(1 + a) if a is positive. In the latter case, the integral is understood 
in principal value. 

3.15 Bi-complex representation of signals 

As recalled in Section [2]T4] electrical engineers usually factor out a term exp(jrit) to represent 
time-harmonic sources. = 2iTf is the angular baseband (or "Fourier") angular frequency. 
The real signal is obtained by taking the real part of the product Vex.p{jnt), where V 
denotes some complex number. Only time-invariant linear causal systems are presently 
considered. It follows that the system response has the same form as the applied source. 
In Physics, it is usual to factor out a term of the form exp(-ia;t) where uj = 2t:v denotes 
the carrier angular frequency. 
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When a source at frequency uj is modulated at frequency the bi-complex representa- 
tion described below proves useful. To avoid bothering with minus signs, it is convenient 
to set ii =-i and i2 =]■ Further, we set pi = iito and p2 = «2^- The algebra of bi-complex 
numbers is associative and commutative. A bi-complex number is written as 

\; = a + bii + ci2 + diii2, (100) 

where if = 12 = —1, ^1^2 = ^2^1 = j, 3^ = 1, and a, 6, c, d are real numbers. The algebra of 
bi-complex numbers was discovered by Segre in 1892. A modern account of the bi-complex 
numbers algebra may be found, for example, in [39]. A bi-complex number is invertible if 
a}' ^ (? ^ £■ / ±2 [ad — 6c) . To prove it, multiply a^hi\^ 0x2 + dj by a — hi\ — 012 + dj 
and obtain A + Bj^ where A = a'^ + + + d'^ and B = 2{ad — bc) are real numbers. Next, 
note that [A + Bj) {^A — Bj) = A^ — B^ is real. The condition for a bi-complex number to 
be invertible is therefore that A / ztS, which is the above condition. 

The real signal v{t) is equal to V(pi,p2) exp ((pi -1-^2)*) +ccc, where "ccc" means that 
one must add 3 terms to the one written out, one with pi changed to — pi, the second with 
P2 changed to —p2, and the third with both pi and p2 changed to —pi and —p2- Let us 
now state the basic theorem. If v{t) denotes a modulated voltage represented by V(pi,p2) 
as said above, Y{p) denotes the usual complex circuit admittance (the ratio of two real 
polynomials in p, and i{t) the real electrical current flowing through the circuit, we have 

i[t) =Y {pi + p2)V{pi,p2)e^v[{Pi + P2)t] + ccc. (101) 
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4 Classical Circuit Theory 



We recall in the present section the basic concepts employed in Circuit Theory, and describe 
devices useful in radio, microwave, and optical frequency ranges, considering both resonators 
and transmission lines. The circuit equations provide currents i{t), viewed as (real, linear, 
causal) responses to specified potential sources v{t). 

We first consider conservative elements (that is, elements that conserve energy) such as 
capacitances, C, and inductances, L. We next consider non-conservative elements such as 
positive or negative conductances. In most practical cases the conductances are obtained 
from separate measurements. We postpone to Section [6^31 a microscopic conductance model, 
namely a single electron located between two parallel conducting plates. A potential v{t) 
oscillating at an optical frequency u is applied to the plates. We suppose exact resonance 
between v{t) an the electron natural oscillatory motion at frequency cOq- The (Quantum- 
Mechanical average) induced current i{t) is proportional to the average electron momentum, 
and, under circumstances to be discussed, the ratio {i{t)) /v{t) may be a real constant G. 
When eU is slightly smaller than huo, the conductance is positive and the optical potential 
delivers energy to the static potential. On the other hand, when eU slightly exceeds fvjJo 
the conductance is negative and the optical field receives energy from the static potential. 
When eU is precisely equal to hujo the conductance vanishes but fluctuations remain. To 
summarize, the complete system including static and optical potentials conserve energy 
aside from an irreversible loss of energy \eU — hijJo\ that can be made as small as one 
wishes. 

Accordingly, the circuit theory may be applied to non-conservative elements involving, 
besides C and L, conductances G. We mainly consider sources (and responses) that vary 
sinusoidally in the course of time at frequency uj = oJq- For a closed linear system, the 
equations have solutions only for discrete complex values a;„ of w. We will be particularly 
interested in circuits that have only one nearly-real frequency, the other ones having large 
negative imaginary parts, corresponding to strongly damped modes. 

We first evaluate the current i{t) induced in a potential source v{t) using the Classical 
Theory of Electron Motion. This is a deterministic problem and no fluctuations are involved 
if we ignore the random electron motion due to heat and the discreteness of the electrons. 
In the realm of Quantum Mechanics we may evaluate the average value and the higher 
moments of the current, but the function i{t) itself, for a single sample of the ensemble, is 
undefined. 

4.1 Classical devices 

Let us recall basic results. The complex notation often employed for describing quantities 
that vary sinusoidally in time is recalled in Section [2.141 For strictly sinusoidal potentials 
and currents represented by the complex numbers V and /, respectively, and linear circuits, 
we have the generalized Ohm law / = Y{uj)V , where the complex constant of proportionality 
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Y(uj), called the admittance, and its inverse the impedance Z(uj), in general depend on the 
frequency w, which may vary from minus to plus infinity. 

If a potential v(t) is applied to a conductance G (a real number), we have by definition 
i{t) = Gv{t), or, using the complex notation, I = GV, where V, and thus /, are in general 
complex numbers. We consider in the major part of this paper ideal conductances defined 
as follows: They are supposed to be independent of the driving potential V and to be 
independent of frequency. Furthermore, they are supposed to have a fixed energy content 
that may be set equal to zero since only energy differences are relevant. A physical model for 
ideal conductances is a piece of metal having a large number of inelastic scattering centers. 
Electrons accelerated by the applied field quickly loose their energy, which is converted into 
heat. Under such circumstances the electron kinetic energy remains negligible, and thus the 
total energy is fixed. In contradistinction, the input conductance of a loss-less transmission 
line of characteristic conductance Gc terminated by an ideal conductance G = Gc (matched 
load) is equal to G at any frequency. But there is in that case a stored energy equal to 
G|F|^r, where r denotes the transit time of a pulse along the transmission line (this is 
power divided by the group velocity times the line length). Thus, a matched transmission 
line does not constitute an ideal conductance in the sense defined above, even though the 
input conductance G is a real constant. 

4.2 Capacitances 

A capacitance G = CoA/dc may consist of two parallel perfectly-conducting plates of area 
A separated by a distance dc <C VA. The constant Cq in this formula is called the free- 
space permittivity. Its exact numerical value is given at the end of Section [2l There are 
two wires connected respectively to the upper and lower plates, so that electrical charges 
may be introduced or removed. If an electrical charge q is displaced from the (say, lower) 
plate to the upper plate a potential v appears between the two plates given hy v = q/C. 
The energy stored in the capacitance is Ec = Gv^/2 = g^/(2C), a result obtained by 
considering elementary charges dq being displaced from the lower to the upper plate of the 
initially-uncharged capacitance until a final charge q is reached. If g is a function of time 
and C is kept constant, we have v{t) = q{t)/C. We may set q = —Ne, where the number 

of electrons is supposed to be so large that q varies almost continuously. As before, e 
denotes the absolute value of the electron charge. 

Let now v{t) be of the sinusoidal form given above. Because the current i{t) represents 
a fiow of electrical charges into one plate or flowing out from the other plate, i{t) is the 
time-derivative of the electrical charge: i{t) = dq{t)/dt. The relation between the complex 
current I and the complex potential V, as defined above, thus reads I = —\GujV. The 
admittance is in the present case Y{(jj) = —xGuj. If we set y = G -|- ii?, we have therefore 
for an ideal capacitance G = and B = —Cuj. The stored energy averaged over a period 
27r/a;, reads {Ec) = C\Vy /2. In following paragraphs, we will consider a light-emitting 
device driven by a very large capacitance (instead of, say, a battery) with a very large initial 
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Figure 4: a) There is a linear relation between the potential V and the current / at some 
frequency oj (generalized Ohm law / = Y(uj)V). b) Illustrates the KirchhofT law: Ii + 
I2 + Id, = 0. c) Represents a potential source with V independent of I. d) Represents 
a current source with / independent of V . e) Inductance-capacitance resonating circuit, 
f) Cavity employed, e.g., in reflex klystrons, g) Fabry-Perot-type optical resonator with 
two curved mirrors facing each other, h) Whispering-gallery mode resonator, i) Low-pass 
filter, j) Parallel conductors, k) Waveguide, 1) Optical fiber. Waves may be split in various 
ways: m) A transmission line is connected to two transmission lines whose characteristic 
conductances sum up to the original line characteristic conductance, n) Directional coupler. 
The two holes are spaced a quarter of a wavelength apart, o) The beam splitter is an optical 
equivalent of the directional coupler, p) The circulator is ideally a loss-less non-reciprocal 
device. 



charge q such that the potential U = q/C across the capacitance has the desired value, for 
example 1 volt. If the light-emitter operation duration is denoted by T, the capacitance 
supplies a current i during that time, and thus loose a charge Aq = iT. Because q is very 
large we have /S.q <C q provided the experiment does not last too long. As a consequence the 
potential U across the capacitance does not vary appreciably. We realize in that manner a 
constant-potential source, that is a source whose potential does not depend appreciably on 
the delivered current. 

A constant-potential source at optical frequency uj may be realized in a similar manner. 
Again the capacitance C and the initial charge q are supposed to be arbitrarily large, but 
we now allow the spacing dc between the capacitance plates to fluctuat^^l at the optical 
frequency uj. This spacing variation entails a fluctuation of the capacitance, and thus of 
the potential across the capacitance since the charge is nearly constant as was discussed 
above. The potential across the capacitance may be written as ?7 + v{t). The important 
point is that the optical potential v{t) as well as the static potential U are independent of 
the current delivered. That is, if atoms are present between the two capacitance plates, 
processes occurring in the atomic collection have no influence on the field. We have just 
described an essential component of our circuit-theory schematic. In contradistinction, the 
potential across a resonating inductance-capacitance circuit modeling a single-mode cavity 
does depend on atomic processes. For that resonator configuration the assumption that the 
optical field is nearly constant holds only in the large potential (or large photon number) 
limit. 

4.3 Inductances 

An inductance L may be constructed from a cylinder of area A and height cLl ^ split 
along its hight, so that an electrical current may flow along the cylinder perimeter. In that 
case L = fiA/di, where fj, denotes the permeability. One may assume that the cylinder 
contains electrons that have magnetic moments. Just above the Curie temperature, 
much exceeds the free-space permeability, so that the latter may be set equal to zero, as 
was discussed earlier. For an inductance L, the magnetic flux (or magnetic charge) is 
(p{t) = Li{t), and the potential across the inductance is v{t) = d(j)/dt. It follows that 
for a constant L, v{t) = Ldi{t)/dt, or, using the complex notation V = —xLujI. Thus 
Y{uj) = I IV = \/[Luj). The energy stored in an inductance with a current i flowing 
through it is El = Li^ /2. For a sinusoidal current represented by the complex number 
I, the time-averaged energy is {El) = L /2. If a large inductance supports a large 
magnetic flux, the current flowing through the inductance is nearly independent of the 
potential across the inductance. In that manner, we may realize constant-current sources, 
either static or oscillating at optical frequencies through a change of L (e.g., by changing 

^^Practically- minded readers may object that mechanical motion may not be feasible at high frequencies. 
Let us recall here that the numerical value of uj is arbitrary. What we call "optical" frequency lo/2tt may 
be as low as IHz provided that the other frequencies considered be much lower, e.g., 1 mHz. 
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the coil length if a coil instead of a simple cylinder is employed) . 

The linear relationships outlined in previous paragraphs are sometimes referred to as the 
"generalized Ohm laws". Let us recall that there is a well-known duality between potentials 
and currents and between electrical charges (expressed in coulombs) and magnetic fluxes 
(expressed in webers), so that expressions obtained for capacitances may be translated into 
expressions relating to inductances. 



4.4 Energy and power 

As an application of the above energy formulas, let us consider a circuit consisting of an 
inductance L and a capacitance C connected in parallel. Since the system is isolated the 
total admittance must vanish and we obtain the resonance formula LCul = 1, where ujo 
denotes the resonant frequency. The sum E of the energy EL{t) located in the inductance 
and the energy Ec{t) located in the capacitance does not vary in the course of time. This 
is twice the time-average energy stored in the capacitance (or inductance). Using above 
formulas we find that the rms (root-mean-square) field across the capacitance is 



E 



(102) 



where V = Adc denotes the capacitance volume. We later show that when a resonator 
such as the one presently considered is in a cold environment it eventually reaches a state 
corresponding to an energy ^o/2, where h denotes the Planck constant (divided by 27r). 

According to the above formula, the so-called " vacuum (rms) field" reads f^jocuum — \/ ^ °y . 
The two oppositely-charged capacitance plates attract one another with an average force 
F = d{hiOo/2)/d{dc) = huJo/i^d,). 

If two sub-systems are connected to one another by two perfectly conducting wires with 
a potential v{t) across them and a current i{t) flowing into one of them (the current —i{t) 
flowing in the other one), the power flowing from one sub-system to the other at some 
instant t is equal to v{t)i{t). For sinusoidal time- variations, the power averaged over an 
oscillation period reads P = ^{VI*}. 

Finally, let us recall that at a node, that is, at the junction between perfectly conducting 
wires, the sum of the currents entering into the node vanishes as a consequence of the fact 
that the electric charge is a conserved quantity. For three wires traversed by currents 
h{t), ^2(^)5 hit), for example, we have at any instant ii{t) + i2(t) + 13(1) = 0. It follows 
that the complex currents sum up to zero, that is Ii + I2 + I3 = 0. Both the real and 
the imaginary parts of the sum vanish. Such relations are sometimes called "generalized 
Kirchhoff laws". The above discussion suffices to treat circuits consisting of conductances, 
capacitances and inductances arbitrarily connected to one another. Some circuits require a 
more complicated description involving for example (non-reciprocal) gyrators. These latter 
components are useful to separate reflected and incident waves. 
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4.5 The tuned circuit 



For the sake of illustration and later use, let us generalize the resonator previously consid- 
ered by introducing in parallel with the capacitance C and the inductance L a conductance 
G. The relation between a complex current source C at frequency u, supposed to be inde- 
pendent of frequency, and the potential V across the circuit reads 

The power dissipated in the conductance G at frequency uj reads 

2 

P{.) =G\V{u)\^ ^ ^ (104) 

in the small-loss approximation. Thus -P(c4;) drops by a factor of 2 from its peak value when 
2C{u}± — ujo) = ±G. The full- width at half power (FWHP) 5uj of the resonance that is, the 
difference of (angular) frequencies at which the dissipated power drops by a factor of two, 
is 

G 1 

(5a; = a;+ — = — = — , (105) 

C Tp 

where Tp = C/G is sometimes called the "photon lifetime". If the resonator is left alone in 
a cold environment (T^=OK) its classical energy decays according to an exp(— t/r^) law. 
For a Fabry-Porot resonator with mirrors of small power transmissions Ti, T2, respectively, 
and spacing L, we have 

where v denotes the group velocity and 2L/v is the round-trip time. 

The energy contained in the resonator is twice the average energy contained in the 
capacitance whose expression was given earlier. We then obtain in the small-loss approxi- 
mation 

where x = 2Tp (u — Uq) ■ 

4.6 Derivative of an admittance with respect to frequency 

For late use, note the expression of the derivative with respect to co of the admittance of a 
linear circuit, submitted to a voltage V 

(108) 
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where the sum is over all the circuit capacitances and inductances. Vk denotes the (com- 
plex) voltage across the capacitance Cfc and Ik the (complex) current flowing through the 
inductance L^. The circuit resistances or conductances do not enter in the sum. This 
relation is readily verified for an inductance in series with a resistance and a capacitance 
in parallel with a conductance. Thus the relation holds for any combination of elements, 
connected in series and in parallel. 

4.7 Matrix formulation 

For an arbitrary circuit, the task is to "extract", figuratively speaking, the (positive or 
negative) conductances from the given circuit, each conductance being connected to the 
conservative circuit that remains after extraction of the conductances. If (positive or 
negative) conductances are involved, the circuit becomes an A^-port conservative device. 
For an A''-port circuit, we define the vectors V = [Vi,V2, ..-VnY and / = [Ii, I2, ...InY 
where the upper t denotes transposition. The linear relation is written in matrix form 
I = Y{uj)V, where Y{uj) is called the circuit admittance matrix. For a conservative circuit 
the total entering power ^{V^*I} = 0. Since this relation must hold for any source this 
imphes that 1"** + 1" = 0. 

It is convenient to view the connections between the conservative circuit and the con- 
ductances as ideal transmission lines of small length and characteristic conductances Gc- 
Supposing that Gc = 1, the potential V across one of the transmission lines and the current 
I flowing through the (say, upper) wire, are combined into an ingoing wave whose ampli- 
tude is defined as a = V + I and an outgoing wave defined as b = V — I. Since under 
our assumptions the circuit elements are linear, there is a linear relationship between the 
a-waves and the 6-waves. The relation between b and a, defined like I and V above, may 
be written in matrix form as 6 = Sa, where the S matrix is called the circuit "scattering 
matrix". Because the circuit is conservative, the outgoing power equals the ingoing power. 
It follows that the S-matrix is unitary, i.e., S**S = 1. We need not assume that the circuit 
is reciprocal, however, that is, the S-matrix needs not be symmetrical. 

4.8 Various circuits 

We have represented a number of important conservative (loss-less, gain-less) components 
in either their circuit form, their microwave form, or their optical form in Fig.[4l The origin 
of the differences is that, as one goes to shorter wavelength (higher frequencies) some circuit 
elements become too small to be fabricated. It should also be noted that metals, such as 
copper, that are excellent electrical conductors up to microwave wavelengths, do not behave 
as electrical conductors any more at optical wavelengths because of electron inertia. On the 
other hand, while it is difficult to find very low-loss dielectrics at microwave frequencies, 
extremely low-loss glasses exist at optical frequencies. Fig. |4] represents four resonating 
circuits, that one may call "0-dimensional" devices. Namely, the inductance-capacitance 
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circuit employed up to about 100 MHz, the cavity employed in reflex klystrons and masers 
for example, the Fabry-Perot resonator consisting of two mirrors facing each other, and the 
whispering-gallery-mode dielectric resonator, first demonstrated in the microwave range 
and now-a-days employed in the optical range. Resonators are primarily characterized by 
their resonant frequency ujo- Small losses may be characterized by the so-called "photon 
life-time" Tp defined earlier. When the resonator size is large compared with wavelength 
many resonating modes may be present. In most applications it is desirable that only one 
of them be loss-less, or nearly so (see, e.g. |40|). 

FigurelHrepresents four one-dimensional devices called "transmission lines". The circuit 
form is a periodic sequence of series inductances and parallel capacitances. The microwave 
form consists of two parallel conductors characterized by a characteristic conductance Gc, 
with waves propagating at the speed of light. Above IGHz one would rather use waveguides. 
The optical form is the now-a-day well-known optical fiber. A glass fiber (core) in vacuum 
may guide optical waves by the mechanism of total reflexion. In order to increase the 
core size without having spurious modes propagating, the core is usually immersed into a 
lower-refractive-index glass. 

Other useful devices are shown in Fig. ID The power carried by a transmission line 
of characteristic conductance Gc may be split into two parts simply by connecting it to 
two transmission lines whose characteristic conductances sum up to Gc- This a three-port 
reciprocal conservative device. Alternatively, when two transmission lines are put side by 
side and coupled at two locations separated by a quarter of a wavelength, some of the power 
incident on a transmission line is transmitted into the other one. This device is called a 
directional coupler. This 4-port device may be reduced to a 3-port device by putting a 
matched load at the end of one of the transmission lines. The optical form of a directional 
coupler is called a beam-splitter, which may simply consist of a flat piece of glass. An 
important non-reciprocal 3-port device is the circulator, which exists in microwave and 
optical versions. It is intrinsically loss-less: a wave entering into port 1 entirely exits from 
port 2, a wave entering into port 2 entirely exits from port 3, and a wave entering into port 
3 entirely exits from port 1. Such a device is convenient to separate reflected waves from 
incident waves without introducing losses. 
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5 Electron motion and induced current 



We first consider a single electron located between two parallel conducting plates pierced 
with two holes. A static potential source U is applied to external plates, as shown in Figl5] 
in (a). An alternating potential source v{t) at frequency u are applied between the inner 
plates. These potentials are supposed to be independent of the current that the electron 
motion may induce. The sign convention is the one given in Fig|4] in (c). In both the 
Classical and Quantum Theories, the induced current i{t) is, to within a constant, equal 
to the electron momentum p{t), although the interpretations of i{t) and p{t) differ. The 
power v{t)i{t) supplied by the alternating potential source, once averaged over a period, is 
denoted P{t). 

The over-all effect of the electron motion is to transfer energy from the static source 
to the alternating source (stimulated emission) or the converse (stimulated absorption). 
Spontaneous emission in the usual sense does not occur. What may occur is that the 
electron is emitted or captured by one of the plates (Classical view-point) or tunnels into 
the plates (Quantum view-point). In the Classical treatment, one first evaluate 1) the 
electron motion under the static field, 2) the perturbation caused by alternating field, and 
3) the electron momentum and the induced current. The same steps are taken in the 
Quantum treatment. Namely, we consider the stationary states of the electron submitted 
to the static field, the perturbation of those states due to the alternating field, and finally 
evaluate the induced current from the electron momentum. 

5.1 Classical Equations of Electron Motion 

The equations of motion of an electron of charge — e and mass m are first established for 
the case of a static (time-independent) potential. As an example consider an anode at zero 
potential and an electron emitted from a cathode at potential —U in vacuum, and look 
for the electron motion and the induced current. If i{t) denotes the current delivered by 
the potential source, the power Ui{t) must be equal at any instant to the power delivered 
to the electron, which is the product of the velocity p{t)/m, where p denotes the electron 
momentum, and the force eU/d exerted upon it, where d denotes the electrode spacing. 
Since the potential U drops out from this equation, the current is 

m = (10!)) 

Solving the equations of motion, we find that the diode current i{t) = {^e^U /md'^^t increases 
linearly with time and drops to zero when the electron reaches the anode. Thus, each elec- 
tron freed from the cathode entails a triangularly-shaped current pulse. If i{t) is integrated 
over time from t = to t = r we obtain the absolute value of the electron charge e. We will 
neglect the pulse duration (or transit time) r, so that triangularly-shaped current pulses are 
approximated by e(5(i)-functions. The above theory is applicable only when few electrons 
are emitted so that the initial electron velocities and space-charge effects are neglected. 
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The expression for the current was estabhshed for the case where the potential generates 
a constant electric field. This is not the case for example if the anode and the cathode are 
coaxial cylinders of radii and Tc with Tc < Va. The electron is submitted to a force 
J^{r) inversely proportional to r, assuming radial motion. The induced current i{t) may 
nonetheless be obtained from the p{t) trajectory by the same argument as above. 

As a second example, consider a one-dimensional square-well, whose potential is equal 
to for |a;| < d/2 and infinite (or nearly so) beyond. This potential may be generated by 
parallel anodes at potential and cathodes at potential — C/, as shown in Fig.O The electron 
space-time trajectories x = x{t) consist of straight lines with slopes dx{t)/dt = ±p/m^ 
where p^/2m = E \s the electron energy, which may be selected arbitrarily from to eU 
so that the electron is not captured by the cathodes. The electron is prevented from being 
captured by the anodes by a strong magnetic field in the x-direction. The quick electron 
incursions between anodes and cathodes are here neglected. We may consider in particular 
a lower electron energy Ei and a higher electron energy £"2, corresponding to small and 
large slopes in the x = x{t) diagram, respectively. 

The electron motion induces an electrical current i{t) in the potential source, which is 
proportional to the electron momentum as said above. In the case of a static potential 
source the induced current does not correspond to any power delivered or received by the 
source on the average, so that the electron motion may go on, in principle, for ever. 

If now the static potential U is supplemented by a sinusoidal potential v{t) of small 
amplitude, whose frequency is resonant with the electron motion described above, the 
unperturbed electron momentum p{t) does cause the alternating potential source to receive 
or generate power, depending of the electron state. However, if we consider a large collection 
of unperturbed electrons, the power averages out again to zero. It follows that a net energy 
transfer may be obtained only if we take into account the fact that the alternating potential 
perturbs the electronic motion. In the present classical picture this amounts to bringing 
all the electrons with the appropriate phase with the alternating potential, an effect called 
"bunching". This name originates from the fact that electrons initially spread out uniformly 
eventually are forced by the alternating field to form periodic "bunches". This desirable 
bunching effect is limited by the velocity spread of electrons originating from the hot cathode 
and the fact that electrons tend to repel each others (space-charge effect). 

This is not however the end of the story. The electron, initially in the lower energy state, 
may gain enough energy to be captured by a cathode. Once in a cathode, the electron flows 
through the static potential source to the anode, delivering an energy eU to that source, and 
is emitted by the anode back into the lower energy state. The net effect of these processes 
is that some power is being transferred from the alternating potential source to the static 
potential source, or the converse, the electrons playing an intermediate role. 

In more general situations, the Classical Equations of Motion of electrons of charge — e, 
mass m, and potential energy —eu{x^ t) are best based on the Hamiltonian formulation. 
The particle total energy E{t) is expressed as a function of position x, momentum p, and 



74 



X 



reflector 



-2q 



quiet 
radiation 




anode 



(a) 



catliode 

(b) 



X 






/ Vi(^) 




(c) 



quiet 
radiation 







L 




T 



AIGaAs QWs 
in AlAs 



(d) 

Figure 5: a) Illustrates the potential created by two anodes (inner electrodes) and two 
cathodes located just outside the anodes. The potential is generated by a large, charged, 
capacitance, shown on the left. According to the Classical Picture the space-time electron 
trajectory is almost a zig-zag path, with slight incursions of the electron between the anodes 
and the cathodes, b) Represents a reflex klystron, which is similar to the previous schematic, 
but with a resonator added to it. The current is regulated by a space-charge limited cathode, 
c) Pictures the wave-functions of the ground state and first excited state of a square potential 
well, d) Represents a "surface-emitting" laser diode. The current is regulated by a large 
cold resistance i?,. 
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time t according to the relation 

2 

H{x,p, t) - E{t) = eu(x, t) - Eit) = 0, (110) 

2m 

where p'^/{2m) represents the kinetic energy. The Hamiltonian equations read 

dxjt) _ dH{x,p,t) _ p{t) 

dt dp m 

dp{t) _ dH{x,p,t) _ ^ du{x,t) 

dt dx dx 

The first equation says that the particle momentum p{t) = mdx{t)/dt, and the second 
equation may be written, with the help of the first equation, in the Newtonian form 
m d'^x{t)/dt^ = e {du{x,t)/dx)^^^(j-y Going back to the first example in this section, 
let us consider two parallel plates located at x = and x = d, and at potentials and 
— respectively. We have u{x,t) = —u{t)x/d, and thus the equation of motion is 
md^x{t) / dt^ = —eu{t)/d. The electron decelerates if u{t) > 0. The electron is repelled by 
the negatively-charged cathode. 

As far as static conditions are concerned, an electron submitted to a static poten- 
tial source U is analogous to an electron submitted to the Coulomb potential created by 
positively-charged nuclei. The potentials, on the order of 1 volt (corresponding to potential 
energies of 1.6 10~^^ joules) are comparable in the two situations. In the case of atoms, 
however, the Bohr radius, which is roughly equivalent to our distance d, is on the order of 
0.05 nanometers while, in the case of two conducting plates, the distance can hardly be less 
than 100 nanometers for practical reasons. As a consequence there exist in the two-plate 
model many states whose energy is comprised between the lower-state energy Ei ^ and 
the upper-state energy E2 ~ eU. In both cases the conductance (ratio of the induced cur- 
rent to the applied potential) is initially equal to zero and grows in time linearly until the 
electron is somehow absorbed or leaves the interaction region. 

5.2 Quantum Equations of Motion 

Before entering into the mathematical details, it is important to understand the significance 
of our schematics, and how these schematics may describe actual devices. The differences 
between schematics and real maser or laser devices, which we will point out below, are 
considered to be of minor importance as far as concepts are concerned. 

The configuration that we have in mind is again the one shown in Fig. [5] in (a), but with 
the electron motion quantized as in (c). In this picture, as was discussed earlier, the electron 
is submitted to a static potential source generated by a charged capacitance of arbitrarily 
large value. This potential is applied between inner electrodes and outer electrodes. The 
electron is constrained to move along the x-axis with the help of a magnetic field (not 
shown on the figure). Classically, the electron performs a zig-zag x{t) path. Prom the 
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time-independent Schrodinger-equation view point, the electron may reside only in a lower 
state 1 and an upper state 2. In real masers or lasers the potential is generated by the static 
potential of fixed, positively-charged point-like nuclei (because of their large mass, plates 
or nuclei recoils may be neglected). However, from our view-point, the two configurations 
differ only in the form of the wave-functions and the value of the transition element later 
on denoted by xu. 

The static potential is supplemented by an alternating potential source, at a frequency 
on the order of 10 GHz for klystrons and 300 THz for lasers, generated in the picture of Fig. 
[5] in (a) on the left by a sinusoidal motion of the inner capacitance plate. In a real klystron, 
the alternating potential is generated by an inductance-capacitance circuit, as shown in 
Fig. [5] in (b). In general, this oscillator, resonant with the electron alternating motion, 
cannot be considered as a source, because the potential depends on the induced current. 
It is only in the limit where the tuned-circuit capacitance would be extremely large and 
the inductance extremely small (so that the resonating frequency remains the same), that 
this tuned circuit could be considered as an alternating potential source. Indeed, in the 
limit considered, for a given alternating field, the tuned circuit energy is extremely large 
and little affected by the electron motion. In Quantum Optics, this limiting situation is 
described by saying that "the number of photons in the cavity is supposed to be extremely 
large, so that a classical treatment of the field is adequate". Let us emphasize that the 
configuration treated in the present section is only one idealized component of a complete 
laser device. We need to characterize this component accurately (in terms of conductances 
and event processes) before going on. 

What is missing in the schematics of Fig. [5] in (a) is the absorber of radiation. In that 
schematics this absorber could be realized by adding on the right a triple-plate capacitance, 
as already shown on the left. Similarly, in Fig. [5] in (b) the wavy line, symbolizing the 
escape of radiation, could be replaced by a potential configuration similar to the one shown 
on the left, but with a slightly different static potential, so that power flows from the 
potential source on the left to the potential source on the right. The electron motion and 
the alternating field may be viewed as playing an intermediate role. If this is the case, 
one may wonder why complicated devices are needed to merely transfer energy from one 
capacitance (or battery) to another. The answer of course is that in the microwave or optical 
forms, energy may be carried over large distances with little (absorption or diffraction) loss. 
High-frequency electromagnetic waves also serve as sensors, e.g., in the radar. Transmission 
lines are not shown in Fig. [H 

The Quantum Equations of Motion of an electron of charge — e and mass m are first 
established for a static (i.e., time-independent) potential source. As an example we consider 
a one-dimensional square-well, whose potential is equal to for |x| < d/2 and infinite (or 
nearly so) beyond. This potential may be generated by parallel anodes at potential and 
cathodes at potential — C/, as shown in Fig. [5l We solve the time-independent Schrodinger 
equation and obtain in particular a state 1 with lower energy Ei and a state 2 with higher 
energy E2. As we shall see, these two states correspond to wave- functions ipiix) = cos(x) 
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and tp2{x) = sin(2x), respectively, leaving aside constants. In the case of a static potential 
there is no energy exchange between the potential source and the electron when the electron 
is initially in a stationary state, so that the electron remains in the stationary state, in 
principle, for ever. There are no energy exchange either if we perform a time averaging when 
the electron is in a superposition of stationary states. This situation may be compared to 
the one discussed classically above. 

Let now the static potential source U be supplemented by a sinusoidal potential source 
v{t) of small amplitude, whose frequency is (in some sense to be defined later) resonant 
with the electron motion described above. A net energy transfer may be obtained only if 
we take into account the fact that the alternating potential perturbs the electronic motion. 
In the Classical picture this amounts to bringing all the electrons with the appropriate 
phase, an effect called "bunching", as said previously. In the Quantum picture (time- 
dependent Schrodinger equation), the electron wave function tp{x,t) is the weighted sum 
of the unperturbed states defined above, with time-dependent weights. The theory leads 
to (Rabi) oscillations between the two states. Initially, the induced current is equal to zero 
and grows in proportion to time, but the conductance vanishes on the average. A non-zero 
positive conductance is obtained if the electron initially in the lower state remains in the 
interaction region for a finite time r. We may then evaluate the average conductance "seen" 
by the optical potential source. A model based on the direct coupling between bands of 
states, assuming that an equilibrium is quickly reached within each band separately, would 
be more realistic for semi-conductors. 

The quantum treatment is based on the Schrodinger equation 

[H{x,p,t)- E]i;{x,t) =0, E = ihd/dt, p=-ihd/dx, (113) 
where the sign "9" denotes partial derivation. ip{x,t) is called the wave- function, whose 

2 

initial value 'ip{x,0) is supposed to be known, and H{x,p,t) = ^ — eu{x,t) as in the Clas- 
sical Equations of Motion, but p and E are now operators of derivation. It is easily shown 
that, provided tp^Xjt) decreases sufficiently fast as x ^ ±oo, the integral over all space of 
\^p{x,t)\'^ does not depend on time. It therefore remains equal to 1 if the initial value is 1, a 
result consistent with the Born interpretation of the wave function. Because of linearity the 
sum of two solutions of the Schrodinger equation is a solution of the Schrodinger equation 
(superposition state). The wave-functions add up, but not in general the probabilities. 

Let us evaluate the time derivative of the average value of x. We have, using the above 
Schrodinger equation to obtain d'4j{x,t)/dt and the mathematical relation in ([43]) with 
p = hk 

dixit)) d , ^ ^ 

= -j dxx V(X, t)r{x, t) 

/•+°° , fdip(x,t) , ,d^*(x,t)\ (pU)) , , 
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This is the first Ehrenfest equation. Thus the classical relation p = rn^ still holds provided 
X and p be replaced by their QM-averaged values. 

5.3 Static potentials 

Let us suppose that u{x, t) = u{x) does not depend on time. In that case solutions of the 
above equation of the form tp{x,t) = ipn{x) exp{—i(xint) may be found, where n = 1,2.... 
The ipnix) are real functions of x and En = hiVn that form a complete orthogonal set of 
functions. 

For n = 1,2 the wave functions obey the differential equations 
h'^ d'^Mx) 



2m dx"^ 
d'^^2{x) 



+ eu{x)ipi{x) + Eiipi{x) = 

+ eu{x)i;2{x) + E2MX) = 0. (115) 



2m dx^ 

with the appropriate boundary conditions. They may be ortho-normalized such that 



dx 'lprn{x)lpn{x) = Smn, (116) 

where 6mn = 1 if m = n and otherwise. 
5.4 Potential well 

As an example consider an electron of mass m moving along the x axis be reflected by 
boundaries at x = —d/2 and x = d/2 where the wave-function is required to vanish, that 
is, tp{±d/2) = 0. The lowest-energy state n = 1 and the first excited state n = 2 are 

'ipi{x,t) = \/2/dcos{-Kx/d) exp(— ioJit) (117) 
tlJ2{x,t) = y/2/dsm{2TTx / d) ey.Y>{-\0J2t) (118) 

Notice that i^iix) is even in x, while '4^2{x) is odd in x. Substituting these expressions in 
the Schrodinger equation (I113|) with m(x, t) = 0, we obtain that 

d^ijJx) , , , , , 

^ +hiOnMx) = (119) 



2m dx"^ 
provided 

En = fuvn = ^^n^ n = 1, 2. (120) 

We will see later on that optical fields at frequency uJq = i02 — = (Svr^/i) / (2md^) may 

cause the system to evolve from state 1 to state 2 and back. Numerically, fkOo ~ 1-12 
electron-volt if d = 1 nano-meter. 
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For later use let us evaluate 



rd/2 

xi2 = dx X 'il;i{x)ip2ix) 

J~d/2 

2 f^/^ 

= - / dx X cos{TTx/d)sm{27rx/d) 
d J-d/2 

16d , , 



where we have used the mathematical relation 

/ tcos(t)sin(2t)(it = -. (122) 
y-7r/2 9 

The parameter xy}. determines the strength of the atom-field coupling. It is convenient to 
define a dimensionless oscillator strength 

2muJo 9 256 
/ - -^^4. = ^ « 0.96. (123) 

The maximum value of / is 1. 
5.5 Perturbed motion 

We next suppose that a potential source v{t) = \/2V cos{u>ot) is applied between the two 
anodes in Fig. [5l Since the potential varies linearly with x the electron is submitted to a 
space-independent optical field £{t) = £ocos{ujot),£o = V^V/d, where lOq = ^^2 — t^i is the 
1-2 transition frequency defined in the previous section. In that case (I113p reads 

Hi; = (— - eSo cos(uJot)x - e] Mx, i) = 0, E = ihd/dt, p = -ihd/dx, (124) 
\2m / 

remembering that for stationary states tpn{x) 

^ - hion] Mx) = 0. (125) 

The wave function may be expressed as an infinite sum of '(/'n(a^) exp(— iw^t) functions 
with slowly time- varying coefficients C„(t), that is 

oo 

ij{x, t) = ^ Cn{t) exp(-iu;„t)'0„(x). (126) 

n=l 
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We first evaluate 



2 X oo 

e£o cos{u!ot)x ] ip(x, t) = y Cnit) exp(— iwnt) (hoJn — e£o cos(u!ot)x) i^nix) 

2m I ^ 

' n=l 

Eil,{x,t) = ^ exp(-ia;„t) f ^„C7„(t) + j (127) 

n=l ^ ^ 

If we subtract the first expression from the second and substitute this expression into the 
Schrodinger equation, taking (|125p into account, we obtain 

= ^ exp(— ia;„t) ( \h — ^ h e£oCos{uJot) x Cn{t) 1 il^n{x)- (128) 



n=l ^ ^ 

If we multiply p28|l throughout by V'mC^;), integrate with respect to x, and take into account 
the ortho-normality of the ipmix) functions, we obtain an infinite number of exact ordinary 
differential equations that can be solved numerically. 
Considering only states 1 and 2, we set 

i;{x,t) = Ci(t)exp(-iwit)V'i(x) + C2(t)exp(-i(j2t)V'2(x). (129) 

Introducing the resonance condition cJq = 0^2 — i^i , we obtain from (|128p 

= ih ^^]^^^ + exp(— iwot) cos(iOot)£oexi2C2(t), 
at 

= i^~^^ + exp(— iwot) cos{L0ot)£oexi2Ci{t) , (130) 

where X12 is given in (|12ip . Because the wave-functions ipi{x),ip2ix) are real xu is real, 
and because of the symmetry of the wave-functions xu = X22 = 0. 

The rotating- wave approximation consists of keeping only the slowly- varying terms p5], 
that is, replacing exp(— itJoi) cos(i^oi) by 1/2. Thus, the complex coefficients Ci{t),C2{t) 
obey the differential equations 

d-9M=,^C2it) ^=^^C,(t) C,it)Cm+C2it)Cm = l, (131) 
where <^ lOq is the Rabi frequency given by 

hil^R = £06X12- (132) 
For the potential considered and the value obtained in (I123|) . the above relation reads 

16 

= -^eV2V w 0.17 eV2V. (133) 
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The pair of first-order differential equations in (1159(1 is easily solved. Assuming that the 
electron is initially {t = 0) in the absorbing state, we have the initial condition C2(0) = 0. 
The wave function thus reads 

ij{x,t) = Ci(t) exp(-iu;it)'0i(a;) + C2(t) exp(-iLJ2t)'02(a;) 

= cos(^^t)^i(x) exp(— icjit) + isin(^^t)'i/'2(2;) exp(— iw2i)- (134) 

It follows that 

P22{t) = C2{t)C2{t) = sm{—t) = 

pUt) - nCi{t)Cm} = (135) 
5.6 Momentum probability law 

Let us now evaluate the momentum probability law P{p, t). The wave functions in momen- 
tum space are defined (IMll with the Planck constant restored and at t = 0, 



1 

^^^^^ ^ V2nK J d/2 '^^^^P'^~^y)'^i'^^^ 

^2(P) = / dxexp(-i— )V'2(2;), (136) 

V27r/i J-d/2 ^ 

where tpi{x),tp2{x) are given in pi7|) . Because ipi{x) is real even and ^p2ix) is real odd, 
tpi{p) is real and tp2{p) is imaginary. We set V'2(p) = ^ip2{p)- The explicit result is 



, , , vvrd cos(pd/2h) 
^i(p) = 



2 {pd/2hy - {TT/2y 

sm(pd/2h) 
{pd/2hf - 7r2 



^2(P) = ^TZJlir^^n (137) 



These expressions, however, will not be needed. 

Since the Fourier transform (with respect to x) is a linear operation, we obtain from 
(I134|) the expressions 

iP{p,t) = Ci{t)Mp) exp(-iwit) + C2{t)Mp) exp(-ia;20 (138) 
= cos{-^t)ipi{p) exp(-ia;it) - sm{-^t)tp2ip) exp(-iu;2i), (139) 
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It follows that 




(140) 



where we have used the expressions in (|134l) of Ci{t), C2{t), and taken into account the fact 
that ipi{p),ip2{p) are real. Of course, the integral of P{p, t) over all p- values is unity at any 
time. 

Without electron-optical field coupling, that is, when 17/? = 0, we have P{p,t) = ipi{p)'^ 
and thus {p{t)) = according to the relation below dSTj. Since the average current induced 
in the potential source is proportional to {p{t)), the average current vanishes. It follows 
that for stationary states there are no energy exchange between the (static) potential source 
and the electron on the average. 

5.7 Average induced current 

To evaluate the average momentum when ^Ir > 0, we notice that the first two terms in 
(1140(1 do not contribute. We thus obtain 




dp p ipi{p)ip2{p) = ism{^Rt) cos{uJot)pi2 



(141) 



where, according to with ei = 2mEi/h^,e2 = 2mE2/h^ 




= — V 



m 



dx X'4ji{x)'4j2{x) = —\'mUJoXi2 



(142) 



J-d/2 
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Thus {p{t)) = sm(i^jit) cos{LJot)mLLioXi2. The QM-averaged induced current is Pl. 

(H*)) = = ^o—j- sm[ilRt) cos{uJot), (143) 

It follows that (i(t)) varies, like the optical potential v{t) = \/^V cos{uJot), essentially 
according to cos(a;ot)-law, but with a slowly-varying factor proportional to sin(r2/jt). Since 
the electrons are not directly coupled to one another, the induced averaged current is 
proportional to the total number of electron^. 

The QM-averaged energy E{t) supplied by the optical potential source to the electron 
from time (when the electron energy is equal to zero) to time r is obtained by integrating 
the (QM averaged) power P{t) = {i{t)) v{t). Averaging over the fast optical variations, that 
is, replacing cos^(a;ot) by 1/2, we obtain 

E{T)=hLOoSm\^T) (144) 

where we have employed the expression in (I132|) of the Rabi frequency in terms of xi2- The 
expression in (I144p shows that, if the time r of interaction of the electron with the field is 
such that OijT = 7r,37r..., the energy supplied by the optical potential source up to that 
time to the electron is equal to the energy ell ~ hwo- In general the electron evolves from 
the absorbing state at t = to a state superposition, receiving from the optical source an 
average energy /iWo sin^(^r), and has some non- unity probability of being in the emitting 
state at the exit time r. If Qrt <C 1 this probability is very small, yet measurable. 

The (QM) average conductance "seen" by the (deterministic) alternating potential 
source v{t) = \/^V cos{u)ot) is initially 

_ (^W) e2 Xi2s2 sin(l^Rt) 



^^The expression = ~ 2^^llP^ (using the first Ehrenfest equation) is often expressed in 

a different but equivalent form in terms of the electron dipole moment ex. Omitting for brevity the QM- 
averaging signs and time arguments, and considering A'^ electrons in a volume V = Ad, the current density 
Jopt = ^ = ^^^dT^- One may view the system as a medium of susceptibility x with the polarization related 
to the optical field by Popt = -y^ea; = eoxS- The expression Jopt ~ dPopt/dt coincides with the previous 
one. The current density Jopt, the field £, and the electric induction D = eo£ + Popt, are the quantities 
that enter into the Maxwell equations. 

An alternative model, more appropriate for semiconductors is that of nearly-resonant coupling between 
two narrow bands of states, assuming that an equilibrium is quickly reached within the two bands separately. 
For example at T=OK, electrons are supposed to decay instantly to the lowest available level of the band 
considered (only one single-spin state electron being allowed in each level, according to the Pauli principle). 
The lower band is called in that case the valence band and the upper band the conduction band, see 
the book by Landau and Lifchitz [41! P- 420]. These authors consider from the outset broadened levels 
and introduce early in their calculation of the electron-optical field coupling population ratios at thermal 
equilibrium. The end result is a relation between the variance of the induced current and the conductance 
(or susceptibility) of the material. 
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In the last expression we have assumed an oscillator strength equal to 1 instead of the 
previously calculated value of 0.96. It is interesting that this expression does not involves 
h. If some external mechanism interrupts the process after a time r much smaller than 
the Rabi period, the time-averaged conductance is given by (11451) with t replaced by r/2. 
Under such conditions the optical potential source "sees" a positive (QM and time)-average 
conductance. A similar discussion applies if the electron is initially in the upper state 2. 
The conductance then has the same absolute value as before, but is negative. 

5.8 Electron energy. 

The electron energy reduces in the present situation to the kinetic energy E = p^/(2m). 
For the non-stationary state previously considered we first evaluate (p^) from (l45l) 

/ + 00 
dpp^ P{p,t) 
-OO 

= cos^(^t) / dp p'^ipi{p)'^ +sin'^{^t) j dp p^i}2{pf 

J OO J OO 

- 2sin(— t) cos(— t) sin(cjot) / dp p'^'ipi{p)ip2{p) 

^ ^ J ~oo 

(146) 

where ijj^Pjt) is given in (11381) . Thus, using previous expressions 

{E) = ^^ = E, cosH^t) + E2 sinH^t). (147) 

If a measurement of the electron energy is performed at some time, the probability that the 
electron be found in state 1 of energy Ei is cos^(^t) and the probability that the electron 
be found in state 2 of energy E2 is sin^(^t). 

To summarize, we have evaluated in the present section the QM-averaged optical cur- 
rent {i{t)) induced by an electron residing between plates submitted to static and optical 
potential sources. The power v{t) (i(t)), when averaged over an optical period, is denoted 
P(t). If the electron is initially in the lower state, the energy E(t) supplied by the optical 
potential up to time t, namely the integral of P{t) from t = to t was found equal to the 
electron QM-average energy at time t. 

If an electron initially in the lower state interacts with the field during a fixed time r, 
with Qrt <^ 1, it has a small probability of being in the upper state at the exit time. If 
this happens, a detection event is recorded. Because the probability is small, the detection 
events form a Poisson point process, even if the electrons are injected regularly, say at 
times t = 0, r, 2t, .... On the other hand, we have shown that a constant (i.e., field - 
independent) positive conductance is obtained. In the next section, we show that such 
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a Poisson distribution may be interpreted by supposing that a current source of spectral 
density huJoG is associated to a conductance G, in the condition previously considered that 
the electron remains most of the time in the absorbing state. This is the Nyquist-like 
formula. This optical noise-current itself is not measurable. But when the complete circuit 
equations are solved, we obtain accurate expressions for measurable photo-currents. 

In subsequent parts of this paper we suppose that G depends on parameters such as 
the number of electrons. The relative variations of these parameters and of the conduc- 
tances are small (e.g., An/n ^ AG/G 0.01), and therefore the spectral density of the 
Nyquist-like noise sources may be supposed to be unaffected. The dependence of G on 
these parameters, however, entails drastic changes in the system behavior. In particular it 
explains the different behaviors of light-emitting diodes and above-threshold lasers. In the 
latter case, intensity fluctuations are very small compared with the mean intensities. 

5.9 C-state 

Light waves are said to be in the C-stat^ if they generate Poissonian photo-electrons 
irrespectively of the carrier phasj^. We show that potential or current sources radiate 
light in the C-state. 

Let us first recall well-known observations. The current emitted by a cathode whose 
emission is temperature-limited consists of independently emitted electrons. Mathemati- 
cally, this electronic emission process is referred to as a Poisson process. Let the emitted 
current be denoted by ^7 = J -|- A J, where J denotes the time-averaged emitted current. 
The spectral density of the fluctuation AJ is given by the formulas 

^AJ = eJ, (148) 
^AD = D, (149) 

where e denotes the absolute value of an electron charge and double-sided spectral densities 
are employed. The letter "D" stands for "detection". The electronic rate is deflned as 
V = J/e = D + AD, with D = J/e and AD = AJ/e. Relation ([148]) says that the average 

C-state beams resemble the so-called "coherent" states of light employed in quantum optics. However, 
C-states are fundamentally states of propagating light while coherent states are primarily states of optical 
resonators. In the context of Quantum Optics, Glauber has shown in 1963 that a classical prescribed current 
(which we call a current source) radiates light in the so-called "coherent state". When coherent states are 
incident on a photo-detector the statistics of the photo-electrons is Poissonian. The results therefore are 
similar. 

^^Concretely, the carrier phase may be changed by inserting on the optical beam (i.e., before detection) 
a second-order all-pass filter, which is a conservative (i.e., lossless, gainless) resonating device that changes 
the carrier phase without changing its amplitude. We suppose that this circuit bandwidth is very small 
compared with the Fourier frequencies of interest. If this is the case, the carrier phase may be changed 
arbitrarily from to 27r simply through a very small detuning. The fluctuations, on the other hand, are 
essentially unaffected by that all-pass filter. The description of second-order all-pass filters may be found 
in Circuit-Theory textbooks. 
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power dissipated in a Ifi resistance following a 1 Hz band-pass filter centered at any low 
frequency (white noise) is given by the shot-noise formula 

2^AJ = 2eJ, (150) 

If a light source, for some reason, is supposed to emit photons independently each photon 
carrying an energy huio, where h denotes the Planck constant and lOq the light frequency, we 
are dealing again with a Poisson process. The spectral density of the light-power fluctuation 



AP is thus 

^AP = fl^oP, (151) 

^AQ = Q, (152) 

where the average photonic rate is Q = P/tkJo and the fluctuation is AQ = AP/hiOo- 



If the optical beam is incident on an ideal photo-detector, the light is converted into an 
electron rate identical to the photon rate. Thus, the electron rate fluctuation spectral 
density S^ad = D = Q. However, the above discussion suggests that the introduction of 
photons is here superfluous. 

In the present theory complex random current sources C{t) = C'{t) +\C"{t) are associ- 
ated with positive conductances G (Remember that G{t) is the complex representation of 
a signal at the carrier optical frequency uj). The real and imaginary parts of C and C" of 
C, respectively, are uncorrelated and have spectral density 

,9'c' = ^c" = f^oG. (153) 

Let us show that this formalism agrees with the previously-quoted shot-noise formulas. We 
have established in the previous section that when a large number of electrons are submit- 
ted to optical potential sources photo-detection events are Poisson distributed. We thus 
consider a potential source applied to a conductance G endowed with its Nyquist-like noise 
source, and show that the photo-electrons are indeed Poisson distributed, irrespectively of 
the potential-source phase. The interest of this representation is that it dispenses us of 
considering real devices involving electrons. We substitute to it a simple circuit schematic. 

Consider indeed a potential source V applied to a conductance G. The current delivered 
by the source consists of two parts. First the current GV flowing through the conductance, 
and secondly the noise current G{t). It follows that the power delivered by the potential 
source (and received by the detector according to the law of average-energy conservation) 
reads 

P{t) =^{V*{GV + G)} = G\V\'^ + V'C' + V"G", (154) 

2 

The first term is the average power P = G\V\ . The second terms are fluctuating terms. 
The spectral density of the fluctuation reads (see Section 13. 6|) 

^AP = V'^^c' + V'^^c" = G\Vfhwo = Phwo. (155) 
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Thus, setting D = P/hwo, AZ) = AP/htOo, we recover the relation S^ad = D. The photo- 
current is Poissonian, irrespectively of the phase of V. A similar calculation can be made 
for current sources. Thus we have shown that potential and current sources radiate light 
in the C-state. To summarize, under the conditions outlined above, detection events are 
Poisson distributed and the conductance is constant. This Poisson behavior is accounted 
for by postulating Nyquist-like currents. 

In the present section we were mostly concerned with electrons that interact with the 
field during some small time r, and reside most of the time either in the absorbing state 
(positive conductance) or in the emitting state (negative conductance). In the subsequent 
section we consider electrons that interact permanently with the field, but may sponta- 
neously transit from one state to the other at some known average rate 27. 
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6 Generalized electron-field interaction 



We consider in the present section an electron present all the time in the interaction region, 
with spontaneous transitions between the two levels considered above. Unlike the spon- 
taneous events usually considered in the Quantum Optics literature, these spontaneous 
transitions do not involve an irreversible loss of energy. They merely convert the electron 
energy into an energy received (or delivered) by the static potential source. 

The waiting-time distribution w{t) is defined as follows: Given that the electron is in the 
lower state at i = (implying that a transition event just occurred), w{t)dt is the probability 
that the next transition occurs between t and t + dt. We first give an approximate form of 
the waiting-time distribution. Then we introduce generalized Rabi oscillations and evaluate 
exactly the waiting-time distribution and the power P{t) supplied by the optical-potential 
source to the electron and to the static-potential source. 

6.1 Waiting time probability for small decay rates 

We considered in the previous section an electron which may reside in either one of two 
states separated in energy by HlOq ~ e.U . If an optical potential source cos{ujot) is 
added, the electron undergoes Rabi oscillations. To wit, if the electron is initially {t = 0) in 
the lower state, the probability that it be in the higher state at time t is sin (%it), where 

denotes the Rabi frequency. The energy delivered by the optical potential source is 
equal to hujo when t = n/Qji. This is the energy that the electron would deliver to the 
static potential source if it were to tunnel through it at that time. 

If an event occurred at i = 0, the probability density of a downward transition is 
approximately the product of the probability sin^(^t) that the electron be in the emitting 
state at time t and some constant 27. The transition events form an inhomogeneous Poisson 
process of rate X{t) = 27sin^(^t) = 7(1 — cos(ili{t)). A change of the time-scale r = T{t) 
with dr(t) = X{t)dt transforms this inhomogeneous Poisson process into a homogeneous 
Poisson process of density 1. For such a process it is known that the probability that an 
the next event occurs between r and t + dr is exp(— r)(iT, see Section [331 Let us set 

fit) = exp(-r(t)) = exp (- dtA(t)) = exp (-7(t - ^^^)) 

(156) 

If an event occurs at i = 0, the probability W(t)dt that the next event occurs in the interval 
[t, t + dt) reads after integration 

W{t)dt = exp{-T{t))dT{t) = -df{t) = f{t)\{t)dt 

= 7(l-cosJ7ijt)exp(-7(t 7^))dt. (157) 
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Figure 6: The figure represents the waiting-time density W{t) according to the approximate 
result in (11571) (plain line) and the waiting-time density w{t) according to the exact result 
in (fToT]) (dashed line) for 7 = 

Prom now on we select a time-scale such that flu = 1. Figure [6] represents W{t) for 
7 = J7/j/7 = 1/7. The expression for W{t) given above almost coincides with the one, w(t), 
known to apply to resonant fluorescence |42i see Figs. 8.4 and 8.5][43] when <C 1. The 
exact expression of the waiting-time probability density w{t) is derived in Section [Ql Note 
that the probability that no transition event occurs before time t is the integral of W{t) 
from t to 00, that is, f{t). If we knew the time at which the last event occurred before time 
i = 0, the electron state at t = would be a pure state. If, however, we only know the 
probability density that this event occurred at time t < 0, the electron state at t = must 
be described by a mixed-state density matrix. This is the situation described in the next 
sections. 

In general, one should allow for the fact that an electron, having a probability sin^(^t) 
of being in the upper state at time t and therefore a probability cos^(^f) of being in the 
lower state, has not only a probability density 272sin^(^t) of performing a transition 
to the lower state as we discussed above, but also a probability density 27icos^(^t) of 
performing a transition to the upper state. Previously, we assumed that 71 = and 72 = 7, 
but this needs not be always the case. 

6.2 Density matrix 

When there is some uncertainty concerning the initial form of the wave-function, it is con- 
venient to employ the 2x2 mixed-state density matrix whose mathematical properties were 
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given in Section [321 We first consider the Rabi equations as given earlier, and subsequently 
introduce the generalized Rabi equations. 

The wave-functions in space and momentum space are of the form given in (|129p and 

(HMD 

'tp{x,t) = Ci(t) exp(-iu;it)V'i(x) + C2it) exp{-\U2t)'4'2{x) 

i;{p,t) = C7i(t)exp(-iLJit)Vi(p) +C2(t)exp(-ia;2t)V'2(p). (158) 

where the complex coefficients Ci{t),C2{t) obey the differential equations 

^=^^C2it) ^=^^C,it) C,(t)Cm+C2{t)Cm = h (159) 

and <^ uJo is the Rabi frequency given by hQji = ex\2'\piV j d. We assume that 
Ci (0)6*2(0) is imaginary by a proper selection of the phase of the wave-function. It then 
follows from the above differential equations that C\{t)C2(t) is imaginary at all times. 
Proceeding as in the previous section, we evaluate 

/+00 
dv p 'il;{p,t)^*{p,t) = Ci{t)C^{t)eM^^ot)pi2 + cc, (160) 
-00 

defining 



Pij= dp p i)i{p)i}*{p) i,j = 1,2. (161) 

J —00 

and assuming as before that pu = P22 = 0. For the situation presently considered pi2 = 
—if}iOoXi2 = P2i- Because Ci{t)C2{t) is imaginary, Ci(t)C|(t)pi2 is real, and we are left 
with a term proportional to cos{uJot), in phase with the driving field. We obtain the average 
induced current {i(t)) = ^^^^ as in (I143|) . The instantaneous power delivered by the optical 
potential source is v{t) {i{t)), where v{t) = V2V cos{u!ot) and hQ^ = ^/llLe£l2_ Averaging 
over an optical period, that is, replacing cos^(a;o^) by 1/2, we may write the power in the 
form 

P{t) = hiOonRQ{Ci{t)CUt)} = fiuJoVLRp'i2 (162) 
Let us now define a 2 x 2 density matrix 



pii(i) pi2{t) \^( c^mtit) c^{t)cm 

P2l{t) P22{t) )-\ C2{t)C{{t) C2(t)Q(t) 



(163) 



The p(f)-matrix is Hermitian, that is, equal to its complex conjugate transpose and trace{p'^{t)} 
trace{p{t)} = 1. 
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We readily obtain the differential equations obeyed by the elements of the p(t)-matrix 
from those in (|159p as 

dpnjt) XIr, 

dp22{t) XIr. ... 

^=^^(m(t)-p.(^)) 

^=^^(p.,(t)-pn(^)). (164) 

The power delivered by the optical potential given in (|162p may be written in terms of the 
density matrix as 

P{t) = trace{p{t)P} (165) 

where 

P = ( I ) (166) 



In view of the mathematical considerations in (|3.2p . it follows that the above expression 
for the power in (|165l) holds also for mixed states, in which case p is the weighted sum of 
pure-state density matrices and is written as 

Pn{t) Pi2{t) \ .^.„. 
P2i(t) P22{t) ) ^''^^ 

without any reference to the Ci(t),C2(t) coefficients anymore. In other words, the density- 
matrix concept has been introduced to generate an expression of the power delivered by 
the optical potential in the case where the wave-function is defined only statistically. We 
are now ready to generalize the Rabi equations. 



6.3 Generalized Rabi equations 

We have previously evaluated the conductance "seen" by an optical potential source under 
the assumption that the electron is submitted to the optical field during some known fixed 
time r. We now go back to the continuous configuration, and introduce instead a known 
average decay rate 27. This leads to a generalized form of the Rabi equations. These 
equations involve a parameter 2a. Prom the fact that the trace of the square of the density 
matrix may not exceed 1 we conclude that 2a may not be less than 1. We justify the 
value 2a = 1 by considering that in the small-7 limit the expression of the waiting-time 
distribution must agree with the expression obtained straightforwardly in a previous section. 
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More generally we assert that the rate of decay of the coherence is as small as is allowed, 
additional possible sources of decoherence (collisions) being not present in our model. 

We again consider an electron with two states, submitted to a resonant field. In the 
absence of spontaneous transitions between states 1 and 2, the 2 x 2-matrix p{t) that 
describes the electron state obeys the differential equation given before in (|164p . If, however, 
an electron in the emitting state may decay spontaneously to the absorbing state at an 
average rate 27, additional terms must be added on the right-hand-side of these differential 
equations. For generality, we consider both a spontaneous-transition rate 72 from 2 to 1, 
and a spontaneous-transition rate 71 from 1 to 2. The previously mentioned equations for 
dp{t)/dt become 

[puit] - P2i{t)) + 272/922(i) - 27i/Oii(i) 
{p2i{t) - puit)) + 27ipii(t) - 272/922(0 
ipu{t) - p22{t)) - 2a(7i + 72)pi2(i) 

{p22{t) - pii{t)) - 2a(7i + 72)p2i(i) (168) 



dpu{t) 


_ i^R 


dt 


2 


dp22{t) 


_ i^R 


dt 


2 


dpi2{t) 


inn 


dt 


2 


dp2i{t) 


_ i^R 


dt 


2 



These equations are unchanged under the interchange of the subscripts 1 and 2. If p is 
Hermitian with trace 1 initially (that is, pn, P22 are real, pu + P22=^, and P21 — P12), 
these conditions are preserved in the course of time. The trace condition expresses the 
fact that the electron resides in either one of the two states considered. The Hermitian 
condition expresses the fact that observables are real. In the low-field limit [Qr = 0), the 
above equations describe accurately the spontaneous decay from level 2 to level 1, and the 
spontaneous promotion from level 1 to level 2. 

The second equation in (|168p (or equivalently the first) may be written as 

rr rr 
{E{t)) = hWo^R / dtp'[^{t) = nWoP22{T) +hWo (272/922 (t) " 27lPll(t)), (169) 

Jo Jo 

where we have assumed that the electron is initially in the absorbing state P22{0) = 0, and 
Puit) denotes the imaginary part of puit). {E{t)) represents the energy supplied by the 
optical potential source from t = to t = t, that is the time integral from to r of the 
power P{t) given in (|165p . This energy split into two parts. First the energy delivered to 
the electron, i.e., the increment of fiuJoP22{T)- Second, the energy transferred to the static 
potential source through the decay rate 2'j2P22{t), minus a similar term for the upward 
spontaneous transition. Thus the first two equations in p68|l are almost obvious, as they 
express the law of average- energy conservation. From now on, for brevity, we set = 1. 

As far of the third equation is concerned (or equivalently the fourth), let us note that if 
initially the trace of equals 1 (pure state), the time derivative of the trace of p^(t) may 
not be positive since trace{p{t)} < 1. For 71 = or 72 = 0, a short calculation shows that 
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this implies that 2a > 1, see Section 13.21 in which we set 7 = 7i + 72, ^7 = 7i — 72- Since 
for a single electron there are no decoherence due to elastic collisions, we may assume that 
the decoherence rate is as small as possible, implying that 2a(7i + 72) = — \/l2) ■ 

After a sufficiently long time a steady-state is reached, obtained by setting the time 
derivatives equal to 0. We have p'l^ioo) = and 



Pl2(00) 



72-71 7 7 



, . 1+4071(71+72) 1 1 



For 71 = 72 we have, as expected because of symmetry, ^"2(00) = 0, /322(oo) = />ii(oo) = |. 
In the second expressions above we have set 71 = 0, 72 = 7 and in the third we have further 
set 2a = 1. 

6.4 Solution for arbitrary times 

Let us now solve the equations for arbitrary times. Using the fact that p2i = P121 Pii = 
1 — P22) P22 real, we only need the second and third equations. Taking the imaginary part 
of the third equation we have 

= 1- P22{t) - 2a(7i + 72)p'i'2(i), (171) 
The second equation may be written as 

P'i2it) = - 271 + 2(71 + 72)p22(t). (172) 

Substituting p'{2{t) from (11721) into (|17ip we obtain for 71 = 0,72 = 7 the second-order 
differential equation for P22{t) 

d^p22{t) , , ^ dp22{t) I 2\ (.\ ^ n7Q^ 

^^2 +2(1 + a) + (I + 407 jp22(t) = -. (173) 

Considering the function y{t) = p22{t) — 1/(2 + 8072), the constant term cancels out. Next, 
we try as usual solutions of the form y{t) = exp(At) and find that A must be a solution of 
the second-degree equation 

+ 2(1 + a)7A + 1 + 407^ = 0, (174) 

whose solution is 



A± = -(1 + 0)7 K= V(l- 0)272-1. (175) 
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With the initial condition P22(0) = the solution is of the form 

1 

2 +"807^ 

Since the condition pui^) = implies that the first derivative of P22{t) with respect to time 
vanishes at f = 0, the expression of P22{t) is 



'°22(i) = iTV^T^i^ + Aexp(A+t) - (1 + ^) exp(A_t)}, (176) 



P22{t) = ^^{1 + ^exp(A^t) - ^exp(A_0} (177) 



Note that this expression is always real. We interpret G{t) = 27^22(0 as the probability 
density of having an event at time t, not to be confused with the conductance G{V), where V 
denotes the applied potential. Since an event occurred at t = 0, the normalized correlation 
reads g{T) = G(r)/G(oo), where G(oo) = 27^22(00) = x+4a7^ probability density of 

the process. The relative noise is obtained straightforwardly from the Fourier transform of 
(7(r) according to ([63]l . 



6.5 Correlation when 2a=l 

When 2a = 1, the expression of G{t) = 27^22 (i) in (11771) reads 



In particular 

7»1 

27 

G(t) Ri7(l-exp(-^)cos(t)) 7<C1 
G(t)=2^-(^ + ^;)-P(-^^) 7 = 2 

6.6 Average conductance and fluctuations when 2a = 1 

Let us rewrite (11781) with the Rabi frequency O/j restored. We set 70 = 'y/^R- 

1 + 27q zk zk 
0M = ^ (180) 
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In the large 7 limit 



Go = G(oo)^^oo = (181) 

Since each event absorbs an energy hujo and Go (not to confused with a conductance) 
denotes the average number of events per unit time, the power delivered by the optical 
potential V to the static potential U is 

P^GV' = Gofu.o=^^. (182) 

Because, according to (I132|) . fi^ = 2 ( ^^^^ ) ^ is proportional to V'^, the conductance G does 
not depend on the potential V, that is, the system is linear in that limit. The conductance 
has the same form as the one given in (|145ll for an electron returning to the absorbing 
level after a fixed time r, if we set r = 2/7. For larger values of V the conductance 
decreases, and tends to 0. This non-linearity may be interpreted as usual from the fact 
that the lower population decreases and the upper population increases as a consequence 
of stimulated absorption. Indeed, the conductance is for arbitrary 7-values equal to the 
previous low-field expression multiplied by pu — P22 = 1 — 2p22 = 27o/(l + 27^). 

Let us now consider the normalized events correlation ^(t) = G{t)/G{oo). For large 
7-values we obtain 

5f(r) = (1 - exp(-7t))^ ^ 1 7^00 

POO 

AA(0) = 2 dT{g{T) - 1) ^ 7 ^ 00. (183) 
Jo 

That is, the relative noise M vanishes at zero Fourier frequencies, so that the point-process 
may be considered Poissonian. For arbitrary 7-values, the spectral density of the point pro- 
cess is viewed as the sum of two terms. One originating from the lower-state population pu 
(positive conductance), and one originating from the upper-state population P22 (negative 
conductance). If they were independent, the corresponding spectral densities should add 
up, that is, be proportional to pu + P22 = 1- We have in fact from (|178p 



- 1 = -^{X- exv{^+^Rt) - A+exp(A_QRt)) 
A± = -^±K ^=M-h (184) 



9{r) 



2 

with g{oo) = 1, g((0) = 0. The relative noise at zero Fourier frequency is 

Af{0) = 2 rdrigir) - 1) = -^(^ - ^) = (185) 







kQr^X- \+^ 1 + 272' 
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It follows that the fluctuations corresponding to the negative and positive conductances are 
correlated. This is to be expected because we are dealing with a single electron. In later 
applications, we consider instead two large collections of electrons, one in the emitting state 
and one in the absorbing state. In that case the fluctuations are independent. 



6.7 Waiting time probability 

The process we are presently considering is an ordinary renewal point process. For such a 
process it is straightforward to go from G{t) to the waiting time probability density 'w{t). 
The concept is that the probability density of an event at t is the sum of the probabilities 
that this occur through 1 jump, 2 jumps, ...and that k-jumps probability densities involve 
k-fold auto-convolutions of w{t). The formula relates the Laplace transform G{p) of G{t) 
and the Laplace transform w{p) of w{t), see Section \TM 

It follows straightforwardly from the fact that G{t) involves a constant term and that, 
consequently, G{p) behaves as 1/p for small p, that 

^ = ^(-) (-) ^[dttwit) = - (^),=o (187) 

implying conservation of the average power. 

To obtain w{p) we must evaluate the Laplace transform of G{t). Since the Laplace 
transform of exp(— at) is l/{p+a) and Laplace transforms are linear operations, the Laplace 
transform of G{t) = 2jp22{t), where P22{t) is given in (I177p . is, with X± = — (l-|-a)7±K, k = 
y^(l — a)^^^ _ \ as before, 

.1 A_/(2k) \+/{2k) 



G(p) = / dtexp(-pt)27p22(t) = ^J. 2 - + 
1 A+)(p- A_) 



p — A_ 



G{p) 7 
From (fTMD and (fTHHIl . we obtain 

w{p) 



7 



p{p-K){p-^-) +7 

7 



p3 + 2(1 + a)7p2 + (4072 + l)p + 7 



(189) 



We will first consider the case where 2a = 1 and subsequently give w{t) for arbitrary 
a- values. 
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6.8 Waiting time probability when 2a = 1 

Let us consider first the case 2a = 1. We have 

w{p) = :^( + 2^^) a = - 1 (190) 

2a^ p + 7 + a p + 'J — a p + j 

Taking the inverse Laplace transform of w{p) in (|19QI) we get 

w{t) = -^{exp(-(7 - a)t) + exp(-(7 + a)t) - 2 exp(-7t)}. (191) 
2a^ 

which is real non-negative for all 7-values and integrates to unity. At small times we have 
G{t) « w{t) « 7*^/2. Special forms are 

w{t) = ^exp{-^) t^O, 7>1 
27 27 

= i{ exp(-(2- \/3)t) +exp(-(2 + \/3)t) -2exp(-2t) } 7 = 2 

= 7(l+72)(l-cos(t))exp(-7t) 7 < 1 (192) 

In the small 7 limit, the expression in (|192p coincides with the expression obtained on 
intuitive grounds in (|157|) . Note that the three expressions in (|192p integrate to unity. 

6.9 Arbitrary a-parameter 

We now have to find the roots of p^ + 2(1 + 0)7^^ + (407^ + l)p + 7 = in (I189p . This is 
done by setting in and dZ]) ao = 7, oi = 407^ + 1, 02 = 2(1 + 0)7. We then obtain w{t; a) 
from the Heaviside theorem in (l48]l . and evaluate it numerically. We define the "distance" 
from W{t) in (|157p and w{t;a) as defined above, in terms of the dimensionless parameter 



00 



A(a) = 10^(T)/ {w{t;a)-W{t)y, (193) 







where (r) denotes the average waiting time. We have obtained for 7 = 0.001 the following 
values 



2a = 0.996, 0.998, 1, 1.04 
A(a) = 1.6, 0.066, 0.1, 3.2 ^ ' 

It follows from these numbers that the expression of w{t) in the 7 — > limit for 2o 7^ 1 
does not coincide with an almost obvious result. This is, it seems to us, a strong argument 
for setting 2a = 1 in the starting Rabi equations. 
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6.10 7i = 72 

The generalized Rabi equations in (I168P involve both a downward transition rate 272/322 
and an upward transition rate 271P11. Plausible expressions for the two parameters are 

/ ell — fluJo\ 

71,2 = ■^^exp ( lb ^— I tz7 <C Wo- (195) 

We will not attempt to justify these expressions here. If ell = strictly, 7i = 72 = 
ro, 6 = 0. The Rabi equations in vector form ([2511 become, setting a = 0, x{t) = 

dy_ _ _ 

dt ~ ^' 
dz 

-^ = y- 2ro^- (196) 

In the large-time limit we have y = z = 0, that is, there is no population inversion and 
the net conductance vanishes. Fluctuations however remain. In a small time interval 0, dt 
we may have zero event, a photo-absorption event (+) or a photo-emission event (-). As 
said above, the average numbers of (+) and of (-) are equal so that the average current 
vanishes. Solving the above equations and supposing that at time t = a photo-absorption 
event occurred {pu{0) = 1 =^ z{0) = — l,y(0) = 0), we obtain the probability density for 
the sum of (+) and (-) events as G(t) = ^z{t), which can easily be obtained explicitly. 

liw<^ = 1, the electron performs many Rabi oscillations before a transition of any 
kind occurs. It follows that we can average the transition probabilities over a Rabi period, 
so that the populations are essentially equal to 1/2. In that limit, the electron population 
follows a telegraphic random motion between state 1 and state 2, with equal upward and 
downward transition probability densities. The spectral density of the photo-current at 
zero Fourier frequency is then the sum of that of the upward and downward processes. If 
the condition <C 1 is not fulfilled a more involved analysis would be required. 
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7 Alternative approaches to Nyquist-like noise sources 



The present theory is based on two fundamental concepts. One is the introduction of 
Nyquist-like current sources, the other is the law of conservation of average energy. In the 
present section we consider alternative approaches to the link that exists between fluctuation 
and dissipation, usually referred to as the "fluctuation-dissipation theorem". Fluctuating 
rate sources are uncorrelated and their spectral densities are equal to the average rates. 

Brown was the first to observe the random motion of small particles in viscous fluids 
(Brownian motion) and Johnson later on measured the electrical noise associated with 
conductances (Johnson noise). The first interpretation of Brownian motion was offered by 
Einstein and the first interpretation of electrical noise by Nyquist. The two phenomena 
are closely related. The Einstein and Nyquist interpretations rest on Classical Statistical 
Mechanics {k-QTm » fiuj). For our purposes we need consider instead the situation where 
teTm ^ The latter "quantum" situation was treated by Callen and Welton in 1951, 
see for example |44|. We first recall the Quantum Optics viewpoint. We consider next 
simpler, but admittedly partly heuristic, explanations. These arguments tend to prove 
that the spectral density of the Nyquist-like noise sources associated with two-level atoms 
exhibiting a peak conductance G at frequency a; = is equal to hwo\G\ when either 
/Oil ~ 1 (positive conductance) or p22 ~ 1 (negative conductance). 

In the first argument we show that there is a (perhaps unique) way of generalizing to the 
quantum domain the classical expression {E) = k^Tm for the average energy of an oscillator 
resonating at frequency uJq that avoids divergences. The expression obtained differs from 
the 1901 Planck formula by an additional energy hjJo/2 called the vacuum energy. This 
expression does not lead to "ultraviolet" divergences if only measurable quantities are being 
considered. 

In the second approach we consider an oscillator in a state of thermal equilibrium 
with a small negative conductance and a slightly larger positive conductance (in absolute 
value). The ratio of these two conductances is equal to the ratio of lower and higher-state 
populations, given by Classical Statistical Mechanics. 

The third approach is quite different from the previous one since we consider an isolated 
oscillator containing two-level atoms in a state of highly-non-thermal equilibrium. Statisti- 
cal mechanics tells us that states of equal energy are equally likely to occur, as long as no 
information concerning these states is available. The noise terms (Langevin forces) must be 
such that the variance of the photon number derived from the above Statistical Mechanical 
law be obtained in the long-time limit. 

7.1 Quantum-optics approach 

In Quantum Optics treatments, loss-less optical oscillators oscillating at frequency uJq are 
viewed as being akin to loss-less mechanical oscillators with quantized energy levels Em = 
hu>o{m + m = 0,1.... If we apply to such oscillators the Boltzmann result that the 
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probability that a level of energy E be occupied is proportional to ex.p{—(3E), (3 = 1/k-QTm 
when the system is in contact with a bath at absolute temperature (subscript "m" 
unrelated to the integer m), the Planck law of black-body radiation results through a 
summation over m. We obtain 

i _ Em=o-^'» exp(-^^n^) _ nwoeyi-p{l3nuJo) + 1 

The fluctuation-dissipation theorem (FDT) applies to systems having a linear causal 
response, in thermal equilibrium. One first evaluates the symmetrized correlation Cj^q(T) = 
(A?(r)Ag(0) + Ag(0)Ai(r)) /2, where the charge q and the current i = dq/dt are conjugate 
operators (as x andp/m = dx/dt in Quantum Mechanics). The end result of the calculation 
is that the spectral density of the current, that is the Fourier transform of the correlation 
Cj,j(T) reads [B] 

where TZ{Y{lj)} = G{uj) denotes the system conductance, a function of the frequency uj. 
In the following, operators are not employed, that is, i{t) is viewed as a classical function 
of time. 

7.2 Heuristic approach 

A simple heuristic derivation of the black-body formula is suggested here. In the classical 
regime, single-mode oscillators at frequency have a probability exp(— /?£') of having an 
energy E, according to Classical Statistical Mechanics. Thus the average energy reads 

, ^ dEEexp{-^E) 1 ^ ^ 

independent of u>o- The action / of an oscillator is the ratio of its average energy and 
frequency, and accordingly /(x) = 1/x, setting x = (3E = w/kBTm- Thus /(x) obeys the 
Riccati equation 

^ + /(x)2 = 0. (200) 

However, as was noted at the end of the 19th century, the expression (E) = k-QTm leads to 
infinite radiated heat since the number of electromagnetic modes is infinite in a cavity with 
perfectly reflecting walls. It apparently did not occur to scientists at the time that this 
difficulty is resolved simply by adding a constant {fi/2) on the right-hand side of (|200p . 
that is, supposing 
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where h is a universal constant with the dimension of action. The solution of this modified 
equation reads 

where Xo denotes an arbitrary constant. For Xo=0 we obtain that 

which coincides with the Planck formula except for the vacuum energy mentionned above. 
The arbitrary constant Xq on the right-hand-side of (|202p must vanish to obtain agreement 
with the classical result. It was noted by Einstein and Stern in 1913 that /(x) — 1/x ^ 
if X ^ 0, that is, the expansion of /(x) is of the form 

fix) = - +ax + (204) 

X 

where a is a constant, without an x-independent term. One may conjecture that there are 
no other differential equation but (|201ll which possesses only one solution satisfying the 
above classical limit. To conclude, the classical expression of the average energy generalizes 
to 

where (3 = 1/k-BTm- If we next consider a cavity with perfectly-conducting walls, solutions 
of the Maxwell equation exist only for a series of real resonating frequencies wi, UJ2-, W3, .... 
Each of these modes of resonance is ascribed an average energy given by the above expression 
with bJo replaced by wi, a;2, ws, .... For a d-dimensional cavity the mode density /o(w), where 
p{oj)duj denotes the number of modes whose frequencies are between uj and oj + duj, grows in 
proportion of w^"^. It follows that the total energy is apparently infinite. This is perhaps 
why, in his original work, Planck subtracted the vacuum energy hjJo/2 from the expression 
given in (|205p . The total energy in a cavity of volume V in thermal equilibrium at absolute 
temperature is then found to be finite and proportional to the fourth power of T^. If 
the cavity is pierced with a small hole that does not perturb much the state of thermal 
equilibrium, the measured output-power spectral density is supposed to be proportional 
to the internal energy spectrum, that is, to the product of {E{oj)) — hijj/2 and the mode 
density p{lo). This is the famous black-body spectrum measured around 1899, and obtained 
by Planck in 1900 on the basis of an "educated guess" of the entropy function. This formula 
agrees very well with measurements. In our model the power is collected by a small external 
conductance g at Tm =0K, having its own noise source. A finite result is then obtained 
from the expression of {E{oj)) as given in (I205|) . 
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7.3 Johnson-Nyquist noise 



When a frequency-independent conductance G is in equilibrium with a bath at absolute 
temperature Tm, there is associated with it a random current source j{t) whose (double- 
sided) spectral density equals IkBT^G, as discovered by Johnson and Nyquist in 1927, 
1928. This expression holds in the so-called "classical regime", that is, at frequencies much 
smaller than k-QTm/h. In the present section currents are denoted by j instead of i because 
we are mostly concerned with low frequencies. 

To justify the Nyquist expression in the classical regime, let us recall that, according to 
Statistical Mechanics, the average energy must equal /cBTm/2 per degree of freedom. As a 
consequence, the energy stored in a capacitance C with a conductance G in parallel must 
equal k^Tm/'i- The modulus square \V/I\ of the impedance of the circuit considered at 
frequency a; is Since the energy stored in a capacitance C submitted to an 

(root-mean-square) voltage V is C|y|^ /2, as we have seen in section \AA\ we must have 



from which it follows that = 2k-QT„iG. What we have given above is essentially the 
Nyquist argument. 

The Nyquist result may be generalized to any (time-invariant, linear) circuit consisting 
of any number of conductances all of them being at the same temperature Tm, capacitances 
and inductances. If the conductance between any two terminals of the circuit is G{oj), the 
circuit may be considered noiseless provided a current source j{t) be applied to the two 
terminals with spectral density ^j(w) = 2k-BTmG{uj). 

For later use note that when a white (i.e., with a frequency-independent spectrum) 
Nyquist current i{t) is applied to an ideal narrow-band filter whose response is centered at 
uj = ibwo, the filter output may be written in the form 



where the real random functions of time C'{t) and C"{t) vary slowly, are uncorrelated, and 
have (double-sided) spectral density = = ^.k^TrnG. 

7.4 Thermal equilibrium approach 

In the present section we generalize the previous Johnson-Nyquist expression to arbitrary 
frequencies. Consider an inductance-capacitance (L — C) oscillator resonating at frequency 
u!o with a small positive conductance Ga and a negative conductance —Ge (subscript "e" 
for "emitting") in parallel, with Ga (subscript "a" for "absorbing") exceeding Ge so that 
there is a small net loss. We set Ga — Ge = G > 0. The conductance Ga is supposed 
to represent the absorption by two-level atoms in the lower state (energy Ea), while Ge 
represents the emission from atoms in the upper state (energy Ee). A near-resonance 




(206) 



i{t) = V2C'{t) cos{uJot) + V2G"{t) sin(cJot) = ^{V2G{t) exp(-ia;ot)} 



(207) 
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Figure 7: a) represents a potential source V applied to a conductance G. C represents 
the Nyquist-like noise current source associated with G. We are mostly concerned with the 
power flowing from left to right through the dashed line, b) represents a conductance driven 
by a current source I, c) is a laser model with a negative conductance — Ge, a resonating 
circuit, and a positive conductance G^, representing the detector of radiation, connected in 
parallel, d) represents a laser model that exhibits a Petermann-like linewidth-enhancement 
factor. 



condition hwo ^ — Ea holds. According to the Schrodinger equation the conductances 
are proportional to the corresponding numbers of atoms Ua and rig, respectively. 

Because the conductances considered are constant the circuit is linear. The potential 
V across the circuit is therefore equal to C/Y{uj), where C denotes the driving current, 
assumed to be independent of frequency, and the resonating circuit admittance Y{uj) = 
G + \B{(jj) where B{ijj) represents the sum of the L and C susceptances. Referring to (1107(1 
the oscillator energy is 

~ J^oo ^G^ + AG^{Lo-iVof ~ 4G^i_^ l + x2 " 4G- ^ ' 

Let this oscillator be in thermal contact with a bath at absolute temperature T^. We 
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replace the deterministic current C by a complex random function of tima^. Because the 
processes are stationary we expect that the statistical properties of the random source 
C{t) = C'{t) + \C"{t) are unaffected by an arbitrary phase change, that is, we require 
that C(t) exp(ii;^) has the same statistical property as C{t) for any phase (f). This entails 
that C'{t) and C"{t) are uncorrelated and have the same statistical density. We thus set 
= = 5^ . We have seen that for n independent atoms in some state both G and 
5^ are proportional to n. We therefore expect that 5^ = aG, where a is a constant to be 
determined. 

Because of the symmetry between stimulated absorption and stimulated emission im- 
plied by the Schrodinger equation the spectral densities have the same form for positive 
conductances Ga and negative conductances — Gg, namely 5^a = ciGa and = aGe, with 
the same constant of proportionality a. Note that C and G" contribute equally and that 
double-sided spectral densities for C", G" in the Fourier O-domain are employed. If the 
conductances Ga and —Ge are connected in parallel the total conductance is G = Gq — Ge, 
as said above, and the total spectral density is 2a(Ga-|-Ge). The average resonator energy 
follows from (12081) 

2a(G, + Ge) « GJGe + l 
= 4G = 2 GJGr^T- (209) 

Classical Statistical Mechanics tells us that at thermal equilibrium ^ = exp(-^^^). 
It follows that 

^ = - = -P (^^) = exp (-^) (210) 

Classical Statistical Mechanics also tells us that in the classical limit the average energy 
equals kBTm/2 per degree of freedom, and thus (E) = k-QTm for the oscillator considered 
when k-QTjn » hwo. According to (|209p this is the case if and only if 

a = tvjJo (211) 

At low-temperatures, k-^Tm ^ the resonator energy is therefore hjJo/2. Since a 
resonator may exchange energy resonantly with atoms only in units of fvjOo, we conclude 
that the energy of a single-mode loss-less resonator may be written as 

Et = {m + l/2)hwo, (212) 

where m is an integer that one may call "number of photons in the resonator". However, we 
view linear loss-less resonators of any kind (mechanical, optical, or other) as abstractions. 
Such oscillators acquire physical meaning only when they are coupled to sources of energy. 
What need to be quantized are those energy-coupling mechanisms. 

^'^Note that in the present linear regime the regulation mechanism at work in above-threshold lasers does 
not occur, and the fluctuations of V are comparable to average values. Supposing that the current source is 
gaussian distributed, this is also the case for the optical potential V and optical current /. Power should in 
general be evaluated as the real part of V* {I + C{t)), but in the linear regime presently considered the term 
C{t), much smaller than the fluctuations of 1, may be neglected. The power 1R{V* 1} is Rayleigh-distributed. 
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7.5 Isolated cavity approach 



A single-mode optical cavity resonating at angular frequency ujo may be modeled as an 
inductance-capacitance {L,C) circuit with LCool = 1. The active atoms, located between 
the capacitor plates, interact with a spatially uniform optical field through their electric 
dipole moment. The 2-level atoms (with the lower level labeled "o" and the upper level 
labeled "e") are resonant with the field. This means that the atomic levels a and e are 
separated in energy by hiOo- The concept of temperature nowhere enters in the present 
section. 



(Highly non-thermal) equilibrium Consider N identical two-level atoms. For each 
atom, the zero of energy is taken at the lower level and the unit of energy at the upper level 
(typically, 1 eV). The atoms are supposed to be at any time in either the upper or lower 
state. The number of atoms that are in the upper state is denoted by n, and the number of 
atoms in the lower level is therefore N — n. With the convention hujo = 1; tlie atomic energy 
is equal to n. Its maximum value N occurs when all the atoms are in the upper state. There 
is population inversion when the atomic energy n > N/2. The atoms are supposed to reach 
a state of equilibrium before other parameters have changed significantly. The strength of 
the atom-atom coupling, however, needs not be specified. 

The statistical weight W{n) of the atomic collection is the number of distinguishable 
configurations corresponding to some total energy n. For two atoms {N = 2), for example, 
W{0) = W{2) = 1 because there is only one possible configuration when both atoms are 
in the lower state (n = 0), or when both are in the upper state (ra = 2). But W{1) = 2 
because the energy n = 1 obtains with either one of the two (distinguishable) atoms in the 
upper state. For N identical atoms, the statistical weight (number of ways of picking up n 
atoms out of N) is 

iV' 

W{n) = — — -. (213) 

^ ^ n\{N-n)\ ^ ' 

Note that T^(0) = W(N) = 1 and that W{n) reaches its maximum value at n = N/2 
(supposing N even), with W{N/2) approximately equal to 2^ ^/2/TrN . Note further that 

N 

Z=^VF(n)=2^. (214) 

11=0 

Consider next an isolated single-mode optical cavity containing N resonant two-level 
atoms, and suppose that initially all of these atoms are in the emitting (upper) state. 
One observes that the number of atoms in the emitting state may vary from to N, and 
thus one presumes that the missing energy is stored in the optical resonator. The atoms 
indeed perform jumps from one state to another in response to the optical field so that 
the number of atoms in the upper state is some function n{t) of time. If m{t) denote the 
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number of light quanta at time t, the sum n{t) + m{t) is a conserved quantity (essentially 
the total atom+field energy). Thus, m jumps to m — 1 when an atom in the lower state 
gets promoted to the upper state, and to m + 1 in the opposite situation. If atoms in 
their upper state are introduced at i = into the empty cavity (m = 0), part of the atomic 
energy gets converted into field energy as a result of the atom-field coupling and eventually 
an equilibrium situation is reached. 

The basic principle of Statistical Mechanics asserts that in isolated systems all states 
of equal energy are equally likely. Accordingly, the probability pr{m) that some m value 
occurs at equilibrium is proportional to W{N — m), where W{n) is the statistical weight 
of the atomic system. As an example, consider two (distinguishable) atoms (A^=2). A 
microstate of the isolated (matter+field) system is specified by telling whether the first 
and second atoms are in their upper (1) or lower (0) states and the value of m. Since the 
total energy is N = 2, the complete collection of microstates (first atom state, second atom 
state, field energy), is: (1,1,0), (1,0,1), (0,1,1) and (0,0,2). Since these four microstates are 
equally likely, the probability that m = is proportional to 1, the probability that m = 1 is 
proportional to 2, and the probability that m = 2 is proportional to 1. This is in agreement 
with the fact stated earlier that pr{m) is proportional to W{N — m). After normalization, 
we obtain for example that pr(0)=l/4. 

The normalized probability reads in general 

, , W(N-m) N\ , , 

= ^T-^ = 2-m^[N-m)^ ^'''^ 

The moments of m are defined as usual as 

N 

{m^) = m^'pr{m) (216) 

where brakets denote averagings. It is easily shown that (m) = N/2 and var(m) = 
(m^) — (m)^ = N/A. Thus the number m of light quanta in the cavity fluctuates, but 
the statistics of m is sub-Poissonian, with a variance less than the mean. For example, one 
may readily deduce from that principle that, if an isolated single-mode cavity initially con- 
tains no photons but N two-level resonating atoms in the upper state, the system evolves to 
an equilibrium state in which the variance of the number of photons in the cavity is half the 
average number of photons, that is, the photon statistics is sub-poissonian. More generally, 
for atoms with B evenly-spaced levels, the variance of the photon number is (B + l)/6 times 
the average photon number, a result that coincides with the previous one if we set B = 2. 

The expression of pr{m) just obtained has physical and practical implications. Sup- 
pose indeed that the equilibrium cavity field is allowed to escape into free space, thereby 
generating an optical pulse containing m quanta. It may happen, however, that no pulse 
is emitted when one is expected, causing a counting error. Prom the expression in (12151) 
and the fact that (m) = N/2, the probability that no quanta be emitted is seen to be 



107 



pr{0) = 4:~^"^K For example, if the average number of light quanta {m) is equal to 20, 
the communication system suffers from one counting error (no pulse received when one is 
expected) on the average over approximately 10^^ pulses. Light pulses of equal energy with 
Poissonian statistics are inferior to the light presently considered in that one counting error 
is recorded on the average over exp((m)) = exp(20) « 0.5 10^ pulses. 

Time evolution of the number of light quanta in isolated cavities Let us now 

evaluate the probability pr{m, t) that the number of light quanta be m at time t. Note that 
here m and t represent two independent variables. A particular realization of the process 
was denoted earlier m{t). It is hoped that this simplified notation will not cause confusion. 

Let Re{m)dt denote the probability that, given that the number of light quanta is m at 
time t, this number jumps to m + 1 during the infinitesimal time interval [t, t + dt], and let 
Ra{m)dt denote the probability that m jumps to m — 1 during that same time interval (the 
letters "e" and "a" stand respectively for "emission" and "absorption"). The probability 
pr{m,t) obeys the relation 

pr(m, t + dt) = pr{m + 1, t)Ra{m + l)dt + pr{m — 1, t)R(,{m — l)dt 

+pr{m,t)[l - Ra{m)dt - Re{m)dt]. (217) 

Indeed, the probability of having m quanta at time t + dt is the sum of the probabilities 
that this occurs via states m + l,m — lormat time t. All other possible states are 
two or more jumps away from m and thus contribute negligibly in the small dt limit. 
After a sufficiently long time, one expects pr{m,t) to be independent of time, that is 
pr{m,t + dt) = pr{m,t) = pr{m). It is easy to see that the "detailed balancing" relation 

pr^m + l)Ra{m + 1) = pr{m)Re{m) (218) 

holds true because m cannot go negative. When the expression of pr{m) obtained in (12151) 
is introduced in (I218|) . one finds that 

Rejm) _ pr{m + I) _ N - m 

-Ra('TT' + 1) pr{m) 771+1 

Re must be proportional to the number n = N — m atoms in the upper state while Ra must 
be proportional to the number N — n = m of atoms in the lower state. We therefore set 
Re{m) = (^N — m) f {m) , Ra{m) = mg{m), where f{m) and g{m) are two functions to be 
determined. Substituting in (|219|) we find that 

/(m)=5(m + l). (220) 

Because we assume that atoms emit or absorb a single light quantum at a time ("1-photon" 
process) the two functions /(m) and g{m) must be of the linear form f{m) = am + b and 
g{m) = cm + d, where a, b, c, d are constants. But Ra is required to vanish for m = since. 
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otherwise, m could go negative, and thus d = 0. Substituting into (|220p . we find the relation 
am+b = c(m+l) which must hold for any m- value. Therefore, a = b = c. Setting for brevity 
a = b = c = 1 amounts to fixing up a time scale. Then Rein, m) = n{m + l) , Ra{n, m) = 
(A^ — n)m. We note here a lack of symmetry between the rate of stimulated emission 
(proportional to m + 1) and the rate of stimulated absorption (proportional to m). Since 
according to the Schrodinger equation the two processes should be similar, we are led to 
define the field energy as m + 1/2, and to express R^ and Ra in terms of the field energy 
at jump time, Ej^mp time-, defined as the arithmetic average of the field energy just before 
and just after the jump. If we do so, we finally obtain, setting n = and N — n = Ua 

Rei'^ejl^) — ^e Eji^^np time (^^1) 
RaiP'ai'^) — Ejn^np timei (222) 

and the symmetry is indeed restored. 

Let us now restrict our attention to the steady-state regime and large values of A^. Since 
m is large, it may be viewed as a continuous function of time with a well-defined time- 
derivative. Because the standard deviation y'N/A of m is much smaller than the average 
value, the so-called "weak-noise approximation" is permissible. Within that approximation, 
the average value of any smooth function /(n, m) may be taken as approximately equal to 
/((n),(m)). 

The evolution in time of a particular realization m{t) of the process obeys the classical 
Langevin equation 

^ = .-a. ,223) 

where 

= Re{m) + e{t) ^ = Raim) +a{t). (224) 

In these expressions, e(t) and a{t) represent uncorrelated white-noise processes whose spec- 
tral densities are set equal to l3Re{{m)) and f3Ra{{m)), respectively, where /5 is a constant 
to be determined. 

Let us show that the variance of m obtained from the above Langevin equation coincides 
with the result obtained directly from Statistical Mechanics only if /? = 1. Without the 
noise sources, the evolution of m in (I223P would be deterministic, with a time-derivative 
equals to the drift term Re{m) — Ra{m). If the expressions (|220p are employed, the Langevin 
equation (|223p reads 

dm „ 2 

= JSm — 2m + e — a 

dt 

^e-a = /3{Re + Ra) = /3N (m) = /3A^V2, (225) 
where the approximation ^ 1 has been made. 
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Let m(t) be expressed as the sum of its average value (m) plus a small deviation Am(t), 
and Nm — 2m? in (|225l) be expanded to first order. A Fourier transformation of Am(t) 
with respect to time amounts to replacing d/dt by jVL. The Langevin equation now reads 

jJlAm = -iVAm + e - a y^-a = f3N^/2, (226) 

where m has been replaced by its average value N/2. 

Since the spectral density of some random function of time z{t) = ax{t), where a = 
a' + ja" is a complex number and x{t) a stationary process, reads J^zi^) = see 
Section [3?6l one finds from (|226p that the spectral density of the Am{t) process is 

yAmm = (227) 

The variance of m is the integral of J^Am(f^) over frequency {Q/2Tr) from minus to plus 
infinity, that is var(m) = PN/A. There is agreement with the previous result derived from 
the basic Satistical Mechanics rule only if /3 = 1. It follows that the spectral density of 
fluctuation rates such as r{t) must be equal to the average rates, say R. 



7.6 Conclusion of the present section 

We have offered independent (but, admittedly, partly heuristic) methods of showing that 
the spectral density of the real part of the fundamental current source associated with a 
conductance G at frequency cOq is equal to hujo \G\, or, equivalently, that the spectral density 
of the real part of the fluctuation-rate equals the absolute value of the average rate, |i2|, 
when only one of the two states is populated. 

In general the admittance Y{uj,n, R,Tm) of a circuit depends on frequency w, on the 
number n of atoms or electrons in the conduction band or upper state, on the emitted rate 
R of electromagnetic radiation, and possibly on temperature Tm- In the above discussion 
we have implicitly assumed that Y does not depend on R, that is, no spectral-hole burning 
occurs. In laser diodes, the number n of electrons in the conduction band is a monotonic 
function of the static electric potential U applied, provided some unessential internal resis- 
tance be neglected. It follows that Y does not depend on n if the applied static potential U 
is held constant. There are other circumstances where n is approximately constant. There 
may exist an equilibrium between the electron injection rate and the spontaneous decay 
rate. As long as the departure from this equilibrium condition remains small, we may as- 
sume that n is approximately constant. If this is the case, the admittance y is a function 
of uj only, and the system is linear. We have assumed that the optical frequency u remains 
close to some average oscillation frequency uJq- 
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8 Conclusion 



We have shown that the photo-electron spectrum originating from a detector submitted to 
non-fluctuating-pump laser light may be understood in semi-classical terms. We did not 
employ the Quantum Optics concept that the optical field should be treated as an operator 
acting on the state of the light, nor the concept that light consists of point particles called 
"photons". The expression "photon rate" was used only as another name for electromag- 
netic power divided by HuJo, where denotes the average laser frequency. Photon rates are 
written in the form R + AR{t) + r{t), where R denotes the average rate, AR{t) is propor- 
tional to the fundamental noise sources in the linearized regime, and the spectral density 
of the fundamental noise source r{t) is equal to R. This conclusion has been reached from 
different approaches, essentially requiring agreement with Classical Statistical Mechanics 
formulas. We restricted ourselves to the non-relativistic (c — > oo), stationary, linear or 
linearized, approximations. Results obtained from the present theory for some particular 
configurations were listed in the first version of this paper. 

Any laser-detector system is viewed as a set of capacitances, inductances, and con- 
ductances connected to one-another, with driving static potentials or currents. Conduc- 
tances were modeled by a single electron located between two conducting plates, see Fig. 
[5] in (a). The electron is submitted to a static potential source Ue, which determines 
the states, 1 (lower), and 2 (upper) the electron may reside in, with an energy difference 
E2 — El = hujo ~ eU, where U is slightly smaller than Ue- The applied (complex) potential 
source V "sees" a conductance G, which under some conditions does not depend on V. The 
laser-detector ensemble is similar to the klystron represented in Fig. [5] in (b). To empha- 
size the symmetry between stimulated emission and stimulated absorption, the absorber of 
radiation schematically shown on the right of that figure should be replaced by a device 
similar to the emitter of radiation, but with a static potential source Ua slightly smaller 
that U. 

The optical-potential source V could be materialized by a tuned circuit with a very large 
capacitance and a very small inductance resonating at the frequency uJq- Then V is almost 
independent of the induced current. In a real device, the capacitance of the tuned circuit 
is not arbitrarily large, that is, the "photon lifetime" is finite. The optical potential V 
across the tuned circuit then depends on the induced current and it is necessary to evaluate 
the complete circuit response to the noise sources, whose statistics is known. The case of 
"quiet" lasers differs from the set-up just described in that the static potential source Ue is 
being replaced by a static current source Jg. This current source could be materialized by 
a very large inductance with a current Jg flowing through it. In the very large inductance 
limit the delivered current is nearly constant, provided the experiment does not last too 
long. Here, by "quiet laser" we mean that the power collected by the static potential Ua 
does not fluctuate much in the course of time. More precisely, the spectral density of the 
collected power vanishes at small Fourier frequencies. This is also what happens for the 
"grand-mother pendulum", an existing purely classical mechanical system, described at the 
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beginning of this paper introduction. 

The over-all effect of the laser action is to slowly transfer energy from one static source 
to the other. Whether the energy flow occurs from left to right, or the converse, in the 
configurations just described, depends on the small difference existing between Ue and Ua- 
There is therefore a slight irreversible loss of energy, perhaps carried away from the structure 
by acoustical waves. This loss may be as small as one wishes, as long as one does not insist 
on having a fast power transfer. This is also what happens in Carnot heat engines, which 
are reversible in the small-power transfer limit. The main features of the (admittedly highly 
idealized) configurations described in this paper may be understood with a single electron 
in the emitting part and a single electron in the absorbing part, as said above. In reality, the 
two elements involve a large number of electrons. We generally neglect the direct electron- 
electron Coulomb interaction (space-charge effects), but need take into account the Pauli 
exclusion principle and the material temperature (Fermi-Dirac statistics). 

In Part II we intend to consider applications to the linear regime. In that regime, the 
circuit considered involves only constant (positive and negative) conductances G, besides 
capacitances and inductances. That is, the conductances are independent of the potential 
applied to them and do not depend on any other parameter. It suffices in principle to 
evaluate the response of a linear system to current sources C{t) = C'{t) + iC"{t) of known 
spectral density. We suppose that the response spectrum is sharply peaked near some 
frequency lOq, so that the spectral density density of the current sources may be written as 
^C' = "^C" = \G\ , '^c'C" = 0. The linear regime is usually applicable when laser- 
diodes are driven by constant electrical potentials, or when the oscillating or amplified 
signals are small. We also consider the noise properties of linear attenuators and amplifiers. 

In Part III of this paper we intend to consider the linearized regime, and in particular 
the relative noise of idealized laser diodes at high power levels. At such high powers, the 
time derivative of the number n of electrons in the conduction band may be neglected (i.e., 
dn/dt = 0). It follows from that assumption that the net rate of photons entering into the 
cavity is equal to the pumping rate. But the detected rate is different because photons may 
be stored for a while in the cavity before exiting. In that limiting situation no relaxation- 
oscillation appear. Spontaneous emission and gain compression are neglected. The more 
general theory will be given later on. Reference to a circuit simulating the lasers Fourier- 
frequency response is cited. We consider the effect of gain compression, and electrical 
feedback. Frequency noise is considered. 

Part IV will discuss the principles of Quantum Statistical Mechanics, beginning with 
the fundamental concepts introduced by Carnot in 1824. 

Part V would provide an evaluation of parameters such as the optical gain, that enter in 
laser-diode theories, essentially on the basis of the Fermi-Dirac distribution (grand-canonical 
ensemble) . 

Part VI would discuss the properties of classical beams and optical fibers. 
Part VII would consider random distortions of optical waveguides. 
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